Revision history [back]

The random_matrix command seems to be the ticket if you use an exact base ring

A=random_matrix(QQ, 4, 4, algorithm='echelonizable', rank=4, upper_bound=60); A

But this is not implemented for RR. Would floating-point issues even allow this?