Numerical simulation of 2D heat conduction equation for steady and unsteady state systems

Objective: The 2D conduction equaiton is given by :- `(del^2T)/(delx^2)+(del^2T)/(dely^2)=alpha(delT)/(delt)` where alpha represents the thermal conductivity Descritising the gicen equation with central difference scheme for space derivatives and forward difference for time derivatives, We get `(delT)/(delt)=(T_(i,j)^(n+1)-T_(i,j)^n)/(/_\\t)`…