voko said:
The only differential thing in your approach is the thickness of the spherical surface ##dr##. Obviously ##r## and ##\theta## are functionally interdependent, so you could find (in theory) ##r = f(\theta)##, then you would be able to do ##dV = A(r) dr = A(f(\theta)) f'(\theta) d\theta##, but I see no reason you would want to do that, as that only seems to complicate matters for no obvious benefit.
Yes, I am not interested in expressing r in terms of theta as it would, as you say, definitely complicate the things. But I can't see why it is Adr instead A(rdθ), can you please explain?
Using dV=Adr, the charge contained in this volume is:
$$dq=\frac{\rho\pi r}{d}\left(R^2-(r-d)^2\right)dr$$
I set up the potential energy integrals in the following manner:
$$U=\int_{d-R}^R \frac{kQ\, dq}{2R}\left(3-\frac{r^2}{R^2}\right)+\int_R^{d+R} \frac{kQ\, dq}{r}$$
I solve the integrals separately.
First integral is:
$$\int_{d-R}^R \frac{kQ}{2R}\left(3-\frac{r^2}{R^2}\right)\frac{\rho\pi r}{d}\left(R^2-(r-d)^2\right)dr$$
I solved the integral using wolfram alpha and got:
$$-\frac{kQ\rho \pi}{2Rd}\frac{d(d-2R)^2(d^3+4d^2R-18dR^2-48R^3)}{60R^2}$$
Second integral is:
$$\int_R^{d+R} \frac{kQ\rho \pi}{d}\left(R^2-(r-d)^2\right)dr$$
Solving:
$$-\frac{kQ\rho \pi}{d}\frac{1}{3}d^2(d-3R)$$
These look dirty, am I going in the right direction?