    init;
    for N = [1000] % N = [50, 100, 500, 1000, 5000, 10000]
        ticBorn = tic;
        seti.nOut = 1;
        seti.pdaN = N;
        seti.k = 70;
        disp(' ');
        fprintf('pdaN = %i.\n',seti.pdaN);
        seti = setData(seti);
        seti = setRecon(seti);
        seti = recon(seti);

        tocBorn = toc(ticBorn);
        fprintf('Elapsed time is %05.1f min.\n',tocBorn/60);

        figure(1); imagesc(real(seti.G(seti.qROIcomp))); colorbar; axis xy;
        set(gca,'FontSize',20); axis square; caxis([-0.15 +1.42]);
        filename = sprintf('guideBornInv_%i.eps',seti.pdaN);
        print(1,'-depsc',filename);

        fprintf('rel. dis. %.4f...\n', seti.dis(seti.iOutStop));
        fprintf('rel. err. %.4f.\n', seti.err(seti.iOutStop));
    end