integrate(sqrt(x*d -x^2))