% SI01INF Symbolic analogue of I01inf.m
%
a=1/6;
dig=28;
%dig=37;
digits(dig);
psi1=feval(symengine,'psi',a,1);
psi1=vpa(psi1); exact=vpa(a*psi1-1)
N=105; 
%N=155;
u=8.4;
%u=11.1;
sab=sr_jacobi01(dig,N);
%for n=80:5:N
%for n=135:5:N
for n=N
  sxw=sgauss(dig,n,sab(1:n,:));
  s=0;
  for k=1:n
    sx=vpa(1+u*sxw(k,1)/a);
    w=swofy0(dig,-sx*exp(-sx));
    sw=(-w/sx)^(vpa(a));
    s=vpa(s+sxw(k,2)*sw);
  end
  int=vpa((u*s+exp(-vpa(a+u)))/a);
  err=vpa(abs(int-exact));
  [int err]'
end