Not yet reviewed by Pith; the record is open.

In this paper, using the approximate particular solutions of Helmholtz equations, we solve the boundary value problems of Helmholtz equations by combining the methods of fundamental solutions (MFS) with the methods of particular solutions (MPS). Then the initial boundary value problems of the time dependent diffusion and wave equations are discretized numerically into a sequence of Helmholtz equations with the appropriate boundary value conditions, which is done by either using the Laplace transform or by using time difference methods. Then Helmholtz problems are solved consequently in an iterative manner, which leads to the solutions of diffusion or wave equations. Several numerical examples are presented to show the efficiency of the proposed methods.