Covariance matrix inequalities for functions of Beta random variables

Based on Jacobi polynomial series expansion and some innovative definitions of high-order differential matrix, we derive lower and upper bounds on covariance matrix for multivariate functions of Beta random variables in the sense of Loewner ordering for matrices. Additionally, corresponding univariate results obtained conveniently along the line of our arguments.