rk4_runge kutta 4_cantidad de iteraciones y paso.sci

function T=rk4(t,y,F,n,h)
    T=[t    y];
    for i=1:n
        k1=F(t,y)*h;
        t=t+h/2;
        k2=F(t,y+k1/2)*h;
        k3=F(t,y+k2/2)*h;
        t=t+h/2;
        k4=F(t,y+k3)*h;
        y=y+(k1+2*k2+2*k3+k4)/6;
        T=[T;t  y];
    end
endfunction

//t es el valor inicial del tiempo
//y es la imagen del tiempo inicial
//F es la funcion a aproximar en cierto valor
//n es la cantidad de iteraciones a realizar
//h es el paso entre dos tiempos consecutivos