Runge-Kutta 4/5 ODE solver

function numerical=RKOD(f,a,b,N,o)
h=(b-a)/N;
t=a;
w=o;
t
w
for i=1:N
    K1=h*f(t,w);
    K2=h*f(t+h/2,w+K1/2);
    K3=h*f(t+h/2,w+K2/2);
    K4=h*f(t+h,w+K3);
 
    w=w+(K1+2*K2+2*K3+K4)/6;
    t=a+i*h;
 
    t
    w
end
end