/*
    RUNGE-KUTTA 2


    t: t inicial
    y: y inicial
    f: funcion
    N: iteraciones
    h: incremento
*/
function T = rk2(t, y, f, N, h)
    T = [t y]
    for i=1 : N
        k1 = f(t, y) * h
        t = t + h
        k2 = f(t, y+k1)
        y = y + (k1 + k2) / 2
        T = [T; t y]
    end
endfunction

/*
    EJEMPLO
    y' = 3ty(t)
    f = 3ty
    y(2) = 6
    N = 100
    h = 1/100

    Aproximar y(2.75)
*/

function y = f1(t, y)
    y = 3 * t * y
endfunction

/*
    T = rk2(2, 6, f1, 100, 1/100)
    busco en T la entrada que tenga 2.75 y esa sera la aproximacion
    de y(2.75)
*/