Я использую Maple для реализации метода Эйлера, он требует 5 параметров:
- func_y (общая форма ODE)
- x_value
- y_value
- n (Точные цифры, я не знаю, как управлять точными цифрами)
- h (Размер шага).
Это моя реализация Maple:
my_euler := proc (func_y, x_value, y_value, n, h)
local point_list, d_x, true_f, true_f_v; point_list := [10^(-n), 0, n];
true_f := unapply(func_y, x, y);
d_x := abs(point_list[2]-point_list[3]);
while d_x < point_list[1] do
true_f_v := true_f(x_value, y_value);
y_value := y_value+h*true_f_v;
x_value := x_value+h
end do;
return y_value
end proc
Итак, я не знаю, как управлять точными цифрами, вы можете научить меня, как закончить sh эту часть?