% INTEXPMINFPINF Integration relative to the weight function w(t)=exp(-1/t^2-t^2)
%
f1='%5.0f %21.13e\n';
fprintf('\n')
disp('    n      n-point Gauss')
load -ascii abexpminfpinf;
ab=abexpminfpinf; 
for n=2:2:14
  xw=gauss(n,ab);
  int=sum(xw(:,2).*exp(-xw(:,1)).*cos(xw(:,1)));
  fprintf(f1,n,int)
end
fprintf('\n')
abl=r_hermite(1000);
for n=100:100:1000
  xwl=gauss(n,abl);
  intl=sum(xwl(:,2).*exp(-xwl(:,1)-1./xwl(:,1).^2)...
    .*cos(xwl(:,1)));
  fprintf(f1,n,intl)
end