How can I get a random positive definite (or positive semi definite) matrix of given order? I know that there is a inbuilt function random_matrix with an additional feature algorithm to set with it. But there we can get some special matrices like 'echelon_form', 'orthogonal' ', 'echelonizable', 'diagonalizable'.... but positive definite command is not in built there. Some other important matrix classes are not there. So how to obtain random matrix of those classes?