void daccum(double*,long,long,double*,long,long,long*,long,long*);
void dcopy(long,double[],long,double[],long);
void dhfti(double*,long,long,long,double*,long,long,double,long*,double[],double[],
      long[]);
double dmpval(double[],long,double);