%%HP: T(1)A(D)F(.); @ FLOW, isentropic gas flow equation solver @ by Denis Hayner DIR W 1.14226035712E-3 T0 1460 T 466.416666667 P0 100 P 5.28281787717 MN 4.74775881895E-4 V 1159.19848414 G 32.174 K 1.4 R 53.3 FLOW 'W=ƒ(G*K/R)* P0/ƒT0*MN/(1+(K-1)/ 2*MN^2)^((K+1)/(2*( K-1)))' TEMP 'T=T0/(1+(K- 1)/2*MN^2)' PRES 'P=P0/(1+(K- 1)/2*MN^2)^(K/(K-1) )' VEL 'V=ƒ(G*K*R*T0 )*MN' F.EQ { FLOW VEL PRES TEMP } EQ TEMP PPAR { (-6.5,-11.216753359) (6.5,75) MN 0 (0,0) FUNCTION Y } END