qbasic代码

#qbasic#
PRINT " q(1), mdot(1) , P(1)";
PRINT USING" ##.#####^^^^"; q(1),mdot(1),p(1)
FOR i=1 TO n
IF (tsat(i-1)-tsat(i))/tsat(i) < epps3 THEN delz=2!*delz
IF i=1 THEN z(i)=0
z(i) =z(i-1)+delz
asurf=pi*dout*delz
delq(i)=-1*hbar*(tsat(i)-tsink)*asurf
GOSUB viscos
hfg(i)=hfg (1)
q(i+1)=q(i)+delq(i)
delmdot(i) = -1!*delq(i)/hfg(i)
mdot(i+1) = mdot(i) - delmdot(i)
IF mdot(i+1) < 0! THEN GOTO endroutine
lmach=mach(i-1)
wdum = mdot(i+1)/mdot(i)
IF mach(i)>limmach THEN mach(i+1)=mach(i): GOTO downstre
CALL mach2((i), (delmdot(i)),(mach(i)),mach(i+1))
downstre:
p(i+1) = p(i)* wdum*(mach(i)/mach(i+1))*SQR((1! + adum*mach(i)^2)/(1!+adum*mach (i+1)^2))
CALL sattemp(p(i+1), (tsat(i)),tsat(i+1))
tsat(i+1)=tsat(i+1)-beta
rho(i+1)=p(i+1)/rg/tsat(i+1)