%Trapetsregeln med extrapolation
n=10;h=2/n;x=h*[0:n];
fx=fq(x);
T=h*(sum(fx)-(fx(1)+fx(n+1))/2)
n=2*n;h=h/2;x=h*[0:n];
fx=fq(x);
Tt=h*(sum(fx)-(fx(1)+fx(n+1))/2)
S=Tt+(Tt-T)/3
n=2*n;h=h/2;x=h*[0:n];
fx=fq(x);
Ttt=h*(sum(fx)-(fx(1)+fx(n+1))/2)
St=Ttt+(Ttt-Tt)/3
richschema=[T 0;
            Tt S;
            Ttt St]
        