solve(x^6-9*x^5+10*x^4+18*x^3+49*x^2+135*x+84 == 0, x )
solve(x^3-5*x^2+8*x-40 == 0,x)
find_root(cos(x)==x, 0, pi/2)
integral(x*sin(x),x)
integral(e^(1/x^2),x,-2,4)
diff(sqrt(x^2*sin(x)+2),x)
diff((x*e^(cos(x)))^3,x)
x,y=var('x,y')
f=x^2+sin(x*y)-y
f.diff(y)
x=var('x')
u=function('u')(x)
rr=diff(u,x,2)-x*diff(u,x)+u
desolve(rr,[u,x],contrib_ode=True)
A=matrix([[1,-4, 7],[14, 6,  8],[-4, -1, 2]])
B=matrix([[2, 3, 9],[-1, 17, 6],[ 0, 5, -7]])
A*B
det(A)
B^-1