A multigrid Method (MGM) is employed to calculate five fixed solutions to the first boundary-value problems of Helmholtz’s Euler equation,which are compared to the super relaxation iteration with variable coefficients,Evidence suggests that the required time of the MGM is less than that of the iteration for the same accuracy with the time/efficiency ratio greatly enhanced with increased grid number,thereby exhibiting advantages in finding numerical solutions of a large-scale difference equation.