کد متلب برای یک صفحه با دمای اولیه که هر وجه ان را ناگهان در دمای مختلف قرار می دهند

[maziiiiii]

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;
offset=5;
eps=.01;
while


eps<offset
offset=0;
TS,: )=Tsheet,: );
for
j=2:m-1
coefficient(1,1)=1;
coefficient(n,n)=1;
matrix1(1,1)=Tsheet(1,j);
matrix1(n,1)=Tsheet(n,j);


for i=2:n-1
coefficient(i,i-1)=1;
coefficient(i,i)=-2*(1+Y);
coefficient(i,i+1)=1;
matrix1(i,1)=-R*Tsheet(i,j+1)-R*Tsheet(i,j-1);


end
solve(1,1)=Tsheet(1,j);
solve(n,1)=Tsheet(n,j);


for di=2:n-1
solve(di,1)=(matrix1(di,1)-coefficient(di,di-1)*Tsheet(di-1,j)-coefficient(di,di+1)*Tsheet(di+1,j))/coefficient(di,di);
Tsheet(di,j)=solve(di,1);


end
Tsheet,j)=solve,1);
end
for


ej=1:m


for ei=1:n
offset=offset+abs(Tsheet(ei,ej)-TS(ei,ej));


end
end
end

​