二维传热Matlab程序

dx=0.01;          %X节点间距(m)
divx=100;         %X分段数
divy=20;         %Y分段数
a=1;              %传热系数
timestep=0.00001; %时间步长(s)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%定义节点和相对应坐标

for i=1:divx+1
    for j=1:divy+1
        node((i-1)*(divy+1)+j,:)=[(i-1)*dx,(j-1)*dx];
        plot(node((i-1)*(divy+1)+j,1),node((i-1)*(divy+1)+j,2),'o');
        hold on
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%定义边界节点

for i=1:divy+1
    nodeBC(i,1)=i;
end
 
for i=1:divx
    nodeBC=[nodeBC;1+i*(divy+1)];
end

for i=1:divx
    nodeBC=[nodeBC;(i+1)*(divy+1)];
end

for i=2:divy