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

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

k2:=h*(zn+l2/2);

l3:=h*f(xn+h,y1[i-1]+k2,zn+l2);

k3:=h*(zn+l3);

k:=(k0+2*k1+2*k2+k3)/6;

l:=(l0+2*l1+2*l2+l3)/6;

y1[i]:=y1[i-1]+k;

zn:=zn+l;

xn:=xn+h;

end;

beta1:=y1[n-1]; beta2:=beta1;

while (beta1=beta2) do begin

alpha2:=alpha2-h;

zn:=alpha2; xn:=0;

for i:=1 to n-1 do begin

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

k0:=h*zn;

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

k1:=h*(zn+l1/2);

l2:=h*f(xn+h/2,y2[i-1]+k1/2,zn+l1/2);

k2:=h*(zn+l2/2);

l3:=h*f(xn+h,y2[i-1]+k2,zn+l2);

k3:=h*(zn+l3);

k:=(k0+2*k1+2*k2+k3)/6;

l:=(l0+2*l1+2*l2+l3)/6;

y2[i]:=y2[i-1]+k;

zn:=zn+l;

xn:=xn+h;

end;