“…Same as the classic implicit continuum models, 10,15,19,20,22,24,36 in order to handle the singularities caused by the Dirac delta distributions, it is a common sense now to use the solution decomposition schemes to isolate the singularities. For the presentation purpose in this paper, we follow 10,15 to decompose the electrostatic potential as follows where is given by is a solution of the following Laplace equation and together with satisfies the following coupled equations …”