clear;
N=10;

% Manifold diminsions
MB=8*N/10;
MT=9*N/10;
ML=4*N/10;
MR=9*N/10;

% C.C. diminsions
CB=2*N/10;
CT=4*N/10;
CL=2*N/10;
CR=3*N/10;

% Initialize
maxChange=1000;
T(1:N,1:N)=0;
newT(1:N,1:N)=0;

% Set boundry values for manifold,  C.C., and walls
T(MB:MT,ML:MR)=400;
T(CB:CT,CL:CR)=1200;
T(1,1:N)=80;
T(N,1:N)=80;
T(1:N,1)=80;
T(1:N,N)=80;


while maxChange>1
        % Your Code goes here

	maxChange=sum(sum(abs(T-newT)));
	T=newT;

end
% Plot the result here using "mesh".

