close all;
clear;
clc;


%Image dimension
n = 128;

%Get filter
sz = 3;
h = fspecial('disk',sz); %-> I suggest to use this filter for the convolution operator

%Parameter
alpha = 0.01;
maxit = 1000;
epsilon = 0.01;

delta = 0.02; %Noise level

%Set stepsize
L = 1 + 8*alpha/epsilon;
tau = 2/(L+1);

%Generate input image
orig = phantom(n);

%Generate data
y = 1 %INPUT HERE

%Iterate
u = zeros(n,n);
energy = zeros(maxit,1);

for k=1:maxit

    u = 1 %INPUT HERE

    energy(k) = 1 % INPUT HERE

end

figure (1)
out = [orig,y,u];
imshow(out)
title(['Original | Data | Reconstruction'])

figure (2)
plot(log10(energy))
title(['Objective value'])




