2D Heat Conduction Equation

Objective: % transient implicit sor method l=1; b=1; nx= 10; ny=10; x= linspace(0,1,nx); y= linspace(0,1,ny); dx= x(2)-x(1); dy= y(2)-y(1); q= 1.4; t= ones(nx,ny); t(:,1)=400; t(1,:)= 600; t(:,end)= 800; t(end,:)= 900; told =t; % the old temperature values in the current time step tol= 1e-4; dt=…