syms x y
p = 1 + x + y + x*y
tic
for i=1:10000
coeffs(p);
end
toc
tic
for i=1:10000
q=expand(p*p);
end
toc

p := 1 + x + y + x*y
for i from 1 to 10000 do coeff(p);monomials(p); end_for
for i from 1 to 10000 do q:=expand(p*p); end_for