a3= –0.00157565; a4=0.00916281; a5= –0.02057706;
a6=0.02635537; a7= –0.01647633; a8=0.00392377;
b0=0.39894228; b1= –0.03988024; b2= –0.00362018;
b3=0.00163801; b4= –0.01031555; b5=0.02282967;
b6= –0.02895312; b7=0.01787654; b8= –0.00420059;
Var x, m, mh, q0p, Sih, Coh, Tgh, Efp, TTp, tp, tt0, Chh, Dsp, P1,
P2, Efpp, Tpp, q0pp, tttp, q0tp, Eftp, ttpf, m2h, fi0, fi1, kfb, y, t1, t2, t3,
t4, t5, t6, t7, t8, t10, t12,fb:real;
Function IB0(x:real):real; Begin if x<3.75 then begin t1:=x/3.75;
t2:=sqr(t1);t4:=sqr(t2); t6:=t2*t6; t8:=t4*t4; t10:=t2*t8; t12:=t6*t6;
IB0:=1+a01*t2+a02*t4+a03*t6+a04*t8+a05*t10+a06*t12
end else begin t1:=x/3.75; t2:=sqr(t1);
t3:=t1*t2; t4:=t2*t2; t5:=t1*t4; t6:=t1*t5; t7:=t1*t6; t8:=t1*t7;
fb:=a0+a1/t1+a2/t2+a3/t3+a4/t4+a5/t5+a6/t6+a7/t7+a8/t8;
IB0:=fb*exp(x)/sqrt(x);end end;
Function IB1(x:real):real; Begin if x<3.75 then begin t1:=x/3.75;
t2:=sqr(t1);t4:=sqr(t2); t6:=t2*t6; t8:=t4*t4; t10:=t2*t8; t12:=t6*t6;
IB1:=(0.5+b01*t2+b02*t4+b03*t6+b04*t8+b05*t10+b06*t12)*x
end else begin t1:=x/3.75; t2:=sqr(t1);
t3:=t1*t2; t4:=t2*t2; t5:=t1*t4; t6:=t1*t5; t7:=t1*t6; t8:=t1*t7;
fb:=b0+b1/t1+b2/t2+b3/t3+b4/t4+b5/t5+b6/t6+b7/t7+b8/t8;
IB1:=fb*exp(x)/sqrt(x) end;end;
Procedure GipF(Var w:real; SGP, CGP, TGP:Real); Begin
y:=exp(w); Sih:=(y–1/y)/2; Coh:=(y+1/y)/2; Tgh:=Sih/Coh; end;
Begin ClrScr; tt0:=ts–tc; m:=sqrt(2*aLf/(Lm*d0)); mh:=m*h;
GipF(mh,Sih,Coh,Tgh); Chh:=Coh; Efp:=Tgh/mh;
q0p:=Lm*d0*tt0*m*Tgh;
mh:=0; GipF(mh,Sih,Coh,Tgh); ttp:= tt0*Coh/Chh; tp:=ttp+tc;
Writeln(' Efp=',Efp:5:3,' tp=',tp:4:1,' q0p=',q0p:5:2);
Dsp:=sqrt(1+sqr(2*m*h)); P1:=(Dsp–1)/2; P2:=(Dsp+1)/2;
q0pp:=Lm*d0*tt0*P1/h; Efpp:=1/P2; tpp:=tt0
Writeln('Efpp=',Efpp:5:3,' q0pp=',q0pp:5:2,' tpp=',tpp:3:1);
mh:=m*h; m2h:=2*mh; fi0:=IB0(m2h); fi1:=IB1(m2h); kfb:=fi1/fi0;
q0tp:= 2*Lm*tt0*kfb/m; Eftp:=kfb/mh; m2h:=0; fi0:= Ib0(m2h);
tttp:=tt0*fi0/fi1; ttpf:= tttp+tc;
Writeln('Eftp=',Eftp:5:3,'q0tp=',q0tp:5:2,'ttpf=',ttpf:3:1); end.
Файл результатов примера 5.5:
m = 17.150; Efp = 0.547; tp = 67.4; q0p = 273.25;
Efpp = 0.437; q0pp = 218.69; tpp = 50.0;
mh = 1.715; Eftp =0.515; q0tp = 175.20; ttpf = 106.6.