taylor.m

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

function [tout,yout] = taylor(f,ft,fy,tspan,y0,N)

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