“…Different ways can be used to obtain solutions to (1), among which the FEM is a representative one with flexibilities to deal with irregular boundaries [14]. In the framework of FEM, the domain Ω is divided into P elements and joined at N vertex nodes.…”