double lagrange_basis_function_1d ( int mx, double xd[], int i, double xi );
double *lagrange_interp_2d ( int mx, int my, double xd_1d[], double yd_1d[], 
  double zd[], int ni, double xi[], double yi[] );