d/dt(sin t cos t / sqrt(a^2sin^2 t + b^2 cos^2 t))