\int (R^2*z-z^3/3)*(z*sqrt(-b^2*(b^2-R^2+z^2)/((R^2-z^2)^2))/(b^2-R^2+z^2)) dz