Determining the Optimal Source Points in MFS by Minimizing an Energy Gap Functional for 3D Laplace Operator. (21st December 2020)