Лабораторные работы по численным методам / Численные методы / Численные методы / лабораторка 7 / рутте-кутта

//Uses crt;
Const
x0=-1;
y0=1.5;
b=0;
eps=0.5;

Var
x1,y1,x2,y2,y2a,y3,h,p:real;
i,N:integer;

Function f(x,y:real):real;
begin
f:=x*x+2*x+y/(x+2);
end;

Begin
y3:=y0;
h:=0.5;
Repeat
//ClrScr;
x1:=x0;
y1:=y0;
h:=h/2;
N:=Round((b-x0)/h);
for i:=1 to N do
begin
x2:=x1+h;
y2a:=y1+h*f(x1,y1);
y2:=y1+h*(f(x1,y1)+f(x2,y2a))/2;
writeln('x= ',x2:6:3,' y= ',y2:6:3,' h= ',h:0:5);
x1:=x2;
y1:=y2;
end;
p:=Abs(y2-y3);
y3:=y2;
Until p<eps;
Readln
End.