Acceleration

speed and position:

I have ``` a = dv / dt; v = dp / dt; ```

I am going to use the second order runga kutta method with a = 0, b = 1, alpha = 1 / 2 and beta = 1 / 2.