qbasic代码7

poode:
p0=7.8617#-2644.1#*(1/ts)-173690#*(1/ts)^2
p0=132.89#*10^p0 ' [N/m^2]
p0prime=LOG(10#)*(2644.1#/ts^2+2#*173690#/ts^3)*p0
GOTO tcont2
tcont1:
ts=tsnew
END SUB
SUB sonic(temp,latheat,qsonic) STATIC
SHARED rho0,k,rg,pi,di
qovera=rho0*SQR(k*rg*temp)*latheat/SQR(2*(k+1))
ac=pi*di^2/4
qsonic=ac*qovera
END SUB
SUB mach2(j,dm,M1 ,M2) STATIC
SHARED k,pi,mdot(),muu(),adum,wdum,delz,dout,epps,lmach,limmach,f
dm=dm/delz : Mfix=M1
Rey=4! *mdot(j)/muu(j)/pi/dout
PRINT "Re No.",Rey
IF Rey<10^5 THEN 302
301 : f=.0791/Rey^(1/4) :GOTO 303
302 : f=SQR((3.6*LOG(Rey/7))^-1)
303 : fstar=f*Rey
bdum=fstar*muu(j)*pi*k/dm
cdum = (1!+adum)/(2!*adum + bdum)
ddum = (1!-bdum)/(2!*adum + bdum)
edum=wdum*M1
edum=edum*(1!+adum*M1^2)^(-1!*cdum)