RK4 routine for rotational model still needs correcting,
currently the solution doesn't preserve normalisation