In this study, by maximizing the effective condition number in the resultant algebraic system via the method of fundamental solutions for 3D Laplace operator, the optimal fictitious boundary with offset between the source and collocation points can be determined. The simplicity, efficiency and high accuracy of this method is illustrated in several numerical examples
