2 - D heat conduction

Objective: Find the code below clear all close all clc %domain %square %dx = 3e-3; %dy = 3e-3; %no of nodes nx = 52; ny = 52; x = linspace(0,1,nx); y = linspace(0,1,ny); dx = 1/nx; dy = 1/ny; T = zeros(nx); %steady state solving to get initial condtn %TP = (((TL+TR)/dx^2) + ((TT+TB)/dy^2)) / K_steady; K_steady = 2*((dx^(2)+dy^(2))/(dx^(2)*dy^(2)));…