But I would like to know which were the different procedures (integration methods, algebraic manipulations, etc.) that were used to arrive at the result. Perhaps SymPy does not offer this option, but do you know of another solution?
I’m not a SymPy user but you might start by working through
the tutorials and documentation to understand what it does
and how to use it. I think that might give you the understanding
to solve your problem/question. It is written in Python so you
could possibly look at the source or trace the operation to see
what is being run…?