| 1 | initial version |
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.