%模拟 TM 波传播的 MATLAB 源程序
%function main
fprintf('please wait patiently ,some files are being created ...');
v=3*10^8;
delta_t=2.5*10^(-11);
delta_s=2*v*delta_t;
steps_t=100;
cicle_N=2*steps_t+1;
sourse_location_x=51;
sourse_location_y=51;
grids_M=100;
grids_N=100;
Ez_M=grids_M+1;
Ez_N=grids_N+1;
Hx_N=grids_N+1;
Hy_M=grids_M+1;
Hy_N=grids_N;
%初始化
temp=zeros(Ez_M,Ez_N);
temp(sourse_location_x,sourse_location_y)=sourse(0);
Ez_data_getfile;
fid=fopen(Ez_data_file{1},'w');
fprintf(fid,'%e',temp);
fclose(fid);
%value_Ez(1)=temp; %cicle=1
Hx_data_getfile;
fid=fopen(Hx_data_file{1},'w');
fprintf(fid,'%e',zeros(Hx_M,Hx_N));
fclose(fid);
%value_Hx(1)=zeros(Hx_M,Hx_N); %cicle=2
Hy_data_getfile;
fid=fopen(Hy_data_file{1},'w');
fprintf(fid,'%e',zeros(Hy_M,Hy_N));
fclose(fid);
%value_Hy(1)=zeros(Hy_M,Hy_N); %cicle=2
for cicle_n=3:cicle_N
%*****************计算 Ez********************
%下面将所有的电场分量都更新
if mod(cicle_n,2)==1
