Abstract
In this paper we develop a fourth order finite difference discretization strategy and the full multigrid method to solve Helmholtz equation in two dimensional. We use Newton interpolation method from the coarse grid to the fine grid. Numerical results are given to illustrate the efficiently of this method.