Эксперименты лаба10,11(2курс)

#вычисление приближенных значений в первых трех точках методом Эйлера

y[2]:=y[1]+h*f(x[1],y[1]);

y[3]:=y[2]+h*f(x[2],y[2]); for i from 3 to n do

#вычисление значений методом Адамса y[i+1]:=y[i]+h/24*(55*f(x[i],y[i])-59*f(x[i-1],y[i-1])+37*f(x[i-2],y[i-2])-9*f(x[i- 3],y[i-3]));

end do;

for i from 0 to n do print(x=x[i], y=y[i]);

end do; end proc:

>f:=(x,y)->y+(1+x)*yˆ 2:

>dy/dx=f(x,y);

Проверим работу процедуры > Adams(f,-1,1,1.5,5);

x = 1., y = -1

x = 1.100000000, y = -0.9000000000 x = 1.200000000, y = -0.8199000000 x = 1.300000000, y = -0.7539980778 x = 1.400000000, y = -0.7031569156 x = 1.500000000, y = -0.6568883609

11