After the successful applications of the combination of the method of fundamental solutions (MFS), the method of particular solutions (MPS), and the eigenfunctions expansion method (EEM) to solve 2D homogeneous and nonhomogeneous diffusion equations by Young et al. (Young et al., Numer Meth Part Differ Equat 22 (2006), 1173), this article intends to extend the same fundamental concepts to calcu...