This function takes a SymEngine object and return its expanded form.
expand(x)
x
A Basic/VecBasic/DenseMatrix S4 object.
Same type as input.
expr <- S(~ (x + y) ^ 3) expand(expr)