I would like to compute the determinant of a matrix involving a symbolic variable:

```
using Symbolics
I = # Identity matrix
A = # Some non-trivial matrix
@variables z
aux = I - z * A
det(aux)
```

This webpage claims that computing determinants is possible, but I cannot find any documentation of that.

Does `Symbolics.jl`

currently offer this feature, or do I have to resort to the `SymPy`

wrapper for this?