clear, clc;

f1=@(x) exp(x)./(1+2*x.^3);
a=0;b=3;
L=b-a;
N=1;
di=inf;
format long

integralen=1;
while abs(di)>1e-4*integralen
    integralenGammal=integralen;
    h=L/N;
    x=a:h:b;
    f=f1(x);
    integralen=h*(sum(f)-(f(1)+f(end))/2);
    disp(integralen)
    di=integralen-integralenGammal;
    N=N*2;
end
