function [xg, wg, xk, wgn, wk] = gaukrondata
% Gausspunkter och vikter
   xg=[0.405845151377397 0.741531185599394 0.949107912342758];  
   xg=[-fliplr(xg)  0  xg];  
   wg=[0.381830050505119 0.279705391489277 0.129484966168870]; 
   wg=[fliplr(wg)  0.417959183673469  wg]; 
 % Kronrodpunkter och nya vikter
   xk=[0.207784955007898 0.586087235467691 0.864864423359769 0.991455371120813];
   xk=[-fliplr(xk)  xk]; 
   wgn=[0.190350578064785 0.140653259715525 0.063092092629979];
   wgn=[fliplr(wgn)  0.209482141084728  wgn];
   wk=[0.204432940075298 0.169004726639267 0.104790010322250 0.022935322010529];
   wk=[fliplr(wk)  wk];


