Решение дифференциального уравнения методом Эйлера

##
  function f(x: real): real := exp(-x) * (cos(x) + sin(x) + x * sin(x));
  var (y, dy, x0, xk, h) := (1.0, 0.0, 0.0, 1.0, 0.05);
  'x    Метод    Tочное'.Println;
  '     Эйлера   значение'.Println;
  $'{x0:N2} {y:N6} {f(x0):N6}'.Println;
  foreach var x in PartitionPoints(x0 + h, xk, round((xk - x0) / h) - 1) do
    begin
      var t := dy;
      dy += 2 * h * (exp(-x) * cos(x) - t - y);
      y += h * t;
      $'{x:N2} {y:N6} {f(x):N6}'.Println
    end

function f(x: real): real;
begin
  f := exp(-x) * (cos(x) + sin(x) + x * sin(x))
end;
 
const
  y0 = 1.0;
  dy0 = 0.0;
  x0 = 0.0;
  xk = 1.0;
  h = 0.05;
var
  x, y, dy, t: real;
  k: integer;
begin
  y := y0;
  dy := dy0;
  writeln('x    Метод    Tочное');
  writeln('     Эйлера   значение');
  writeln(x0:4:2, y0:9:6, f(x0):9:6);
  for k := 1 to round((xk - x0) / h) do
    begin
      x := x0 + k * h;
      t := dy;
      dy := t + 2 * h * (exp(-x) * cos(x) - t - y);
      y := y + h * t;
      writeln(x:4:2, y:9:6, f(x):9:6);
    end;
  readln
end.

function f(x: real): real;
begin
  f := exp(-x) * (cos(x) + sin(x) + x * sin(x))
end;
 
const
  y0 = 1.0;
  dy0 = 0.0;
  x0 = 0.0;
  xk = 1.0;
  h =  0.05;
var  x, y, dy, t, fx, x2, x4, y1, dy1 : real;
  k: integer;
begin
  y := y0;
  dy := dy0;
  writeln('x    Метод    Tочное    погрешность ');
  writeln('     Эйлера   значение       %');
  writeln(x0:4:2, y0:9:6, f(x0):9:6,x:7:2);
  for k := 1 to round((xk - x0) / h) do
    begin
      x := x0 + k * h;
      t := dy; 
      //{ Модифицированный метод
      x2 := x-h/2.0; //вводим промежуточные точки 
      x4 := x2-h/4.0;
      dy1:= t +  h * (exp(-x4) * cos(x4) - t - (y+t*h/4.0));
      y1 := y + h * (t+dy1)/4.0;
      dy := t + 2 * h * (exp(-x2) * cos(x2) - dy1 - y1);
      y := y + h * (t+dy)/2.0;
      //}
      { обычный метод
     dy := t + 2 * h * (exp(-x) * cos(x) - t - y); 
     y := y + h * t;
     }
      fx:=f(x);
      writeln(x:4:2, y:9:6, fx:9:6,   abs(200*(fx-y)/(fx+y)):11:6);
    end;
 
end.

  Float128 = record
    private a,b,c,d: integer;
    
    //TODO Все методы и операторы
    
  end;