Хотя бы minimax аппроксимацию возьмите и вообще exp/log нужно кусочно аппоксипировать и на бесконечности и около нуля там будую скорее дробнорациональные аппроксимации
С log(1+x)/x при х-> -1 тоже чепуха. Может они там модуль потеряли в условии?
А с другой стороны, их результат со сдвигом на единицу - для log(2 + x)/(1 + x) качественно отличается от предыдущего. Хотя там для IEEE754 разницы быть не должно. Так?