милн (хуйзнат)

Program miln;
Uses crt;
type real=double;
const
  n=10;
  h=0.01;
function F(x, y: real): real;
begin
  F:= (-1*(2*x*y)/(y*y-3*x*x));
end;
Function ff(x: real): real;
begin
  ff:=(y*y*y-y*y);
end;
var
  i: integer;
  k: word;
  x: array[0..n] of real;
  y: array[0..n, 1..5]  of real; 
  yy: real;
  w, a1, a2, b1, b2: real;
begin
  clrscr;
  GOTOXY(20,2);

  WRITELN;
  y[0,1]:=ff(2); y[0,2]:=ff(2); y[0,3]:=ff(2); y[0, 4]:=ff(2);
  x[0]:=(2);
  for k:=1 to n do x[k]:=x[0]+h*k;
  {Метод Эйлера (простой)}
  for k:=0 to 2 do
    y[k+1, 1]:=y[k, 1]+h*F(x[k], y[k, 1]);
  {Метод Эйлера (модифиц)}
  for k:=0 to 2 do begin
    y[k+1, 2]:=y[k, 2]+h*F(x[k], y[k, 2]);
    y[k+1, 2]:=y[k, 2]+h/2*(f(x[k], y[k, 2])+f(x[k+1], y[k+1, 2]));
  end;
  { Метод Эйлера (уточнен)}
  for k:=0 to 2 do begin
    y[k+1, 3]:=y[k, 2]+h*F(x[k], y[k, 2]);
  {Уточняем}
    repeat
      yy:=y[k+1, 3];
      y[k+1, 3]:=y[k, 3]+h/2*(f(x[k], y[k, 3])+f(x[k+1], y[k+1, 3]));
    until y[k+1, 3]-yy<=0.000001
  end;
  {с наилучшей w}
  for k:=0 to 2 do begin
    w:=0.75;
    y[k+1, 4]:=y[k, 4]+h*((1-w)*f(x[k], y[k, 4])+w*f(x[k]+1/(2*w)*h, y[k, 4]+1/(2*w)*h*f(x[k], y[k, 4])));
 end;
  {точное решение}
  for k:=0 to n do
    y[k, 5]:=ff(x[k]);
  {Метод Милна (для всех столбиков)}
  for i:=1 to 4 do {четыре столбика}
    for k:=4 to n do begin
      {Предположительное значение:}
      y[k, i]:=y[k-4, i]+4*h/3*(2*f(x[k-3], y[k-3,  i]) -
                                  f(x[k-2], y[k-2,  i]) +
                                2*f(x[k-1], y[k-1,  i]));
      {Уточняем:}
      y[k, i]:=y[k-2, i] + h/3*(  f(x[k-2], y[k-2,  i]) +
                                4*f(x[k-1], y[k-1,  i]) +
                                  f(x[k],   y[k,  i]));
    end;
  WriteLn('N':3, '¦ ', 'Простой':12, '¦ ', 'Модифицир':12, '¦ ', 'Уточненный':12, '¦ ', 'С наилучшим':12, '¦ ','Точное':10);
  WriteLn(' ':3, '¦ ', 'Эйлер':12, '¦ ', 'Эйлер':12, '¦ ', 'Эйлер':12, '¦ ', 'w':12, '¦ ', 'значение':10);
  Writeln('===+=============+=============+=============+=============+===============');
  for i:=0 to n do begin
    WriteLn(i:3, '¦ ', y[i, 1]:12:7, '¦ ', y[i, 2]:12:7,  '¦ ', y[i, 3]:12:7, '¦ ',
            y[i, 4]:12:7, '¦ ', y[i, 5]:10:7);
        end;
 readln;
end.