Mattias Engdegård <mattiase@acm.org> writes:
> No doubt there are better ways of computing the determinant for
> non-numeric matrices. In particular, division shouldn't be necessary
> at all, and the current algorithm results in some untidy
> expressions. For example, the determinant of
>
> [a 0 0 1]
> [0 b 0 0]
> [0 0 c 0]
> [1 0 0 d]
>
> is given as abc*(d-1/a), falsely giving the impression of being
> undefined for a=0, rather than bc(ad-1).
I see. Yes, using a fraction-free algorithm would be good. But for the
moment, getting bc(ad-1) as the result is one `j M' away.