Решает ли эта программа данную задачу

program Eiler;
uses crt,graphABC;
const 
h=0.001;{Точность}
var 
  i,j,p:integer;
  n,q:real;
  s:string;   { Строка } 
  r1,r2,r3,r4,yi,xi,dy,xm,ym:real;
  {Одномерные массивы решения данного дифф. уравнения}
  a:array[0..2000] of real;
  b:array[0..2000] of real;
  a1:array[0..2000] of real;
  b1:array[0..2000] of real;
{Данное дифф. уравнение} 
function f(x,y:real):real;
  begin 
    f:=sqr(x)+sqr(y); 
  end;
begin 
clrscr;
 {Начальное значение интервала}
xi:=1; 
{Начальное условие}
yi:=-1; 
{Начальные номера массивов}
i:=0; j:=0; 
 {Начальное значение интервала}
xm:=1; 
{Начальное условие}
ym:=-1;
q:=-2.5;
{Метод Рунге-Кутта}
  writeln('Метод Рунге-Кутта');
  writeln('№=',i, '     | ','x=',xi:5:3, ' | ','y=',yi:5:3);
  repeat
    r1:=h*f(xi,yi);
    r2:=h*f(xi+h/2,yi+r1/2);
    r3:=h*f(xi+h/2,yi+r2/2);
    r4:=h*f(xi+h,yi+r3);
    dy:=(r1+2*r2+2*r3+r4)/6;
    yi:=yi+dy;
    xi:=xi+h;
    a[i]:=xi;{Заполнение массива элементом x}
    b[i]:=yi;{Заполнение массива элементом y}
    i:=i+1;{Счётчик}
    if i mod 200=0 then
    { Вывод результатов }
    writeln('№=',i, ' | ','x=',xi:5:3, ' | ','y=',yi:5:3);
  until xi>2;
readln;
{Метод Эйлера}
writeln('Метод Эйлера:');
writeln('№=',j, '      | ','x=',xm:5:3, ' | ','y=',ym:5:3);
  repeat
    ym:=ym+h*f(xm,ym);
    xm:=xm+h;
    a1[j]:=xm;{Заполнение массива элементом x}
    b1[j]:=ym;{Заполнение массива элементом y}
    j:=j+1;{Счётчик}
    if j mod 200=0 then
    { Вывод результатов }
    writeln('№=',j, '  | ','x=',xm:5:3, ' | ','y=',ym:5:3);
  until xm>2;
{ Построение графика функции }
SetWindowHeight(700);{Высота экрана в пикселях}
  SetWindowWidth(900);{Ширина экрана в пикселях}
  SetPenColor(clBlack);{Цвет линий}
    SetPenWidth(3);{Толщина линий в пикселях}
    {Построение осей графика 0xy}
    line(180,400,900,400);
    line(280,0,280,700);
    line(280,0,280-10,10);
    line(280,0,280+10,10);
    line(900,400,900-10,400-10);
    line(900,400,900-10,400+10);
{Построение шкал графика}
for p:=0 to 7 do
   
    line(180+100*p,395,180+100*p,405);
   
for p:=0 to 4 do
   
    line(275,100+100*p,285,100+100*p);
  
{Вывод на экран текста,заключённого в ' '}
TextOut(885,370,'х');
TextOut(295,0,'у');
{Вывод на экран цифрового значения шкал}
q:=-0.5;
for p:=0 to 7 do
  begin
    str(q:4:1,s);{Преобразует элемент q в строковый s}
    TextOut(166+100*p,410,s);{Вывод элемента s на экран}
    q:=q+0.5;{Счётчик}
  end;
 q:=-1;
for p:=0 to 5 do
  begin
    str(q:4:2,s);{Преобразует элемент q в строковый s}
    TextOut(280,600-100*p,s);{Вывод элемента s на экран}
    q:=q+0.5;{Счётчик}
  end;
  {Построение сетки графика}
for p:=0 to 6 do
  begin
   SetPenColor(clGreen);{Цвет сетки}
   SetPenWidth(1);{Толщина линий сетки в пискелях}
   line(280+100*p,0,280+100*p,700);{Построение линий сетки по вертикали}
  end;
  for p:=0 to 5 do
  begin
   SetPenColor(clGreen);{Цвет сетки}
   SetPenWidth(1);{Толщина линий сетки в пискелях}
   line(190,100+100*p,900,100+100*p);{Построение линий сетки по вертикали}
  end;
  {Построение графика функции}
  while p<=2000 do {Кол-во точек графика}
    begin
    n:=200;{коэф.}
    Setpixel(280+round(a[p]*n),round(400-b[p]*n),clBlue); {Вывод на экран пикселей с указанным цветом 1-го графика}
    Setpixel(round(280+a1[p]*n),round(400-b1[p]*n-10),clRed);{Вывод на экран пикселей с указанным цветом 2-го графика}
    p:=p+1;{Счётчик}
    end;
readln;
end.