integrate_0^s (a^2*cos(t)*sin(t))/(pi*((a^2−1)*cos(t)^2+1)^2) dt