julia> using AbstractAlgebra
julia> S, (a,b,c) = PolynomialRing(QQ, ["a","b","c"]);
julia> T, (x,y) = PolynomialRing(S, ["x","y"]);
julia> p = a*x + b*x^2*y + c*x*y
b*x^2*y + c*x*y + a*x
julia> p+p
2*b*x^2*y + 2*c*x*y + 2*a*x
4 Likes