I'm trying to solve the following boundary value problem:
$$\begin{cases}
\Delta u =1 & \{(x,y,z)\in\Bbb R^3:1<x^2+y^2+z^2<9\}\\
u = 0 & \{(x,y,z)\in\Bbb R^3: x^2+y^2+z^2 = 1\}\\
\nabla u\cdot n = 0 & \{(x,y,z)\in\Bbb R^3:x^2+y^2+z^2 = 9\}\\
\end{cases}$$
Here, $n$ is the outward unit normal vector field on $\{(x,y,z)\in\Bbb R^3:x^2+y^2+z^2 = 9\}$.
First change into spherical coordinate $u(x,y,z) = v(\rho,\phi,\theta)$. Since boundary conditions are independent of $\phi,\theta$, $v(\rho,\phi,\theta) = v(\rho)$. $\Delta v=1\iff {\partial^2 v\over \partial \rho^2}+{2\over\rho}{\partial u\over…