%Trapetsregeln med en extrapolation: I=int_0^2 exp(-x^2)dx
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);
errT(1)=T(1)-int;
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;
    k=k+1;
end
[T' errT']
errT(1:4)./errT(2:5)

    
        