function v = rk2(V,t,k,F)

    k1 = feval(F,V,t);                                 
    k2 = feval(F,V + k*k1,t+k);
                            
    v = V + 0.5*k*(k1 + k2);