Solving the steady and unsteady 2D heat conduction problem

Objective: %code to solve second order steady state two dimensional equation %d^2t/dx^2+d^2t/dy^2=0; nx=10;ny=10; dx=(1/nx); dy=(1/ny); T=ones(nx,ny); for i=1:nx T(i,1)=400; T(i,ny)=800; T(1,i)=600; T(nx,i)=900; end for i=2:9 for j=2:9 T(i,j)=400; end end Told=T; k=(2*(dx^2+dy^2)/(dx^2*dy^2)); eor=100; tol=1e-4; iteration=3; %code…