Mike Zabrocki writes:
> sage: a = var('a')
> sage: M = MacdonaldPolynomialsH(QQ[a])
> sage: s = SFASchur(M.base_ring())
> sage: q = M.q
> sage: s[2,1].plethysm(s[1]/(1-q))
> (((-9)*q^2+9*q)/(9*q^6+(-27)*q^5+27*q^4+(-18)*q^3+27*q^2+(-27)*q+9))*s[1,
> 1, 1] +
> (((-9)*q^3+9*q^2+(-9)*q+9)/(9*q^6+(-
I want to use symmetric functions with an extra parameter in the
coefficient ring.
Can anyone tell me what I am doing wrong?
sage: a = var('a')
sage: M = MacdonaldPolynomialsH(QQ[a])
sage: s = SFASchur(M.base_ring())
sage: q = M.q
sage: s[2,1].plethysm(s[1]/(1-q))
(((-9)*q^2+9*q)/(9*q^6+(-27)*q^
