default(realprecision,38);
default(seriesprecision,6);
gamma(2+I+x)
gamma(1+x)
gamma(-2+x)
lngamma(2+I+x)
lngamma(1+x)
lngamma(-2+x)
psi(2+I+x)
psi(1+x)
psi(-2+x)
gamma(-1/2+x)

psi(x)
gamma(x)
gamma(1000)
psi(2^400)
psi(-1.5)

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)
