gaosi.m

来自「高 斯 迭 代」· M 代码 · 共 27 行

u=zeros(7,7);
m=7;
for n=2:(m-1)
u(7,n)=100;
end
for i=2:(m-1)
    for j=2:(m-1)
        u(i,j)=0;
    end
end
u
N=input('please input N(N<1=Z):N=');
for k=1:N
    v=u;
    for i=2:(m-1)
        for j=2:(m-1)
          u(i,j)=1/4*(u(i-1,j)+u(i,j-1)+v(i,j+1)+v(i+1,j));          
        end
    end
end    
end
contour(u,'b')
xlabel('x');
ylabel('y');
u
t=abs(u-v)