with ε ε being a diagonal D × D tensor, this differential equation is reformulated as an integral equation. Integration the above equation over a D-dimensional region Ω Ω , and using Gauss theorem, then the integral formulation is obtained