решение задачи коши методом рунге-кутта 4 порядка в паскале

program rungekut 4;
uses crt;
const NMax=2;
type
Mas=array[1..NMax] of real;
var
Y, Yo: mas;
X, Xo, Xn, h: real;
i, j, M: integer;
ch: char;
procedure p(x: real; y: mas; var f: mas);
begin
f[1]:= y[1]*y[1]/x*x+8*y[1]/x+4;
end;
 
procedure Rgk(x, h: real; n: integer; var y: mas);
var
k1, k2, k3, k4, y1: mas;
i: integer;
begin
p(x, y, k1);
for i:=1 to n  do
  y1[i]:=y[i]+h*k1[1]*0.5;
p(x+h*0.5, y1, k2);
for i:=1 to n do
  y1[i]:=y[i]+ h*k2[i]*0.5;
p(x+h*0.5, y1, k3);
for i:=1 to n do
  y1[i]:=y[i]+h*k3[i];
  p(x+h*0.5, y1, k4);
  for i:=1 to n do
  y[i]:=y[i]+h*0.166657*(k1[i]+2*k2[i]+2*k3[i]+k4[i]);
  end;
  
  begin
    ClrScr;
    writeln( 'Ïðèáëèæåííîå ðåøåíèå çàäà÷è Êîøè');
    writeln( '________________________________');
    writeln( 'Ìåòîä Ðóíãå-Êóòòà 4 ïîðÿäêà');
    writeln( '_______________________________');
    writeln;
    write('Ââåäèòå íà÷àëüíîå çíà÷åíèå Xo: '); readln(Xo);
    write('Ââåäèòå êîíåö èíòåðâàëà Xn: '); readln(Xn);
    write('Ââåäèòå øàã h: '); readln(h);
    write('Ââåäèòå íà÷àëüíûå óñëîâèÿ: ');
    for i:=1 to NMax do
    begin
      write(' Yo[',i,']='); readln(Yo[i]);
      end;
      M:=Trunc((Xn-Xo)/h);
      X:=Xo;
      for i:=1 to NMax do Y[i]:=Yo[i];
      writeln;
      writeln('Ðåçóëüòàòû: ');
      for i:=1 to M do
      begin
        Rgk(X, h, NMax, Y);
        X:=X+h;
        write('X=', X:10:7);
        for j:=1 to NMax do write('Y[',j,'] =', Y[j]:10:7);
        writeln;
        end;
        writeln('Íàæìèòå ëþáóþ êëàâèøó...');
        ch:=readkey;
        end.