% INTEXP0 Integration relative to the weight function w(t)=exp(1/t)
%
f1='%5.0f %21.13e\n';
fprintf('\n')
disp('    n      n-point Gauss')
load -ascii abexp0;
ab=abexp0; 
for n=2:2:12
  xw=gauss(n,ab);
  int=sum(xw(:,2).*log(1+xw(:,1)));
  fprintf(f1,n,int)
end
fprintf('\n')
abL=r_laguerre(200);
for N=95:105
  xwL=gauss(N,abL);
  intL=exp(-1)*sum(xwL(:,2).*(1+xwL(:,1)).^(-2).*log(1+...
    (1+xwL(:,1)).^(-1)));
  fprintf(f1,N,intL)
end