M=8;
D=zeros(M);
C=zeros(M);
for i=2:M-1
D(1,i)=35;
end
J=ones(M);
J(:,M)=0;
J(M,:)=0;
J(1,:)=0;
J(:,1)=0;
while sum(abs(J(:)))~=0
C=D;
w=1e-5;
for i=M-1:-1:2
for j=2:M-1
D(i,j)=C(i,j)+1.45/4*(C(i-1,j)+D(i+1,j)+C(i,j+1)+D(i,j-1)-4*C(i,j));
if abs(D(i,j)-C(i,j))<w
J(i,j)=0;
end
end
end
end
contour(D);
grid on
- 1
- 2
前往页