# Series expansion
x = Symbol("x")
print((1/cos(x)).series(x, 0, 10))

# Integration
print(integrate(exp(-x**2)*erf(x), x))

# Differentiation
y = Symbol("y")
f = x**2/y + 2*x - ln(y)
print(diff(f, x))