double *vandermonde_coef_1d ( int nd, double xd[], double yd[] );
double *vandermonde_matrix_1d ( int nd, double xd[] );
double *vandermonde_value_1d ( int nd, double cd[], int ni, double xi[] );