side4_aa.m

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

function sd4p=side4_point(mat,i,j,dt,d,db,alf,h)
sd4p=dt.*(...
    d(i,j)./(h^2.*(1-alf(i,j).*mat(i,j))).*(mat(i-1,j)+mat(i+1,j)+mat(i,j+1)-3*mat(i,j)+alf(i,j)./(1-alf(i,j).*mat(i,j)).*((mat(i+1,j)-mat(i,j)).^2+(mat(i,j+1)-mat(i,j)).^2))...
    +...
    (d(i+1,j)-d(i,j))./h.*(1-mat(i,j))./(1-alf(i,j).*mat(i,j)).*((mat(i+1,j)-mat(i,j))./h.*(1+alf(i,j)./(1-alf(i,j).*mat(i,j))))...
    +...
    (db(i+1,j)-db(i,j))./h.*mat(i,j).*(1-alf(i,j))./(1-alf(i,j).*mat(i,j)).*((mat(i+1,j)-mat(i,j))./h.*(1-alf(i,j).*(1-mat(i,j))./(1-alf(i,j).*mat(i,j))))...
    );