p1l1 function

Numerically Stable p1l1(t) = (t+1)*log(1+t) - t