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)