First time here? Check out the FAQ!

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?

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.