In this paper, we propose a matrix-variate generalization of the Gauss hypergeometric distribution and study several of its properties. We also derive probability density functions of the product of two independent random matrices when one of them is Gauss hypergeometric. These densities are expressed in terms of Appell’s first hypergeometric function F1 and Humbert’s confluent hypergeometric function Φ1of matrix arguments.