методом Рунге-Кутта с точностью 0,01 решить дифференциальное уравнение: k'=f(x,y) на интервале [0,1] с заданным начальным условием k(a)=1.

rogram kurs;
uses crt,graph;
const
xmin=0;
xmax=1;
h=0.01;
n=10;
var
   st:string;
 x0,y0,gd,gm,mk,i,a,b,px,py:integer;
k0,k1,k2,k3,mx,my,z,x,y,dx,dy,yi,xi:real;
 kx,ky:real;
 ymin,ymax:real;
 s:string;
 function f(x,y:real):real;
begin
 f:=((y/(x+1))-y*y);
end;
begin
  xi:=0;
  yi:=1;
  writeln(' k|  xk |  yk');
  for i:=0 to n do
  begin
    k0:=h*f(xi,yi);
    k1:=h*f(xi+h/2,yi+k0/2);
    k2:=h*f(xi+h/2,yi+k1/2);
    k3:=h*f(xi+h,yi+k2);
    dy:=(k0+2*k1+2*k2+k3)/6;
    yi:=yi+dy;
    xi:=xi+h;
     writeln(i:2, '|',xi:5:2, '|',yi:7:4);
     end;
   readkey;
    begin
     gd:=detect;
     initgraph(gd,gm,'egavga.bgi');
     x:=xmin;
     z:=f(x,y);
     repeat
      if y<ymin then ymin:=y;
      if y>ymax then ymax:=y;
      x:=x+h;
      until(x>=1);
     x0:=getmaxx div 10;
     y0:=getmaxy div 2;
     mx:=500/abs(1);
     my:=400/abs(1.5);
       setbkcolor(black);
     setlinestyle(0,0,1);
 
 
 
         setcolor(green);
     line(0,y0,getmaxx,y0);
     line(x0,getmaxy,x0,0);
     setcolor(red);
     for i:=1 to 10 do
      begin
       line(10,y0-round(i*my),630,y0-round(i*my));
       line(10,y0+round(i*my),630,y0+round(i*my));
       line(x0+round(i*mx),0,x0+round(i*mx),getmaxy);
       line(x0-round(i*mx),0,x0-round(i*mx),getmaxy);
       str(i,s);
       outtextxy(x0-10,y0-round(i*my),s);
       outtextxy(x0-15,y0+round(i*my),'-'+s);
      end;
      setcolor(white);
      outtextxy(x0+5,y0+10,'0');
      outtextxy(getmaxx-10,y0-10,'X');
      outtextxy(x0-10,10,'Y');
      setcolor(red);
      x:=xmin;
      dx:=0.01;
        for i:=0 to n do
       begin
         z:=f(x,y);
        px:=x0+round(xi*mx);
        py:=y0-round(yi*my);
        putpixel(px,py,red);
        xi:=xi+dx;
        yi:=yi+dy;
      end;
      setcolor(red);
      setlinestyle(0,0,1);
      outtextxy(230,30,'F=(y/x+1)-y^2');
      readln;
      closegraph;
   end;
   end.