A solver for Helmholtz equation using Jacobi in 2d