Hello Maple insiders, I'm still learning to get along with Maple (version 10.01). Currently, I would like to have a small procedure which creates a random matrix of arbitrary dimension with complex floats. However, the generated matix, say A, has to fulfill the following conditions (because it should be a valid physical density operator): 1.) A is hermitian 2.) A is positive semidefinite, together with 1.) this should mean that A has only real and nonnegative eigenvalues 3.) A should have trace equal 1 With the randomly generated matrices I would like to cover (reasonably uniformly) the set of all possible matrices with the above constraints.