int sqrt(a^2+b^2+2*a*b*cos(x)) dx