integrate y^2 * (pi - arccos(1 - y/a) + sqrt(2a*y - y^2)) wrt y