Table-driven implementation of the Expm1 function in IEEE floating-point arithmetic

Algorithms and implementation details for the function <italic>e</italic><supscrpt><italic>x</italic> - 1</supscrpt> in both single and double precision of IEEE 754 arithmetic are presented here. With a table of moderate size, the implementations need only working-precision arithmetic and are provably accurate to within 0.58 ulp.