int triangle_num ( int n );
double *vandermonde_interp_2d_matrix ( int n, int m, double x[], double y[] );