%Trapetsregeln med två extrapolationer
int=0.8820813907624;%13 korr decimaler
a=0;b=2;
n=10;h=(b-a)/n;
x=[a:h:b];
f=fexp(x);
T(1)=h*(sum(f)-(f(1)+f(n+1))/2);
S(1)=0;
errS(1)=0;
k=2;
for n=[20 40 80 160]
    h=h/2;
    x=[a:h:b];
    f=fexp(x);
    T(k)=h*(sum(f)-(f(1)+f(n+1))/2);
    errT(k)=T(k)-int;
    korr(k)=(T(k)-T(k-1))/3;
    S(k)=T(k)+korr(k);
    errS(k)=S(k)-int;
    k=k+1;
end
disp('           T(h)           errT     S(h)=T(h)+(T(h)-T(2*h))/3       errS')
[T' errT' S' errS']
errS(2:4)./errS(3:5)
pause
disp('           T(h)           errT        korr=(T(h)-T(2h))/3')
[T' errT' korr']

    
        