(mu / 2) * integrate integrate integrate (((e * h) / (m * rho^3))^2 * (cos^2(theta) + sin^2(theta)) * r^2 * sin(theta)) dr dtheta dphi over 0 to rho over 0 to 2 * pi over 0 to pi