Numerical solution of the Helmholtz equation for two-dimensional polygonal regions