А почему дивиргенция равна 3Х, скорее 3*(R^2) или вы так и считали? Потом пределы для R, 0=<R=<1/2
Еще, угол фи будет меняться от нуля до двух пи... Попробуйте...

так ведь по условию X^2 + Y^2 + Z^2 = X, а X в сферических координатах равен R*cosU*sinO
0 ≤ fi ≤ 360° — угол между осью X и проекцией отрезка, соединяющего начало координат с точкой, на плоскость. поэтому от -pi/2 до pi/2 (исправьте, если где-то ошибся)