sum1.m

用fdtd方法 模拟集成光学离子交换中扩散离子浓度的分布 进而求得折射率的分布。更创新性加入了 可变扩散系数 条件（绝对原创）

%一次交换主程序
%=======================================================
%先输入计算参数
halfwindow=5;t=3000;high=30;halfwidth=20;
%调用程序进行计算
diffus=moon_diffu(halfwindow,t,high,halfwidth);
fdiffus=fliplr(diffus);
sumdiffus=[fdiffus,diffus];


x=-19:20;y=(0:29)';[X,Y]=meshgrid(x,y);subplot(1,2,1),mesh(X,Y,sumdiffus);subplot(1,2,2),contour(X,Y,sumdiffus,5);axis square;axis([-20,20,0,30])