import sympy as sym
x=[Link]('x')
y=[Link]('y')
#differentiation
print([Link](x**5))
#integration
print([Link]((3*x**4+4*x**3+4*x),(x,1,3)))
#solve expression
sol=[Link](x**2+x-12,x)
print(sol)
sol=[Link]((x+5*y-2,-3*x+6*y-15),(x,y))
print(sol)
#expand
expr = (x + 1) ** 2
print([Link](expr))
#simplify
print([Link]((x + x * y) / x))
#factors
f=x**4-3*x**2+1
print([Link](f))
print([Link](x & y))
#matrix
m1 = [Link]([[1,2,3],[4,5,6],[7,8,9]])
m2 = [Link]([[9,8,7],[6,5,4],[3,2,1]])
print(m1)
#addition of two matrix
[Link](m1+m2)
#multiplication
[Link](m1*m2)
#eigen vectors
eg=[Link](m1)
print(eg)