[求助]平面波程序调试遇到问题了～～急，急，急需各位高手师兄指点

我拿出我写的程序的一部分说明一下 x#SE%j?  
根据网格对应关系 ~Sb)i f  

图片:4.jpg

!l#aq\:}~e  
根据下面这个公式 3S_H&>K  

图片:5.jpg

AlDp+"|  
我写的程序如下 ~tfd9,t  
首先定义引入磁场分量的存储空间 }X8P5c!\  
Eincxd_Left(ie_jmin:ie_jmax-1,ke_jmin:ke_jmax)=0； I$sXbM;z=  
然后根据图中的迭代公式写出程序 0U%tjYk(  
for xhz=ie_jmin:ie_jmax-1   FYFlh^}  
       for zhz=ke_jmin:ke_jmax QlB9m2XB  
       hz(xhz,je_jmin-1,zhz)= hz(xhz,je_jmin-1,zhz)-db(1)*Eincxd_Left(xhz,zhz); t=}]4&Yp  
       end _XP}fx7$C  
   end A/a=)su  
中间一维入射波场分量的引入我没贴出来，我认为错误处在这2段程序内 BhAT@%  
请师兄指点一下

那就针对那两个角点进行调试，第一步，第二步。。。，它们的值你先手算出来好了，然后看在哪里和你理论值不一样。

这部分程序，一个学生也在做，做通了才能毕业，就这两天可能做完 &'i.W}Ib!  
到时我让他讨论这个问题

谢谢cem-uestc 老师，我的qq是24360120，目前我自己在这个东西，老师也不懂，周围也没人懂，感觉很吃力，希望能跟您的学生交流一下