using Symata [a,b,c] {a,b,c} [ a c + d "cat" Expand((x+y)^2) ] f(x) f % list f % [a,b,c] x .% y f .% g(1,2) [Rule(a,b), a => b , a ⇒ b] [RuleDelayed(a,b), a .> b] [a, x^2, b^3, (a+b)^3] ./ ( x_^n_ => g([n],x))