In this article, we implement a new analytical techniques, He's variational iteration method for solving the linear Helmholtz partial differential equation. In this method, general Lagrange multipliers are introduced to construct correction functionals for the problems. The multipliers in the functionals can be identified optimally via the variational theory. The initial approximations can be freely chosen with possible unknown constants, which can be determined by imposing the boundary/initial conditions. The results compare well with those obtained by the Adomian's decomposition method.