[求助]一个问题

在看葛德彪、闫玉波的《电磁波是与有限差分方法》的有关UPML 的递推时遇到 zg)-RCG  
一个问题： G+Gd;`4  
在对电场分量ex进行递推时 公式为 yc ize2>q  
ex(n+1)=c1*ex(n)+c2*dx(n+1)-c3*dx(n)------------------------------括号里的n和n+1表示时刻，c1、c2、c3表示系数 Buazm3q8H  
!I3_KuJ5  
在对ex（n+1）的递推时用到了dx（n+1），既然两个的时刻相同，那该怎么去进行递推呢。如果先算出了dx（n+1），则dx（n）会被dx（n+1）替代了啊 MBlhlMyI  
;Ohabbj*  
纠 .. E=P
mOw7b  
'#JC 6#X   

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

用一个中间变量给它存储起来不就行了吗？

问题是 RZKdh}B?\  
直接这样做就对了吗？

ex(n+1)=c1*ex(n)+c2*dx(n+1)-c3*dx(n)------------------------------括号里的n和n+1表示时刻，c1、c2、c3表示系数 mNc?`G_R  
F!
N D  
从这个公式上来看，ex在推导n+1步的时候，需要n时刻，以及n+1时刻的dx值，所以你需要两个数组来保存dx，比如dx[ ] ,以及pre_dx[ ]。 NU]+ {7  
+-s$Htx  
步骤差不多这样吧 [9hslk  
。。。 E:_m6 m  
pre_dx = dx ; //save pre_dx lKtA.{(  
dx[ ] =.... //update dx to n+1 0.MD_s0)>  
ex[ ] = ...dx[ ] + ...pre_dx[ ] BVJ6U[h`  
。。。 ( o(,;  
<./r%3$;7  
   手边没书，只从公式上来看，应该是这样的吧。

谢谢，我这样做了。可是还是有发散的问题。我自己都不知道该怎么办好 了 .=kXO{>  
不知道版主，你能不能帮我看看我的程序呢？ r bfIH":  
我在这里附上我的程序和葛老的书

这是剩下的部分