performIntegralUpper2: Apply this function over all m, j not equal (order matters)...
In DBpower: Finite Sample Power Calculations for Detection Boundary Tests

Apply this function over all m, j not equal (order matters) to calculate each portion of the integral we need for the lower bound for calc_b2.

performIntegralUpper2(
  x,
  muVec,
  sigMat,
  lBounds1,
  uBounds1,
  lBounds2,
  uBounds2,
  lBounds3,
  uBounds3,
  lBounds4,
  uBounds4
)

`x`	Apply over this 2*1 vector, the elements that will be the largest and second largest in magnitude.
`muVec`	Mean vector of test statistics under the alternative (assuming it's MVN).
`sigMat`	Covariance matrix of test statistics under the alternative (assuming it's MVN).
`lBounds1`	A 3J-2 vector of lower bounds for the first integral (see paper).
`uBounds1`	A 3J-2 vector of upper bounds for the second integral (see paper).
`lBounds2`	A 3J-2 vector of lower bounds for the first integral (see paper).
`uBounds2`	A J3J-2 vector of upper bounds for the second integral (see paper).
`lBounds3`	A 3J-2 vector of lower bounds for the third integral (see paper).
`uBounds3`	A 3J-2 vector of upper bounds for the third integral (see paper).
`lBounds4`	A 3J-2 vector of lower bounds for the fourth integral (see paper).
`uBounds4`	A 3J-2 vector of upper bounds for the fourth integral (see paper).