d/dt [t*sqrt(t^2-a^2)/sqrt(1+t^2) + EllipticE[ArcCot[t], 1+a^2] + a^2 EllipticF[ArcCot[t], 1+a^2]]