for k=1:K+1
t(k)=(k-1)*tau; %时间网格,tau为时间步长
end

for ii=1:101
for jj=1:101
for m=1:21
x(ii)=-X/2+(ii-1)*h; %x,y,z为空间坐标，h为网格尺寸。X，Y为常数
y(jj)=-Y/2+(jj-1)*h;
z(m)=(m-1)*h1;
T(ii,jj,m,1)=0; %初始温度设为环境温度T0=0
g(ii,jj,m)=200;
end
end
end

for k=1:20
for jj=1:101
for m=1:21
D(1)=(3*T(1,jj,m,k)+g(1,jj,m))./(b*(1-h*He)+d); % b、d为已知常数
for ii=2:100
D(ii)=(3*T(ii,jj,m,k)+g(ii,jj,m)-b.*D(ii-1))./(b.*C(ii-1)+d);% C（ii）为已定义函数
end
end
end
end
运行时出现错误信息：Index exceeds matrix dimensions! 在''D(1)=(3*T(1,jj,m,k)+g(1,jj,m))./(b*(1-h*He)+d); ''
我不知道这怎么回事，请哪位高手帮忙看一下，谢谢！