%Create a script that:
%  ● reads in this image: AndreKertesz_Paris_ManOnBicycle_part.jpg,
%  ● converts it to grayscale,
%  ● calls your anisotropic diffusion function on it (parameters eg.
%  iternum: 200, K: 4, ᶝ: 0.1),
%  ● displays the input image and the diffused image.

im = imread("AndreKertesz_Paris_ManOnBicycle_part.jpg");
if (size(im,3) ~= 1) 
    im = rgb2gray(im);
end

figure
subplot(1,3,1)
imshow(im)
title("Original image")

subplot(1,3,2)
imshow(anisotropic_exp(im,200,4,0.1))
title("Anisotropic diffusion filter (exponential g function)")

subplot(1,3,3)
imshow(anisotropic_quotient(im,200,4,0.1))
title("Anisotropic diffusion filter (qoutient-type of g function)")