function y=sgam_upinc(dig,a,x)
%SGAM_UPINC Upper incomplete gamma function.
%   Y=SGAM_UPINC(DIG,A,X) computes the integral Y from X to 
%   Inf of t^{A-1}*exp(-t) in DIG-digit arithmetic for real A
%   and X>0.

digits(dig); sa=vpa(a); sx=vpa(x);
y=feval(symengine,'igamma',sa,sx);
