X=input('X(Width)=');
Y=input('Y(lenght)=');
m=100*Y;
n=100*X;
dx=X/n;
dy=Y/m;
R=dx^(2)/dy^(2);
Tsheet = zeros(n,m);
TS = zeros(n,m);
matrix1 = zeros(n,1);
coefficient= zeros(n,n);
solve = zeros(n,1);
Tsheet(:,: )=220;
Tsheet(1,: )=100;
Tsheet(n,: )=200;
Tsheet(:,1)=300;
Tsheet(:,m)=400...