I have the necessity to compute approximate solutions to the equation (pedices denote derivatives) $$\tag{W} u_{tt}-u_{rr}-\frac{n-1}{r}u_r=f(t, r)\quad t\in\mathbb R,\ r\in [0, R]$$ subject to the initial conditions $$ u(0, r)=0,\quad u_t(0, r)=0$$ and the boundary conditions $$ u_r(t, 0)=u(t, R)=0.$$ The equation (W) is obtained by writing the wave equation $u_{tt}-\Delta u = f$ in spherical polar coordinates in $\mathbb R^n$, under the assumption that the source term $f$ is radially symmetric at all times. I am especially interested in the value $n=5$ for the spatial dimension.
Can you give me to some numerical scheme to compute such a solution numerically?
Unfortunately, I am not familiar with numerical analysis. I would be glad to receive an answer or literature pointer containing some details on how to implement such a method.
Thank you for reading.