In [1]: f = (x + y)**2

   In [2]: expand(f)
   Out[2]: 
            2    2
   2*x*y + x  + y 

   In [1]: diff( sin(x)*cos(x), x )
   Out[1]: 
      2         2   
   cos (x) - sin (x)

   In [2]: integrate( cos(x)**2, x )
   Out[2]: 
   x   cos(x)⋅sin(x)
   ─ + ─────────────
   2         2