steady state analysis

Objective: jacobi method:- function[out] = jm;% definig a function with time lx = 1;ly = 1;nx = 10;ny = 10; x = linspace(0,lx,nx);y = linspace(0,ly,ny);dx = x(2)-x(1);dy = y(2)-y(1);[xx yy] = meshgrid(x,y); % initilizing the temparature field t = zeros(nx,ny); % boundary conditions% temparature the top surfacet(1,:)=900;%…