from sympy import *
a = symbols('tau10:15', integer=True) #,cls=sy.Function)
sum(a)
x,y,z,pi = symbols('x y z pi')
expr = y*(x**2 + y * sin(x) + 3* exp(z) + pi*x)
expr = $\displaystyle y \left(\pi x + x^{2} + y \sin{\left(x \right)} + 3 e^{z}\right)$
expr * y = $\displaystyle y^{2} \left(\pi x + x^{2} + y \sin{\left(x \right)} + 3 e^{z}\right)$
expr - x² = $\displaystyle - x^{2} + y \left(\pi x + x^{2} + y \sin{\left(x \right)} + 3 e^{z}\right)$
expand(expr)
factor(expr)
expr.subs(x,2)
expr.subs({x:2,z:20})
expr.subs([(x,2),(z,20)])
expr.evalf(subs={x:2,z:20,y:2,pi:3},)
Eq(expr,1024)
Eq(2+2,4)
string = 'x/5 + e**y + z + 6**3 + 0.12 + mu'
sympify(string)
expr = sin(x)**2 + cos(x)**2 + ln(y)
expr = ((x + 2) * (x + 3) )
factor(expr)
expand(expr)
expr = (tan(x)*sec(x)*Pow(cos(x),2) * Rational(5,3))
expr
trigsimp(expr)
sec(x).rewrite(cos)
sec(x).rewrite(sin)
powsimp(x**z * x**y)
expand_power_exp(x**z * x**y)
binomial(10,x)
log(5,4)
factorial(y*x)
gamma(y)
beta(x,10)
hyper([2,8],[10],x)