  syms a b
  f=(a+b)^2
  expand(f)
