    init;                   % Initialization (addpath...)
    seti = struct;
    seti = setData(seti);   % Set data: experimental set-up, contrast, simulate data
    seti = setRecon(seti);  % Settings for the variational reconstruction
    seti = recon(seti);     % Variational reconstruction (process)

    figure(1); imagesc(real(seti.G(seti.qROIexact))); axis xy; colorbar;
    %set(gca,'FontSize',20); axis square; caxis([-0.15 +1.42]); print(1,'-depsc','guideRecon1.eps');

    figure(2); imagesc(real(seti.G(seti.qROIcomp))); axis xy; colorbar;
    %set(gca,'FontSize',20); axis square; caxis([-0.15 +1.42]); print(2,'-depsc','guideRecon2.eps');

    fprintf('Reconstruction stopped after %i outer iterations...\n', seti.iOutStop);
    fprintf('... with relative discrepancy %.4f...\n', seti.dis(seti.iOutStop));
    fprintf('... and relative error %.4f.\n', seti.err(seti.iOutStop));