function [Igk,absfel,x15,y15]=g7k15p(F,a,b);
% Gauss-Kronrods 15-punktsmetod, integrand F(x), intervall [a,b]
% Förutom integralvärde och felskattning sparas de 15 punkterna
% Gausspunkter och vikter
   xg=[0.405845151377397  0.741531185599394  0.949107912342758];  
   xg=[-fliplr(xg)  0  xg];  
   w= [0.381830050505119  0.279705391489277  0.129484966168870]; 
   w= [fliplr(w)  0.417959183673469  w]; 
% Kronrodpunkter och nya vikter
   xk=[0.207784955007898 
       0.586087235467691 
       0.864864423359769 
       0.991455371120813]';   xk=[-fliplr(xk)  xk]; 
   u= [0.190350578064785  0.140653259715525  0.063092092629979];
   u= [fliplr(u)  0.209482141084728  u];
   v= [0.204432940075298 
       0.169004726639267 
       0.104790010322250 
       0.022935322010529]';   v=[fliplr(v)  v];
% Variabelsubstitution till standardintervallet [-1,1]
   hh=(b-a)/2;  abm=(a+b)/2;
   tg=hh*xg+abm;  fg=feval(F,tg); 
   tk=hh*xk+abm;  fk=feval(F,tk);  
   x15=[tg tk];   y15=[fg fk];   
   G7=hh*sum(w.*fg);                  % Sjupunktsformeln
   K15=hh*(sum(u.*fg)+sum(v.*fk));    % Femtonpunktsformeln 
   d=abs(K15-G7);                     % Avvikelse
   fel=10*d*sqrt(d/abs(K15));
   Igk=K15;                           % Integralapproximation
   absfel=max([eps*abs(K15)  fel]);   % Felskattning 
  

