Курсовая работа: Метод Стрілянини

ybeg, yend, t: Extended;

alpha1, alpha2, beta1, beta2 : extended;

b:array of Extended;

a:array of array of Extended;

y,y1,y2:array of Extended;

x:array of Extended;

M0,M1,M2,err:Extended;

k,k0,k1,k2,k3,l,l0,l1,l2,l3,zn,xn:extended;

function f(x,y,dy:extended):extended;

begin

result:=100*y;

end;

begin

try

h:=StrToFloat(LabeledEdit1.Text);

n:=round(1/h)+1;

ybeg:=1; yend:=1;

SetLength(y,n);

SetLength(y1,n);

SetLength(y2,n);

y1[0]:=ybeg;

y2[0]:=ybeg;

alpha1:=yend;

alpha2:=yend;

zn:=alpha1; xn:=0;

for i:=1 to n-1 do begin

l0:=h*f(xn,y1[i-1],zn);

k0:=h*zn;

l1:=h*f(xn+h/2,y1[i-1]+k0/2,zn+l0/2);