ME685 Homework4
ME685 Homework4
end
TempDistribution=T;
% Create a new figure
figure;
subplot(2, 3, 1);
% At t = 0.05
contourf(TempDistribution(:,:,Nt/5));
title('t = 0.05');
subplot(2, 3, 2);
% At t = 0.10
contourf(TempDistribution(:,:,(2*Nt)/5));
title('t = 0.10');
subplot(2, 3, 3);
% At t = 0.15
contourf(TempDistribution(:,:,(3*Nt)/5));
title('t = 0.15');
subplot(2, 3, 4);
% At t = 0.20
contourf(TempDistribution(:,:,(4*Nt)/5));
title('t = 20');
subplot(2, 3, 5);
% At t = 0.25
contourf(TempDistribution(:,:,end));
title('t = 0.25');
for i=2:Ny-1
b(2:end-1)=(2*(1-lam)*Temp(i,2:end-1)+lam*Temp(i-
1,2:end-1)+lam*Temp(i+1,2:end-1)+q)';
b(1)=Temp(i,1);b(end)=2*(1-
lam)*Temp(i,end)+lam*Temp(i-1,end)+lam*Temp(i+1,end)-
2*lam*a*hx+q;
Tempf(i,:)=(Ax\b)';
end
Tempf(Ny,:)=T(end,:,k-1);
T(:,:,k)=Tempf;
end
TempDistribution=T;
% Create a new figure
figure;
subplot(2, 3, 1);
% At t = 0.05
contourf(TempDistribution(:,:,Nt/5));
title('t = 0.05');
subplot(2, 3, 2);
% At t = 0.10
contourf(TempDistribution(:,:,(2*Nt)/5));
title('t = 0.10');
subplot(2, 3, 3);
% At t = 0.15
contourf(TempDistribution(:,:,(3*Nt)/5));
title('t = 0.15');
subplot(2, 3, 4);
% At t = 0.20
contourf(TempDistribution(:,:,(4*Nt)/5));
title('t = 20');
subplot(2, 3, 5);
% At t = 0.25
contourf(TempDistribution(:,:,end));
title('t = 0.25');