solving the 2-D heat conduction equation

Objective: matlab code as plots are atteched below: clear all close all clc % steady state heat conduction n = 11; T = zeros(n,n); % initial temperature condition T(:,1)=400; % boundry condition left T(:,n)=800; % boundry condition right T(1,:)=600; % boundry condition top T(n,:)=900; % boundry condition bottom T_old1=T;T_old2=T;T_old3=T;…