The Method of Fundamental Solutions (MFS) is a boundary-typemeshless method for the solution of certain elliptic boundaryvalue problems. In this work, we investigate the properties of thematrices that arise when the MFS is applied to theDirichlet problem for Laplace's equation in a disk. In particular,we study the behaviour of the eigenvalues of these matrices andthe cases in which they vanish. Based on this, we propose amodified efficient numerical algorithm for the solution of theproblem which is applicable even in the cases when the MFS matrixmight be singular. We prove the convergence of the method foranalytic boundary data and perform a stability analysis of the methodwith respect to the distance of the singularities from the originand the number of degrees of freedom. Finally, wetest the algorithm numerically.
