2D thermal conduction equation solver

Objective: Transient Implicit clc clear all close all sqn=10; % boundary conditions u=ones(sqn); for a=1:sqn u(a,1)=500; u(a,sqn)=800; u(1,a)=600; u(sqn,a)=900; end uold=u; uprevdt=u; uprime=u; % assigning constant k info and time loop specs dt=0.1; alpha=0.1; %inversely proportional to no. of iterations %end time tend=2; dttend=tend/dt;…