int triangle_num ( int n );
double *vandermonde_approx_2d_coef ( int n, int m, double x[], double y[], 
  double z[] );
double *vandermonde_approx_2d_matrix ( int n, int m, int tm, double x[], 
  double y[] );