import org.apache.commons.math3.analysis.solvers.UnivariateSolverUtils; //导入依赖的package包/类
/**
 * Prepare for computation.
 * Subclasses must call this method if they override any of the
 * {@code solve} methods.
 *
 * @param maxEval Maximum number of evaluations.
 * @param f the integrand function
 * @param lower the min bound for the interval
 * @param upper the upper bound for the interval
 * @throws NullArgumentException if {@code f} is {@code null}.
 * @throws MathIllegalArgumentException if {@code min >= max}.
 */
protected void setup(final int maxEval,
                     final UnivariateFunction f,
                     final double lower, final double upper)
    throws NullArgumentException, MathIllegalArgumentException {

    // Checks.
    MathUtils.checkNotNull(f);
    UnivariateSolverUtils.verifyInterval(lower, upper);

    // Reset.
    min = lower;
    max = upper;
    function = f;
    evaluations = evaluations.withMaximalCount(maxEval).withStart(0);
    count       = count.withStart(0);

}
 

import org.apache.commons.math3.analysis.solvers.UnivariateSolverUtils; //导入依赖的package包/类
/**
 * Prepare for computation.
 * Subclasses must call this method if they override any of the
 * {@code solve} methods.
 *
 * @param maxEval Maximum number of evaluations.
 * @param f the integrand function
 * @param lower the min bound for the interval
 * @param upper the upper bound for the interval
 * @throws NullArgumentException if {@code f} is {@code null}.
 * @throws MathIllegalArgumentException if {@code min >= max}.
 */
protected void setup(final int maxEval,
                     final UnivariateFunction f,
                     final double lower, final double upper)
    throws NullArgumentException, MathIllegalArgumentException {

    // Checks.
    MathUtils.checkNotNull(f);
    UnivariateSolverUtils.verifyInterval(lower, upper);

    // Reset.
    min = lower;
    max = upper;
    function = f;
    evaluations.setMaximalCount(maxEval);
    evaluations.resetCount();
    iterations.resetCount();

}
 

import org.apache.commons.math3.analysis.solvers.UnivariateSolverUtils; //导入依赖的package包/类
/**
 * Solves the method-of-moments estimate numerically. We use a cache
 * on the same observed instances in the sample as q is constant and
 * min/max are chosen conservatively.
 * 
 */
private static double getMethodOfMomentsEstimate(int nj, double q, double min, double max, 
	HashMap<Integer, Double> solveCache) {
	if( solveCache.containsKey(nj) )
		return solveCache.get(nj);
	
	double est = UnivariateSolverUtils
		.solve(new MethodOfMomentsFunction(nj, q), min, max, 1e-9);
	
	if( solveCache.size()<MAX_SOLVE_CACHE_SIZE )
		solveCache.put(nj, est);
	
	return est;
}