A function whose direct evaluation gives large roundoff error

Let's actually evaluate this function and see what we get ...

Numpy arrays

Direct evaluation of log(1+x)/x

Taylor polynomial approximation of log(1+x)/x and its truncation error