from sympy import *
init_printing()
x, y, z = symbols('x y z')
expr = exp(x)/y + z**2
expr
expand((x + y)**3)
factor(x**2 + 2*x*y + y**2)
sympify((x + 2*x + x**3)*4 + 3 + x + 4*x)
expr = x * y + 2/y
solve(expr, x)
# 1/x を x について無限大の極限をとる
limit(1/x, x, oo)
diff(sin(x), x)
integrate(cos(x), x)