Computing the hypergeometric function of a matrix argument