Well, that may be because scipy and numpy are not intended for symbolic computation. If you want to do this, use a piece of Software that is intended for those things. SageMath would be an option.

The true values of Mathematica are it’s Solve-Algorithms and it’s Simplification-Algorithms. They did a heap fine job there.