Abstract.The relationships between various numerical methods for obtaining polynomial approximations to the first derivative of a known function are investigated, and their computational advantages discussed. Optimum sequences of interpolation points are then selected with the objective of minimising the relative contribution of rounding errors to the total error, and geometric sequences, though non-optimal in this sense, are considered for computational reasons.