gammaln() 在处理上溢和下溢方面要比 log(gamma()) 更好;
log1p() (1 plus)对于较小的 x 值时,log1p(x) 要比 log(1+x) 更为精确;
log1p(-.99) ⇒ log(.01)
>> eps ans = 2.2204e-16