Distribution of the ratio of two Wishart matrices and cumulative probability evaluation by the holonomic gradient method

We study the distribution of the ratio of two central Wishart matrices with different covariance matrices. We first derive the density function of a particular matrix form of the ratio and show that its cumulative distribution function can be expressed in terms of the hypergeometric function 2F1 of a matrix argument. Then we apply the holonomic gradient method for numerical evaluation of the hypergeometric function. This approach enables us to compute the power function of Roy's maximum root test for testing the equality of two covariance matrices.