Implement characteristic polynomial

Implement characteristic polynomial using `berkowitz`'s algorithm. Only the coefficients of the characteristic polynomial are returned, not the characteristic polynomial itself.
Consider adding a EyeMatrix class/function to the csympy namespace (inspired by numpy.eye) if it's not already there.