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

Series4: TLineSeries;

Series5: TLineSeries;

procedure FormCreate(Sender: TObject);

procedure Button1Click(Sender: TObject);

type

vector = array[0..101]of real;

var

Form1: TForm1;

n:integer;

m,p,m1,p1:Vector;

x,r,e:vector;

implementation

procedure DoCount(const Minf,Xk,Xmid,Xm,Rmid:real);

var a1,pr,sm:real;

k,i:Integer;

begin

a1:= sqrt( sqr(Minf)-1);

x[0]:=0;

for i := 1 to n do

x[i]:=x[i-1]+Xk/n;

for i:=0 to n do

if x[i]<=Xmid then r[i]:= (Rmid/Xmid)*(2-x[i]/Xmid)*x[i]

else r[i]:= Rmid ;

for i := 1 to n do

begin

e[i]:=x[i]-a1*r[i];

end;

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

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