# FH HN   kr
set dummy V
set grid
set key Left rev
set nolabel
set noarrow
set nolog
set title "Carnot-Prozess" 0,0
set xlabel "Volumen V  / l" 0,0
set xrange [0 : 25]
set ylabel "Druck p / bar" 0,0
set yrange [0 : 30]
#
Rm=8.3144 * 1000.
#
p(V,kap)=p0*(V0/V)**kap
isot(V,T)=1.E-5*Rm*T/V
kap = 1.4
eps = 2.
V1 = 22.4
T1 = 273.0
p1 = isot(V1,T1)
T3 = 573.0
#
V2= V1/eps
p2= isot(V2,T1)
#
V3= V2*(T1/T3)**(1./(kap-1))
p3= isot(V3,T3)
#
V4= eps*V3
p4= isot(V4,T3)
#
set label 1 "(1)" at V1,p1
set label 2 "(2)" at V2,p2
set label 3 "(3)" at V3,p3
set label 4 "(4)" at V4,p4
#
plot isot(V,T1) tit "273 K",isot(V,T3) tit "573 K",  V0=V2,p0=p2, p(V,1.4) tit "Adiab. 2->3  ", V0=V1,p0=p1, p(V,1.4) tit "Adiab. 4->1  "