% jzheats
clear,clc
nx=10;ny=10;hx=1/nx;hy=1/ny;
x=0:hx:nx*hx;
by0=100;
byn=100;
y=0:hy:ny*hy;
bx0=0;
bxn=0;
[T,iter]=HeatS(nx,hx,ny,hy,bx0,bxn,by0,byn);
T=rot90(T,1);
iter
% Graph
T=flipud(T);
[c,h]=contour(T);clabel(c,h),colorbar
title('Temperature Distribution')...