double *shepard_basis_1d ( int nd, double xd[], int k, double p, int ni, 
  double xi[] );
double *shepard_interp_1d ( int nd, double xd[], double yd[], double p, int ni, 
  double xi[] );