from sympy import *
x, y = symbols('x y')
expr = expand((x+y)**6)
expr
diff(expr, x)
diff(expr, y)
diff(expr, x, 2)
diff(expr, x, x)
diff(expr, x, x, y, x)
diff(expr, y, x, x, x)
diff(expr, y, x, 2, y, x)
diff(sin(x),x)
diff(sin(x),x, 4)
integrate(x**2, (x, 0, 1))
integrate(E ** (-x**2), (x, -oo, oo))
integrate(E ** (-x**3), (x, -oo, oo))
4 * integrate(sqrt(1 - x**2), (x, 0, 1))