% I1ETA1 The integral I_{1,[eta,1]}(a)
%    for a=1/6, eta=.01
%
f0='%8.0f %18.14f\n';
a=1/6; eta=.01;
N=80;
ab=r_jacobi01(N);
for n=10:10:N
  xw=gauss(n,ab);
  s=0;
  for k=1:n
    t=eta+(1-eta)*xw(k,1);
    y=-t*exp(-t);
    w1=wofy1(y);
    err=abs((w1*exp(w1)-y)/y);
    w=(-wofy1(y)/t)^a;
    s=s+xw(k,2)*w; 
  end
  s=(1-eta)*s;
  fprintf(f0,n,s)
end