% qrombg, Integralberäkning med Rombergs metod
  clear 
  flops(0)
  n=40,  h=2/n; x=h*(0:n); f=fq(x);  
  T=h*(sum(f)-(f(1)+f(n+1))/2);      % trapetsregeln 
  n=2*n, h=h/2; x=h*(0:n); f=fq(x);  % halverat steg
  Tt=h*(sum(f)-(f(1)+f(n+1))/2);     % trapetsregeln  
  S=Tt+(Tt-T)/3;                     % extrapolera till simpsonvärde S
  n=2*n, h=h/2; x=h*(0:n); f=fq(x);  % halverat steg igen
  Ttt=h*(sum(f)-(f(1)+f(n+1))/2);    % trapetsregeln 
  St=Ttt+(Ttt-Tt)/3;                 % extrapolera fram nytt S
  d=St-S;  B=St+d/15;                % extrapolera på S-värdena 
  format long
  richschema=[T    0    0
              Tt   S    0
              Ttt  St   B]
  rmbgflops=flops




