% qg7k15, integral med Gauss-Kronrods formel
% Anropar funktionen g7k15
  clear, clf 
%      function f=fq(x);
%            f=sqrt(0.5+2*exp(-x).*sin(2*x.^2));
% ------------------------------------------------------------------
  t=0:0.02:2;  plot(t,fq(t)), hold on, drawnow
  a1=0.8; a2=1.5;
  plot([a1 a1],[0 fq(a1)],'--', [a2 a2],[0 fq(a2)],'--')
  xlabel('Antal punkter: 15+15+15')   
  flops(0) 
  [I1,fel1]=g7k15('fq',0, a1); fel1
  [I2,fel2]=g7k15('fq',a1,a2); fel2
  [I3,fel3]=g7k15('fq',a2, 2); fel3
  format long,  Ig7k15=I1+I2+I3
  format short, absfel=fel1+fel2+fel3,  GKflops=flops




