adapt_integrate()

int adapt_integrate	(	unsigned	fdim,
		integrand	f,
		void *	fdata,
		unsigned	dim,
		const double *	xmin,
		const double *	xmax,
		unsigned	maxEval,
		double	reqAbsError,
		double	reqRelError,
		double *	val,
		double *	err
	)

{
     int ret;
     fv_data d;

     if (fdim == 0) return SUCCESS; /* nothing to do */     
     
     d.f = f; d.fdata = fdata;
     d.fval1 = (double *) malloc(sizeof(double) * fdim);
     if (!d.fval1) {
          unsigned i;
          for (i = 0; i < fdim; ++i) {
               val[i] = 0;
               err[i] = HUGE_VAL; 
          }
          return -2; /* ERROR */
     }
     ret = integrate(fdim, fv, &d, dim, xmin, xmax, 
                     maxEval, reqAbsError, reqRelError, val, err, 0);
     free(d.fval1);
     return ret;
}