dv/dt=sqrt(2+a^2(1+ sin^2v))