Make determinant work for singular matrices

Note that the current implementation of the determinant is rather
subobtimal.  While it is rather efficient, it is neither the most
numerically stable way to compute the determinant, nor does it work
for integers.  But I think it is much better than no implementation at
all.

[[Imported from SVN: r5549]]