Implementing Heun's first order ODE

This what I have so far
function[t,y]=heun(f,a,b,h,y0)
n=(b-a)/h;
t=[a];y=[y0];
for k=2:n+1,
t=[t; a+(k-1)*h];
y=[y; y(k-1)+h*f(t(k-1),y(k-1))];
end
But I need it to do the corrector cycle equations to make the y values more accurate. This is the formula for corrector cycle
%

回答(0 个)

类别

帮助中心File Exchange 中查找有关 Programming 的更多信息

提问:

2017-4-13

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by