integrate (1+a*cos u)^(-3/2) from 0 to pi