Implement characterstic polynomial

Implement characteristic polynomial using `berkowitz`'s algorithm. Only the coefficients of the characteristic polynomial are returned, not the characteristic polynomial itself.
For the moment, this can't be implemented as we don't have a solver yet.