> As someone who does research […] mainly Wolfram Mathematica […] it took me ages to figure out how to implement this in SymPy, whereas I in Mathematica I can write this in a minute.

Hardly surprising, is it? As someone with more experience in Python, it would likely be the opposite for me.