where ε ε x and ε ε y are the two components of the diagonal of the tensor ϵ ϵ . Discretizing the right-hand side of the Poisson equation is fairly simple. ρ ρ is discretized on the same grid nodes, as it's been done for φ φ .