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

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

дом Эйлера

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