plot[(-13.8629-log(1-x))/log(1-x), {x, 0, 1}]