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

program rk;

var
x,x0,y,y0,y1,y2,a,b,h,k1,k2,s,r: real;

function f(x,y: real): real;
begin
f:=x*sin(x)+y/x;
end;

begin
//writeln('vvedite y, a, b, h');
x0:=pi/2;
h:=(pi-pi/2)/8;
y0:=1;
repeat
//writeln(x0:1:5);
writeln( y0:1:5);
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
{writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( );
writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); writeln( ); }
k1:=f(x0,y0);
k2:=f(x0+h,y0+k1*h);
r:= (k1+k2)/2;
y1:=y0+h*r;
x0:=x0+h;
y0:=y1;
until not (x0<=pi+0.0001);
writeln(h);
end.