I aim to write a multidimensional Taylor approximation using sympy, which
sympy
You can use expr.removeO() to remove the big O from an expression.
expr.removeO()
Oneliner: expr.series(x, 0, 3).removeO().series(y, 0, 3).removeO()
expr.series(x, 0, 3).removeO().series(y, 0, 3).removeO()