double *vandermonde_interp_1d_coef ( int n, double x[], double y[] );
double *vandermonde_interp_1d_matrix ( int n, double x[] );