double *vandermonde_approx_1d_coef ( int n, int m, double x[], double y[] );
double *vandermonde_approx_1d_matrix ( int n, int m, double x[] );