I am having a problem working in sagemath with the following series. In Maple, I could define $$n(d):= \sum_{k=1}^{d} schur((k,1^{d-k}))\prod_{\square \in (k,1^{d-k})}G(c(\square)h) \tag{*}$$ where schur is function is usual schur function I have created using character formula. $G(h)$ is a series in $h$. The $c(\square)$ is the content of the young tableaux. For example in case of $d=3$ it can be read as follows $$schur([1,1,1])G(h)G(2h) - schur([2,1])G(h)G(-h)+ schur([3])G(2h)G(h)$$
I wanted to do the similar thing in sagemath with jack polynomials instead of schur polynomial. Hence I defined
~~~
Sym = SymmetricFunctions(QQ)
JJ = Sym.jack(t=1).J()
s = Sym.schur() ~~~
I can see that I have defined the symmetric function over field of rational, even when I define a new variable in sage var('h') I cannot multiply $JJ[2,1] $ and $h$ it gives error hence I cannot compute $*$ with jack polynomial in sagemath. I am sure there is a way around. Please let me know