It seems the key is to define an auxiliary independent variable (in attachment the variable "x") and be able to define the Dirac delta function. In an imaginary perspective for a 1D problem we have 2 axis. In example variables "y" and "x" share a same domain.
I would appreciate for any comments, :-)