/* nuderiv.h  non uniformly spaced derivative coefficients */
void nuderiv(int order, int npoint, int point, double x[], double c[]);