Apply this function over 1:J to calculate each portion of the integral we need for the upper bound.
performIntegralUpper1(j, muVec, sigMat, lBounds1, uBounds1, lBounds2, uBounds2)
Apply over this integer, the element that will be the largest in magnitude.
Mean vector of test statistics under the alternative (assuming it's MVN).
Covariance matrix of test statistics under the alternative (assuming it's MVN).
A 3J-2 vector of lower bounds for the first integral (see paper), bounds will be longer than for performIntegralLower1.
A 3J-2 vector of upper bounds for the second integral (see paper).
A 3J-2 vector of lower bounds for the first integral (see paper).
A 3J-2 vector of upper bounds for the second integral (see paper).
The value of the integration.