clear all, close all, clc;
fel=inf;
n=1;
T=litenTrapets(@fEXS64, 0,1,n);
n=2;
Tny=litenTrapets(@fEXS64, 0,1,n);
Ttak=Tny+(Tny-T)/3;
iter=1;
while abs(fel)>Ttak*1e-2
        n=n*2;
        T=Tny;
        Tny=litenTrapets(@fEXS64, 0,1,n);
        Ttakny=Tny+(Tny-T)/3
        fel=Ttak-Ttakny;
        Ttak=Ttakny;
        iter=iter+1;
end
disp([Tny fel iter])

%% b-del

n=1;
T=litenTrapets(@fEXS64, 0,100,n);
n=2;
Tny=litenTrapets(@fEXS64, 0,100,n);
Ttak=Tny+(Tny-T)/3;
iter=1;
while abs(fel)>Ttak*1e-5
        n=n*2;
        T=Tny;
        Tny=litenTrapets(@fEXS64, 0,100,n);
        Ttakny100=Tny+(Tny-T)/3;
        fel=Ttak-Ttakny100;
        Ttak=Ttakny100;
        iter=iter+1;
end
disp([Tny fel iter])
