|   | 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.
 
                
                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.