it requires to write the Matlab program for the question.

clc
clear all
close all

I0=double(imread('heart1.jpg')); %%% Reading image
I1=I0;
I2=I0;
[m,n,~]=size(I0);
n1 = 10; %%% No. of iterations
n2 = 100;
n3 = 1000;
dt=0.1;
gamma = 1;
subplot(2,2,1), imshow(uint8(I0))
title('Original image');

for t = 0:dt:n1
u_x = I0(:,[2:n n],:) - 2*I0 + I0(:,[1 1:n-1],:); % FD approximation for u_x
...