The triangularize function isn't idempotent. Is there a good reason for this? (%i13) matrix([2,4],[0,6]); (%o13) matrix([2,4],[0,6]) (%i14) triangularize(%); (%o14) matrix([2,4],[0,12]) (%i15) triangularize(%); (%o15) matrix([2,4],[0,24]) Barton