F.9.5.3 The lgamma functions
lgamma(1) returns +0.
lgamma(2) returns +0.
lgamma(x) returns +(inf) and raises the “divide-by-zero” floating-point exception for x a negative integer or zero.
lgamma(-(inf)) returns +(inf).
lgamma(+(inf)) returns +(inf).