I have run your code "Helmholtz.ipynb", the result is here  which seems to be wrong. It is the parameter you set up is incorrent?