plot of f(p)=p(1-exp(-(1-p)/p))/(1-p) for p\in (0,1)