求帮找错

function T
for i=1:16
T(i,1)=30;
end
for j=1:12
T(1,j)=30;
end
for p=6:16
for q=6:12
T(p,q)=0;
end
end
T(16,2:5)=[24,18,12,6];
T(2:5,12)=[24,18,12,6];
while 1
for i=2:15
for j=2:5
eps=1;
temp=T(i,j);
T(i,j)=(T(i-1,j)+T(i+1,j)+T(i,j-1)+T(i,j+1))./4;
eps=T(i,j)-temp;
end
end
if eps<1e-8
break;
end
end
while 1
for i=2:5
for j=2:11
eps=1;
temp=T(16,j);
T(16,j)=(2.*T(15,j)+T(16,j-1)+T(16,j+1))./4;
eps=T(16,j)-temp;
end
end
if eps<1e-8
break;
end
end
while 1
for i=2:5
temp=T(i,12);
T(i,12)=(2.*T(i,11)+T(i-1,12)+T(i+1,12))./4;
eps=T(i,12)-temp;
end
if eps<1e-8
break;
end
end
pcolor(T');shading interp
figure,contour(T',10);
end

解释一下 这是个二维墙角导热的数值模拟程序，学能源或者热工的同学应该有很多做过的吧。 然而意义其实不是很重要了， 重要的是为何这个程序运行之后没有任何数据输出，我想要的是这个12*16的矩阵每个点上都应该有一个数 代表该点的温度。求大神帮忙解答，错出在哪里

自顶 勿沉

30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000
30.0000 28.0681 25.7243 22.2202 15.6323 0 0 0 0 0 0 12.0000
30.0000 26.5480 22.6091 17.5240 10.3091 0 0 0 0 0 0 6.0000
30.0000 25.5148 20.6401 14.9576 8.0802 0 0 0 0 0 0 3.0000
30.0000 24.8710 19.4788 13.5861 7.0541 0 0 0 0 0 0 0.7500
30.0000 24.4904 18.8181 12.8539 6.5500 0 0 0 0 0 0 0
30.0000 24.2725 18.4495 12.4613 6.2919 0 0 0 0 0 0 0
30.0000 24.1502 18.2459 12.2498 6.1566 0 0 0 0 0 0 0
30.0000 24.0824 18.1342 12.1355 6.0844 0 0 0 0 0 0 0
30.0000 24.0450 18.0732 12.0736 6.0457 0 0 0 0 0 0 0
30.0000 24.0245 18.0398 12.0399 6.0247 0 0 0 0 0 0 0
30.0000 24.0133 18.0215 12.0216 6.0134 0 0 0 0 0 0 0
30.0000 24.0071 18.0115 12.0115 6.0071 0 0 0 0 0 0 0
30.0000 24.0036 18.0059 12.0059 6.0036 0 0 0 0 0 0 0
30.0000 24.0015 18.0025 12.0025 6.0015 0 0 0 0 0 0 0
30.0000 24.0096 18.0352 12.1264 6.4655 1.7324 0.4642 0.1244 0.0333 0.0089 0.0022 0

倒数第二行加一个T

figure,contour(T',10);
T
end