double *newton_coef_1d ( int nd, double xd[], double yd[] );
double *newton_value_1d ( int nd, double xd[], double cd[], int ni, double xi[] );