Need help with FDFD (Frequency domain) Method

Hi Dear all
I have wrote a FDFD code in MATLAB. But, it doesn't work fine and there is reflections from the right and upper PML. The program is very simple and can be understood with a glance. Please, if you can, have a look to my code to see what is the problem. I have checked my code many times but I don't know what is the problem. Tell me if you can guess the problem.
Regards