The first step in the work was to solve the theoretical basis for the fixed method. This method previously required the formulation of an expression for the potential distribution in a multilayer medium. The medium was defined as multilayer where the first layer was air. The basic expression that had to be determined was the potential distribution, so the electric potential could be given by the expression of the Helmholtz differential equation. Immediately in the first calculations, an assumption for the wave constant was introduced in order to bypass the Sommerfeld integrals, and the Poisson differential equation was obtained. The retardation factor was included in the expression later. Then, a solution to the Poisson differential equation was derived and an expression was obtained that contained two kernel functions. Further calculations were actually based on the adaptation and approximation of these two kernel functions at the boundary values, and 25 points were taken for the approximation, and finally the retardation factor was introduced. The second step was the complex method and one of its forms, which is called the matrix bundle least squares method. This method focused on the Green’s function, whose expression was written in a general form. The key point here is that it was necessary to take into account that this is a cylindrical coordinate system and the Green’s function was written for such a case. The expression for the closed Green’s function was obtained using the Lipschitz integral and introducing the above substitution. In the end, an expression was obtained in which almost every part of this expression represents a quasi-static image and the electromagnetic field was considered quasi-static, which is why this method was called the quasi-static method of complex images. The matrix bundle least squares method is a method that gives approximately accurate results despite the occurrence of numerical errors and noise. By using matrices and substitution, the error and noise are nullified and a more precise method is obtained than the ordinary matrix bundle method. Finally, a numerical solution is given, as well as graphical examples and the procedure for calculating the expression for a three-layer and four-layer medium, and a conclusion is given on the relationship between these two methods.