Support computing determinants

[Jacalz]

We should add support for computing the determinant of square matrices. Using a brute force algorithm (just recursively expanding along a row or column until we reach a 2x2 matrix) should be a good enough start but will likely be too slow in the long run.