# FH HN kr   5/2002
#
#
set dummy V    
set grid 
set title "Isothermen des idealen Gases (n=1 mol)" 0.000000,0.000000  ""
set xlabel "Volumen  V /l" 
set xrange [ 0 : 100 ] 
set ylabel "Druck p / bar" 
set yrange [ 0 : 5 ] 
set nokey
#
p1 = 1.
V1 = 22.4
T1 = 273.
#
## Polytrope ZÄ
#
p(V,n)=p1*(V1/V)**n
#
## Isotherme
#
isot(V,T)=p1*V1*T/(T1*V)
#
# Beschriftung ...
#
T0 = 2000.
V0 = 40.
p0 = p1 * (V1/V0) * (T0/T1) 
#
T =   50.; set label 1 "  50 K" at sqrt(T/T0)*V0,sqrt(T/T0)*p0 left
T =  100.; set label 2 " 100 K" at sqrt(T/T0)*V0,sqrt(T/T0)*p0 left
T =  200.; set label 3 " 200 K" at sqrt(T/T0)*V0,sqrt(T/T0)*p0 left
T =  273.; set label 4 " 273 K" at 1.3*sqrt(T/T0)*V0,sqrt(T/T0)*p0/1.3 left
T =  300.; set label 5 " 300 K" at sqrt(T/T0)*V0,sqrt(T/T0)*p0 left
T =  500.; set label 6 " 500 K" at sqrt(T/T0)*V0,sqrt(T/T0)*p0 left
T = 1000.; set label 7 "1000 K" at sqrt(T/T0)*V0,sqrt(T/T0)*p0 left
T = 2000.; set label 8 "2000 K" at sqrt(T/T0)*V0,sqrt(T/T0)*p0 left
#
plot isot(V,50),isot(V,100),isot(V,200),isot(V,273),isot(V,300),isot(V,500),isot(V,1000),isot(V,2000)
#    EOF