2d heat conduction

Objective: it has been asked to solve the 2d heat conduction proplem the main equation governing this heat transfer phenomen is: `(dT)/dt -alpha((d^2T)/dx^2 + (d^2T)/dy^2) = 0` first the equatiom need to be descretsized in order to be able to solve it which can be sumarize in the following steps: `(T_i^(n+1)-T_i^(n))/dt =alpha*((T_(i-1,j)-2*T_(i,j)+T_(i+i,j))/(dx^2))^n+alpha*((T_(i,j-1)-2*T_(i,j)+T_(i,j+1))/(dy^2))^n`…