lgamma1p function

Accurate log(gamma(a+1))