The combinatorial species of rooted trees can be defined recursively via $A = X\cdot E(A)$, where $X$ denoted the singleton species and $E$ the set species. How can I make the following code work?
sage: X = species.SingletonSpecies()
sage: E = species.SetSpecies()
sage: A = CombinatorialSpecies()
sage: A.define(X*E(A))
sage: A.generating_series()