  Provided a special function of one variable and some of its derivatives can
be accurately computed over a finite range, a method is presented to build a
series of polynomial approximations of the function with a defined relative
error over the whole range. This method is easy to implement and makes possible
fast computation of special functions.
