default(realprecision,38);
default(seriesprecision,6);
gamma(2+I+x)
gamma(1+x)
gamma(2+x)
gamma(-2+x)
lngamma(2+I+x)
lngamma(1+x)
lngamma(-2+x)
lngamma(2+x)
lngamma(1)
lngamma(1.)
psi(2+I+x)
psi(1+x)
psi(-2+x)
gamma(-1/2+x)
gamma(1+a*x+O(x^2))
lngamma(1+a*x+O(x^2))
gamma(O(x))
lngamma(O(x))
psi(O(x))
gamma(0.+1/2*x+O(x^2))

psi(x)
gamma(x)
gamma(1000)
psi(2^400)
psi(-1.5)
psi(x+O(x^2))
real(psi((3/2+0.E-38*I)+x/2+O(x^2)))
gamma(x+O(x^2))
gamma(-1-10^-16)
gamma(10^-16)
lngamma(-1-10^-16)
lngamma(10^-16)

binomial(2^64+1,2)
binomial(1001.,1000)

lngamma(2^301)
lngamma(2^301 + 2*I)
lngamma(2^64 + 2*I)
lngamma(-200.5)
gammah(400)
gammah(400+I)
default(realprecision,115);
lngamma(2^301 + 2*I)
lngamma(2^301 + 200*I)
default(realprecision,154);
lngamma(2^301 + 2000*I)

gamma(1+10^-30)
gamma(-1e-6)
