void dc2fit(double[],double[],double[],long,double[],long,double*,long,double[],
      double[],double*,long*);
double dhint(double,long,long,double[],double[],double[]);