Polar - axisymetric poisson formulation with PointSource at origin

I would recommend having a look at this 1 D axisymmetric heat conduction with neumann boundaries and the link within thus post