gamln1 function

Compute log( Gamma(x+1) ) Accurately in [-0.2, 1.25]