График к алгоритму Рунге-Кутта для системы дифференциальных уравнений

program modelirovanie_mayatnika;
uses crt, GraphABC;
type
mas = array [1..2] of real;
matrix = array [1..4,1..4] of real;
var
x, h, b:real;
i, j, n, m:integer;
Y:mas;
F:mas;
procedure PRAV(m:integer; var x:real; var Y, F:mas);//podprogramma PRAV
begin
F[1]:=Y[2];
F[2]:=Y[2]-((9.8/0.10)*sin(Y[1]));
end;
procedure RGK(m:integer; h:real; var x:real; var Y:mas);// podprogramma runge-kutta
var
y1:mas;
F:mas;
T:matrix;
j: integer;
begin
PRAV(m,x,Y,F);
for j:=1 to m do
begin
T[1,j]:=h*F[j];
y1[j]:=Y[j]+T[1,j]/2;
end;
x:=x+h/2;
PRAV(m,x,y1,F);
for j:=1 to m do
begin
T[2,j]:=h*F[j];
y1[j]:=Y[j]+T[2,j]/2;
end;
PRAV(m,x,y1,F);
for j:=1 to m do
begin
T[3,j]:=h*F[j];
y1[j]:=Y[j]+T[3,j]/2;
end;
x:=x+h/2;
PRAV(m,x,y1,F);
for j:=1 to m do
begin
T[4,j]:=h*F[j];
y1[j]:=Y[j]+(T[1,j]+2*T[2,j]+2*T[3,j]+T[4,j])/6;
end;
Y:=y1;
end;
//osnovnaya programma
begin
writeln('Введите начальные значения:');
write('Порядок системы M:');
readln(m);
write('Введите значение X:');
readln(x);
writeln('Введите значение функции и производные:');
write('y='); readln(y[1]);
write('y`='); readln(y[2]);
write('Введите h:');
readln(h);
write('Введите b:');
readln(b);
writeln;
writeln('Решение:');
writeln(' i    ','x        ','y          ','y`');
n:=round((b-x)/h);
for i:=1 to n do
begin
        RGK(m,h,x,y);
        if i=1 then
        writeln('  ',i, x:6:2,'   ',y[1]:6:5,'   ',y[2]:6:5,'  ');
        if(i mod 10)=0 then begin
        if(i<100) and(i>=10)then write(' ')
        else if(i<10)then write('  ');
        write(i, x:6:2,'  ');
        if(y[1]<0) then write(y[1]:6:5,'  ')
        else write(' ',y[1]:6:5,'  ');
        if(y[2]<0) then write(y[2]:6:5,'  ')
        else write(' ',y[2]:6:5);
        writeln;
        end;
end;
readln;
end.