Implement characteristic polynomial

Implement characteristic polynomial using `berkowitz`'s algorithm. Only the coefficients of the characteristic polynomial are returned, not the characteristic polynomial itself.
That's a nice addition. It's not a part of the CSymPy namespace yet. Thanks for the suggestion.