Implement characterstic polynomial

Implement characteristic polynomial using `berkowitz`'s algorithm. Only the coefficients of the characteristic polynomial are returned, not the characteristic polynomial itself.
Can you document how to obtain the polynomial from `B`? Is the algorithm that you need to multiply out all the elements of the `(n, 1)` vector `B`?