#!/usr/bin/gnuplot # hugonita.plt # Hugoniota v porovnani s adiabatou. # Miroslav Broz (miroslav.broz@email.cz), May 6th 2009 # tj. fce p1/p2 hugoniota(rho2_rho1) = ((gamma+1)-(gamma-1)*rho2_rho1) / ((gamma+1)*rho2_rho1 - (gamma-1)) # adiabata: # p V^gamma = konst. # p rho^-gamma = konst. # p1 rho1^-gamma = p2 rho2^-gamma # p1/p2 = (rho2/rho1)^-gamma adiabata(rho2_rho1) = rho2_rho1**(-gamma) gamma = 1.4 set xl "{/Symbol-Oblique r}_1/{/Symbol-Oblique r}_2" set yl "{/Helvetica-Oblique p}_1/{/Helvetica-Oblique p}_2" offset +1,0 set xr [0:6] y1 = -1 y2 = 20 set yr [y1:y2] set ytics 5 set mytics 5 set zeroaxis set samples 1000 set nokey set label "1/6" at 1./6., graph 1.04 center tmp=4.25 set label "{/Helvetica-Bold hugoniota} " at tmp, hugoniota(1./tmp) right tmp=4.0 set label " {/Helvetica-Bold adiabata}" at tmp, adiabata(1./tmp) left set style line 1 lt 1 lw 5 set style line 2 lt 7 lw 3 set style line 3 lt 0 lw 1 set lmargin 5.9 set rmargin 0.7 set bmargin 3.1 set tmargin 1.3 p adiabata(1./x) ls 2,\ hugoniota(1./x) ls 1,\ "vline.dat" u ($1*1./6.):(y1+(y2-y1)*$2) not w l ls 3 #pa -1 set term post eps enh color solid "Helvetica" 18 set out "hugoniota.eps" set size 0.75,0.8 rep