Add support for Berkowitz algorithm in SymEngine Python wrapper

[thurnerp]

The Berkowitz algorithm is a well-known and widely-used method for computing the determinant of a matrix. It would be great if this algorithm could be made accessible through the Python wrapper of SymEngine.

Currently, the Berkowitz algorithm is implemented in the C++ core of SymEngine, but it is not exposed in the Python wrapper. This makes it difficult for Python users to take advantage of this functionality without writing C code or using a low-level library like ctypes.

I would like to request that the berkowitz function be added to the Python wrapper of SymEngine, so that users can more easily compute determinants of symbolic matrices in their Python scripts.

Thank you for your consideration of this request!