double *lagrange_approx_1d ( int m, int nd, double xd[], double yd[], 
  int ni, double xi[] );
double *lagrange_basis_1d ( int nd, double xd[], int ni, double xi[] );