qbasic代码9

machfunc:
IF mg > Iimmach THEN GOTO endloop2
IF (2!-bdum*mg^2)/( 2!-bdum*Mfix^2) < 0 THEN STOP
xdum = mg*(1!+adum*mg^2 )^(-1!*cdum)*( (2!-bdum*mg^2)/(2!-bdum*Mfix^2))^ddum
xdum =xdum-edum
IF jj=0 GOTO it3cont1
IF jj=1 GOTO it3cont2
IF jj>1 AND jjj=0 GOTO it1cont1
IF jj>1 AND jjj= 1 GOTO it1cont2
endloop2:
PRINT" Solution becoming asymptotic i"
PRINT" Upstream Mach No.(M1 ) and downstream search interval [a.b] for M2"
PRINT USING"########.# ";Mfix,m1 ,m2
PRINT"Continue ? M2 and all subsequent math numbers will be set "
INPUT"at INPUT value. (negative number = NO )";mg
IF mg<0 THEN STOP
endloop3:
M2 = mg
END SUB