ДУ методом Эйлера, Рунге-Кутта (Паскаль)

[Gebbellz]

yx'+x=4x exp3 (в 3 степени) + 3y exp 2(во 2 степени).
отрезок интегрирования - [2.3]
y(2)=1;
точное решение - x=y exp2+y exp2

Помогите с решением задачи методом эйлера или рунге-кутта

[Gebbellz]

upppppppppppp

[DoDge_VipeR]

метод эйлера)
замените только функцию и значения переменных

Uses Crt;
Var
 xn,xk,yn,h,x,y:real;
 i:integer;

Function f(x,y:real):real;
begin
 f:=2*x*x+2*y;
end;

Begin
 ClrScr;
 Writeln(' dy/dx=2x^2+2y ');
 xn:=0; yn:=1; xk:=1; h:=0.1;
 Writeln('--------------------');
 Writeln('| i|   x  |   y   |');
 Writeln('--------------------');
 x:=xn; y:=yn; i:=1;
 Writeln('|', i:2, ' |', x:5:2, ' |', y:7:4, ' |');
 repeat
   y:=y+h*f(x,y);
   Writeln('|', i:2, ' |', x:5:2, ' |', y:7:4, ' |');
   x:=x+h;
   i:=i+1;
 until x>xk;
 Writeln('--------------------');
 Readln;
End.

[DoDge_VipeR]

..................

[Gebbellz]

xm....а куда вписать точность решения x=y exp2+y exp2 ?

[Gebbellz]

модифицированный метод эйлера - где-то графическая ошибка - помогите найти(((

var
a,b,c,h,y2,y3,y22:real;
i:integer;
x,y,yt:array [1...100] of real;
function y1(xx,yy:real):real;
begin
y1:=4*yy*yy+3*yy*yy-xx/yy;end;
function y11(yy:real):real;
begin
y11:=yy*yy*yy+yy*yy;end;
begin
clrscr;
x[1]:=2;y[1]:=1;
write('a=');readln(a);
write('b=');readln(b);
write('h=');readln(h);
for i:=1 to n+1 do begin
x[i+1]:=x[i]+h;
y[i+1]:=y[i]+y1(x[i],y[i]);
y22:=(y1(x[i+1],y[i+1])-y1(x[i],y[i]))/h;
y[i+1]:=y[i]+y1(x[i],y[i])*h+y22*h*h/2;
yt[i]:=y11(x[i]);
write(x[i]:2:1);
write(y[i]:6:4);
write(yt[i]:6:4);
end;readkey;end.

[Gebbellz]

up..........