This chapter focuses on the project of finding the potential for a given distribution of charges in a two-dimensional system, which does not possess any symmetrical properties, an extension of the cylindrical potential problem discussed in the previous chapter. Using a method of minimising a functional, specifically the Gauss–Seidel method of iterative minimisation, the Poisson’s equation is adjusted to a 2D case, neglecting one partial derivative in Cartesian coordinates. We subsequently derive a discretised form of the functional, leading to a multi-variable function, following which the problem can be solved using the Gauss–Seidel iterative method. The numerical method discussed here is the finite elements method (FEM), with an emphasis on the need for a specific sequence for updating values to optimise computation efficiency. The discussion sheds light on the importance of the uniqueness of solutions in electrostatic systems, thereby exploring a fundamental question in electrostatics. The concluding part of the chapter provides an outline of a numerical algorithm for the problem, suggesting potential modifications and points for further exploration.