GN = input('What is the name of the gas?','s')
holdon
IP = input('What is the initial pressure in the tank measured at 270K, in [atm]')
holdon
GT = input('the temperature of interest, in a range of 270 – 480K')
holdon
Nitrogen_P=2.5
Nitrogen_V=12
Nitrogen_T=270
Nitrogen_R=0.0821
Nitrogen_N=(Nitrogen_P*Nitrogen_V)/(Nitrogen_R*Nitrogen_T)
Oxygen_P=4
Oxygen_V=15
Oxygen_T=270
Oxygen_R=0.0821
Oxygen_n=(Oxygen_P*Oxygen_V)/(Oxygen_R*Oxygen_T)
plot (GT, IP,'r--')
holdon
plot (Nitrogen_T,Nitrogen_P,'k--')
holdon
plot (Oxygen_T,Oxygen_P,'b-')
holdon
xlim([250 500])
ylim([0 inf])
gridon
legend('GN','Nitrogen','Oxygen')
xlabel ('Temperature (T) [K]')
ylabel ('Pressure (P) [atm]')
title ('Pressure Change by Temperature')
