Solving Poisson's equation with Dirac delta function

Hello everybody!

I clearly understand how to solve a standard Poisson’s equation in FF++. But can anyone suggest me how to solve it having the Dirac delta function as a RHS? I’m pretty newbie with FF+ so I’d appreciate any help you can give.

Here is an example: FreeFem-sources/Laplace-RHS-Dirac.edp at master · FreeFem/FreeFem-sources · GitHub.

Thank you for suggestion! It really helped me