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(1)
lngamma(1.)
[psi(2+I+x,d)|d<-[0..4]]
[psi(1+x,d)|d<-[0..4]]
[psi(-2+x,d)|d<-[0..4]]
psi(-1/2,1)
psi(2^-1024,1)
gamma(-1/2+x)
gamma(1+a*x+O(x^2))
lngamma(1+a*x+O(x^2))
gamma(0.+1/2*x+O(x^2))

psi(x)
gamma(x)
gamma(1000)
[psi(2^400,d)|d<-[0..9]]
[psi(2^400*I,d)|d<-[0..9]]
[psi(-1.5,d)|d<-[0..9]]
[psi(x+O(x^2),d)|d<-[0..9]]
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(0)
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)

\p38
gamma(1/3)
gamma(2/3)
gamma(1/4)
gamma(3/4)
gamma(1/6)
gamma(5/6)
gamma(1/8)
gamma(3/8)
gamma(5/8)
gamma(7/8)
gamma(1/12)
gamma(5/12)
gamma(7/12)
gamma(11/12)
gamma(1/24)
gamma(5/24)
gamma(7/24)
gamma(11/24)
gamma(13/24)
gamma(17/24)
gamma(19/24)
gamma(23/24)
gamma(1/24+100)
gamma(1/24+99)
gamma(1/24-99)
gamma(1/24-100)
gamma(1/24+10^5)
gamma(1/24-10^5)

\\ ERRORS
lngamma(-2+x)
gamma(O(x))
lngamma(O(x))
psi(O(x))
psi(-1)
