выч методы не разобраны! / Лаба№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.