Численные методы на Паскале

[Veina]

Я написала программу на Паскале для формулы yn+1=yn+h*fn
Вот код:

program eiler;

var x,y,a,b,h:real;

function f(x,y:real):real;

begin f:= y*y*(8-20*x); end;

begin

writeln('vvedite y, a, b, h');

readln(y,a,b,h); x:=a;

repeat

writeln(x:0:1,' ',y:0:3);

y:=y+h*f(x,y);

x:=x+h;

until not (x<b);

readln;

end.

А как написать переделать для формулы yn+1=yn+h*fn+1
и для формулы yn+1=yn+(h/2)*(3fn-fn-1)
Подскажите, пожалуйста, кто-нибудь. Заранее спасибо

[Veina]

Ну откликнитесь же уже кто-нибудь....пожалуйста.......

[phomm]

если функция та же самая (которая y*y(8-20*x)) там в коде) то по идее
для второго случая :
вместо y:=y+h*f(x,y); писать y:=y+(h/2)*(3*f(x,y)-t); где t - сохраненное в переменною старое значение функции, для чего надо после расчета дописать t := f(x,y); при чем надо ещё рассмотреть условие того, что на первом шаге это значение неизвестно(ну я не знаю что у вас, может там и 0), конечно , это t надо ещё объявить

для первого случая мне влом раскидывать свой мозг, чтобы въехать в то , как это у вас в функции расчета y[n+1] используется функция от y[n+1] , пока первая моя мысль - вы что-то напутали, ну хотя может и я уже всё забыл (((