clear all, clc;

fel=inf;
a=0;
b=3;
delta=b-a;
n=1;
h=delta/n;
x=a:h:b;
I=exp(x)./(1+2*x.^3);
kanter=sum(I)/2;
T=h*(sum(I)-kanter);
while abs(fel)>1e-4
    n=2*n;
    h=delta/n;
    x=a:h:b;
    I=exp(x)./(1+2*x.^3);
    Told=T;
    T=h*(sum(I)-kanter);
    fel=T-Told;
end
T
n
n=1;
h=delta/n;
x=a:h:b;
I=exp(x)./(1+2*x.^3);
T=h*(sum(I)-kanter);
n=2*n;
h=delta/n;
x=a:h:b;
I=exp(x)./(1+2*x.^3);
Told=T;
T=h*(sum(I)-kanter);
Ttak=T+(T-Told)/3;
fel=inf;
while abs(fel)>1e-4
    n=2*n;
    h=delta/n;
    x=a:h:b;
    I=exp(x)./(1+2*x.^3);
    Told=T;
    Ttakold=Ttak;
    T=h*(sum(I)-kanter);
    Ttak=T+(T-Told)/3;
    fel=Ttak-Ttakold;
end
Ttak
n
