Решение задачи Коши для обыкновенного дифференциального уравнения методом Эйлера

{Решение задачи Коши методом Эйлера}
program Eyler;
uses crt; {модуль управления экраном}
var i,n:integer;
    a,b,h,x,x0,y,y0,eps,emax:real;
 
function f(x,y:real):real; {функция уравнения}
begin
 f:=x*x+y*y;
end;
 
(*function ft(x:real):real; {функция точная}
begin
 ft:=2/x;
end;*)
 
begin
 clrscr; {очистить экран}
 writeln ('Решение задачи Коши методом Эйлера');
 writeln ('Уравнение dy/dx=x^2+y^2)');
 write ('Введите интервал поиска решения [A,B]: ');
 read (a,b);
 write ('Введите число шагов N: ');
 readln (n);
 h:=(b-a)/n;
 x0:=a;
 write ('Введите начальное значение Y(0): ');
 read(y0);
 y:=y0;
 x:=x0;
 emax:=0; {макс. отклонение}
 writeln ('X':19,'Y':19);
 for i:=1 to n+1 do begin
  y:=y+h*f(x,y); {делаем шаг метода}
  (*eps:=abs(y-ft(x));
  if eps>emax then emax:=eps; {оцениваем макс. погрешность}*)
  writeln (x:19:8,y:19:8);
  x:=x+h;
 end;
 {writeln ('Pmax=',emax:19:8);}
 readln;
end.

{Решение задачи Коши методом Эйлера}
program Eyler;
uses crt; {модуль управления экраном}
var i,n:integer;
    a,b,h,x,x0,y,y0,eps,emax:real;
 
function f(x,y:real):real; {функция уравнения}
begin
 f:=-y+exp(-x);
end;
 
function ft(x:real):real; {функция точная}
begin
 ft:=x*exp(-x);
end;
 
begin
 clrscr; {очистить экран}
 writeln ('Решение задачи Коши методом Эйлера');
 writeln ('Уравнение dy/dx=y/x-4/sqr(x)');
 write ('Введите интервал поиска решения [A,B]: ');
 read (a,b);
 write ('Введите число шагов N: ');
 readln (n);
 h:=(b-a)/n;
 x0:=a;
 write ('Введите начальное значение Y(0): ');
 read(y0);
 y:=y0;
 x:=x0;
 emax:=0; {макс. отклонение}
 writeln ('X':19,'Y':19,'Eps':19);
 for i:=1 to n+1 do begin
  y:=y+h*f(x,y); {делаем шаг метода}
  eps:=abs(y-ft(x));
  if eps>emax then emax:=eps; {оцениваем макс. погрешность}
  writeln (x:19:8,y:19:8,eps:19:8);
  x:=x+h;
 end;
 writeln ('Pmax=',emax:19:8);
 reset (input); readln;
end.