Steady and unsteady 2D heat conduction problem

Objective: A. Explicit Solver - Unsteady (Gauss Jacobi) % Explicit Solver - Unsteady clear all close all clc % Initial Variables nx = 40; ny = nx; L = 1; % unit square domain x = linspace(0,L,nx); y = linspace(0,L,ny); dx = 1/(nx-1); dy = 1/(ny-1); dt = 0.01; alpha = 1e-4; % CFL CFL = alpha*dt/dx^2; % Temperature Matrix and BC T…