euler.m

matlab codes for solving ODEs, including euler method, RK methods etc, and some examples.

function [tout,yout]=euler(f,tspan,y0,N)

h=(tspan(2)-tspan(1))/N;
t = tspan(1);
tout = t;
y=y0(:);
yout = y.';
for n=1:N
    y = y+h*feval(f,t,y);
    t = t+h;
    yout=[yout;y.'];
    tout=[tout;t];
end