from symbolica import *

x = S('x')

# Symbolic derivative with chain and product rules applied automatically
e = E('exp(x^2) * sin(x)').derivative(x).expand()
print(e)  # ➜  2*x*exp(x^2)*sin(x) + cos(x)*exp(x^2)

from symbolica import *

# Factor a multivariate polynomial over the rationals
p = E('x^3 - y^3')
for factor, exp in p.factor():
    print(factor)
# ➜  x - y
# ➜  x^2 + x*y + y^2

from symbolica import *

x, y, f = S('x', 'y', 'f')

# Mix built-in and user-defined functions in one expression
e = E('cos(x)*3*x^2 + f(x)')
ev = e.evaluator({}, {(f, "f", [y]): y**2}, [x])
print(ev.evaluate([0.3]))  # ➜  0.3479408520639137