MATLAB有限差分法程序

MATLAB有限差分法求解拉普拉斯（Laplace）方程，长直接地金属矩形槽内部电位分布

d=zeros（110）;
h=0;
for e=1:0.1:1.9
a=zeros（2141）;
a（1:）=100;
b=zeros（1939）;
c=eye（1939）;
count=1;
g=0;
while（count==1）
    m=0;
    for i=2:1:20
        for j=2:1:40
            b（i-1j-1）=a（ij）;
            a（ij）=a（ij）+e*0.25*（a（i-1j）+a（i+1j）+a（ij-1）+a（ij+1）-4*a（ij））;
            c（i-1j-1）=abs（a（ij）-b（i-1j-1））;
            if（c（i-1j-1）<0.001）
                for k=1:1:19
                    for n=1:1:39
                        if（c（kn）<0.001）
                            m=m+1;
                        else
                            m=0;
                            break;
                        end
                    end
                    if（m==0）
                        break;
                    end
                    if（m==741）
                        count=0;
                    end
                end
            end
            if（count==0）
                break;
            end
        end
        if（count==0）
            break;
        end
    end
   g=g+1;
end
h=h+1;
e
d（1h）=g;
end
d


                        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1110  2018-04-21 14:40  有限差分法求解PISSION方程（高斯与超松弛）\chaosongchi.m
     文件        1028  2018-04-21 16:05  有限差分法求解PISSION方程（高斯与超松弛）\gaosi.m
     目录           0  2018-04-22 19:13  有限差分法求解PISSION方程（高斯与超松弛）\