clear?all;

u=4.*pi.*10.^（-7）;
e=8.854187817.*10.^（-12）;
cig=0;
c=3*10^8;

Imax=200;
Jmax=200;
Kmax=200;
Nmax=800;
d=0.001;
dt=d/2/c;

s=（c*dt-d）/（c*dt+d）;???????????????????????????????%邊界計算公式系數（相鄰）
s1=（c*dt-sqrt（2）*d）/（c*dt+sqrt（2）*d）;??????????????%邊界計算公式系數（斜對角）

Ca=1;
Cb=dt/e/d;
Da=1;
Db=dt/u/d;

Ex=zeros（ImaxJmax）;
Ey=zeros（ImaxJmax）;
Ez=zeros（ImaxJmax）;
Hx=zeros（ImaxJmax）;
Hy=zeros（ImaxJmax）;
Hz=zeros（ImaxJmax）;

A=zeros（ImaxJmax）;?????????????????????????????????%三組記錄上一時刻的數組
B=zeros（ImaxJmax）;
C=zeros（ImaxJmax）;
for?n=1:Nmax
????
????
?????A=Ex;??????????????????????????????????????????%記錄上一時刻Ex????
????for?x=1:60??????????????????????????????????????%計算Ex（1-60）?
????????for?y=2:Jmax
?????????Ex（xy）=Ca*Ex（xy）+Cb*（Hz（xy）-Hz（xy-1））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3481??2017-04-26?14:57??FDTD2withM.m

-----------?---------??----------?-----??----

?????????????????3481????????????????????1