qbasic代码

#qbasic#
executor:
INPUT" ",ss$
OPEN "zdat" FOR OUTPUT AS #1
OPEN "mdotdat" FOR OUTPUT AS #2
OPEN "tdat" FOR OUTPUT AS #3
OPEN "pdat" FOR OUTPUT AS #4
OPEN "rhodat" FOR OUTPUT AS #5
OPEN "machdat" FOR OUTPUT AS #6
OPEN "bigxdat" FOR OUTPUT AS #7
OPEN "titlesdat" FOR OUTPUT AS #8
OPEN "miscdat" FOR OUTPUT AS #9
OPEN "t.profile" FOR OUTPUT AS #10
OPEN "tmach" FOR OUTPUT AS #11
solution:
i=0 : jcount=0
INPUT" Enter stagnation temperature, [K]",t(0)
tsat(0) = t(0)
GOSUB satpress
rho0=p0/rg/tsat(i) ' [kg/m^3]
p(i+1)=p0/(1+k)
t(i+1)=2!*t(0)/( k+1!)
i=1
CALL sattemp(p(i),t(0),tsat(i))
GOSUB latheat
hfg(0)=hfg(1)
CALL sonic(t(0),hfg(1),q(1))
tsat(i)=tsat(i)
delq(i)=0
mdot(i)=q(i)/hfg(i)
mach(i)=1! :mach(i-1)=0!
rho(i) =p(i)/rg/tsat(i)
z(0)=0