Improved Euler's method
显示 更早的评论
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
2 个评论
James Tursa
2021-3-2
Are you getting errors or an answer you didn't expect? What behavior are you seeing that prompted you to ask a question?
Aaron Millan
2021-3-2
回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Linear Algebra 的更多信息
2021-3-2
2021-3-2
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
