Improved Euler's method

15 次查看（过去 30 天）

Aaron Millan 2021-3-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/760006-improved-euler-s-method

评论： Aaron Millan 2021-3-2

Hello! I am working on a code I tried to run last year for Euler's improved method, which uses the height of a trapezoid to find a definite integral value. My code is written below.

the first line of my code is exactly what we went over in class, so I know that part is correct; do you have any suggestions or see anything wrong with my code? Many thanks!

function y = ieuler(f,y0,t0,tn,N)
h = (tn-t0)/N; 
t = t0:h:tn;
y = zeros(length(y0), N+1);
y(:,1) = y0;
for k = 1:N
    e = y(:,k) + h*f(t(k), y(:,k)) ;
    y(:,k+1) = y(:,k) + .5*h*(f(t(k), y(:,k))+ f(t(k+1),e));
    
end