Accurate polynomial interpolations of special functions
  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. 
C. Semay