Контрольная работа: Расчет сверхзвукового обтекания заостренных тел вращения

begin

sm:=0;

pr:=0;

for k := 1 to i do

begin

pr:=(x[i+1]-e[k-1])/(a1*r[i+1]);

sm:=sm+ m[k]*( pr* sqrt(pr*pr - 1)+ln( pr + sqrt(pr-1) ))

end;

pr:=(x[i+1]-e[i])/(a1*r[i+1]);

m[i+1]:=(1-sm)/( pr* sqrt(pr*pr-1)+ln( pr + sqrt(pr-1) ));

end;

for i:= 1 to n-1 do

begin

sm:=0;

for k:=1 to i do

begin

pr:=(x[i]-e[k-1])/(a1*r[i]);

sm:= sm+ m[i]*sqrt( pr*pr-1)

end;

p[i]:=(-4/a1)*sm;

end;

pr:=(x[1])/(a1*r[1]);

m1[1]:=((x[1]-Xm)/Xk)/(pr*sqrt(pr*pr-1)+arccosh(pr));

for i := 1 to n-1 do

begin

sm:=0;

for k := 1 to i do

begin

pr:=(x[i+1]-e[k-1])/(a1*r[i+1]);