Método de RK3:      

q1 = h f( un  ,   tn )

q2 = h f( un + q1 / 2         ,   tn+h/2 )

q3 = h f( un - q1 + 2*q2   ,   tn+h/2)

u(n+1) =  un  +  1/6 * (q1+4*q2+q3 )