Stablity analysis

Objective: Here plots at different dt is shown for both implicit and explicit methods reducing the dt value from 1 to 0.001. At around 0.45 CFL number the solution starts becoming unstable. close all clear all clc nx=10;ny=nx;dt=0.01;nt=2000;s=1e-4; x=linspace(0,1,nx);y=linspace(0,1,ny); z=fliplr(y);dx=(x(2)-x(1));dy=dx;f=(s*dt)/(dx^2);…