On Jun 19, 6:50 am, Pierre <pierre.guil...@gmail.com> wrote: > is there an easy way of creating a subdiagonal matrix, that is with 1 > just under the diagonal and 0's elsewhere ?
This is a hack (especially since I'm no mathematician and no expert in Sage or any of its components), but it seems to work. A formal proof would be an interesting read. Here's an example: n = 10 m = matrix(ZZ,n) m[range(n-1),range(n-1)] = identity_matrix(n-1) m_subd = m.LLL() Works with some large n I played with. Richard --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
