выч методы не разобраны! / Лаба№3 / Мишины / MILN

Program Miln;
Uses Crt;
Const
      y0=0.5; x0=1;
      h=0.1 ; a=0;
      b=1   ; n=10;

var
     m     : text;
     i     : integer;
     x,y   : array [0..10] of real;
Function f(xx,yy:real):real;
  begin
    f:=2*xx+yy
  end;
begin
  Assign(m,'Miln.inf');
  ReWrite(m);
  clrscr;
  writeln(' Miln method  : ');
  writeln(m,' Miln method  : ');

  x[0]:=x0;
  y[0]:=y0;
 {=======================ЌЂ—Ђ‹њЌ›‰ Ћ’ђ…‡ЋЉ==================================}
  for i:=1 to n do   x[i]:=x[i-1]+h;
  for i:=0 to 3 do
	y[i+1]:=y[i]+h*0.5*(f(x[i],y[i])+f(x[i+1],y[i])+h*f(x[i],y[i]));
{========================METOD Њ€‹ЌЂ========================================}

 for i:=4 to n do
   begin
 y[i]:=y[i-4]+4*h/3*(2*f(x[i-3],y[i-3])-f(x[i-2],y[i-2])+2*f(x[i-1],y[i-1]));
 y[i]:=y[i-2]+h/3*(f(x[i-2],y[i-2])+4*f(x[i-1],y[i-1])+f(x[i],y[i]))
   end;
 {=======================‚›‚Ћ„ ЌЂ ќЉђЂЌ € ‚ ”Ђ‰‹==========================}
  writeln('y"(x,y)=2*x+y');
  writeln('h=',h:1:1,'  x0=',x0,'  y0=',y0:1:2);
  writeln;
  writeln(m,'y"(x,y)=2*x+y');
  writeln(m,'h=',h:1:1,'  x0=',x0,'  y0=',y0:1:2);
  writeln(m,'           ');
  for i:=0 to n do
    begin
      writeln('y(',x[i]:1:1,')=',y[i]:1:10);
      writeln(m,'y(',x[i]:1:1,')=',y[i]:1:10)
    end;
  close(m);
  readln
 end.