%
% sI101AA  The integral I_{1,[0,1]}(a,a) in
%    variable-precision arithmetic
%
function y=sI101aa(dig,a)
digits(dig);
if abs(subs(a))>=1
  error('a not in range')
  return
end
p=feval(symengine,'psi',1-a,1);
psi=vpa(p);
y=a*psi+1;
