For a spherically symmetric function , the Laplacian in spherical coordinates is
Integrating from zero to and using regularity at the origin gives
For the stated density, integration and the boundary condition give
The two pieces and their first derivatives agree at . If another solution existed, its difference from would be harmonic with zero boundary data, so part (a) proves Uniqueness of the Dirichlet problem.
On the shell , normalize the outer radial part to obtain the harmonic function
By the Dirichlet principle, this function minimizes the energy among all with the same boundary values. Since
its energy is
Therefore