Implement characteristic polynomial

Implement characteristic polynomial using `berkowitz`'s algorithm. Only the coefficients of the characteristic polynomial are returned, not the characteristic polynomial itself.
I defined this function so that we can use this to find eigenvalues, but we can't do it with CSymPy for the time being as we don't have a solver yet. I'll remove the definition and open an issue for implementing a solver.