Heun Method in Matlab


If  we want to implement Heun method in matlab, we can write down like the code below:

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

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s