Hi Fabrizio, you wrote in part:
> I can't find a Maxima command to triangularize
> a matrix (possibly non-square).
Perhaps triangularize(M) is what you need. ? triangularize
claims triangularize can handle non-square matrices.
Looking at the code (src/matrix.lisp and src/mat.lisp),
$triangularize calls triang, which calls tfgeli, which
says it is a "TWO-STEP FRACTION-FREE GAUSSIAN ELIMINATION ROUTINE".
I don't know anything more about this.
Hope this helps,
Robert Dodier
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com