%u0 = imread('saturn_noise.bmp');
u0=imread('fic.noise.bmp');



% coefficient of the TV norm (needs to be adapted for each image) 
lambda=0.028; 
%  number of iterations (depends on the image) 
IterMax=100;
% call ROF TV_L2 model to denoise
u = TV_L2(u0,lambda, IterMax);