Advertisement
mikhail_dvorkin

SymPy Example 2

May 5th, 2016
263
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 1.08 KB | None | 0 0
  1. import sympy
  2.  
  3. x = sympy.Symbol("x")
  4. y = sympy.Symbol("y")
  5.  
  6. m = sympy.Matrix([[sympy.cos(x), sympy.sin(x)], [-sympy.sin(x), sympy.cos(x)]])
  7. m = sympy.Matrix([[1, x], [y, 1]])
  8. m = sympy.Matrix([[0, 1], [-1, 0]])
  9. print(m)
  10. print(repr(m))
  11. print(sympy.simplify(m.det()))
  12. print(m.inv())
  13. print(m ** -1)
  14. f = sympy.Matrix(m.inv()).row(0)[0]
  15. print(sympy.simplify(f))
  16. print(m ** 2)
  17. print(sympy.simplify(m ** 4))
  18. print(sympy.simplify(m.eigenvals()))
  19. print(sympy.simplify(m.eigenvects()))
  20. p, d = m.diagonalize()
  21. print(m.charpoly())
  22. print(p)
  23. print(d)
  24. print(sympy.simplify(p * d * p.inv()))
  25.  
  26. i = sympy.solve(x**4 - 1, x)[3]
  27. print(type(sympy.solve(x**4 - 1, x)[3]))
  28.  
  29. e = sympy.Eq((x-2)**3, x - y)
  30. print(sympy.solve(e, x))
  31. print(sympy.solve(x ** 4 - 1, x))
  32. print(sympy.factor(x ** 4 - 3 * x**2 + 1, modulus=11))
  33. print(sympy.solve(x * x + x ** 3 < 3))
  34. print(sympy.solve(x**3 + x**2 - 3, x))
  35. print(sympy.satisfiable((x | y) & (x | ~y) & y))
  36.  
  37. f = sympy.Symbol("f")
  38. print(sympy.latex(f(x, y).diff(*[x, 3])))
  39.  
  40. f = x ** 4 + x ** 2 + sympy.exp(x ** 4)
  41. g = f.subs(x ** 2, y)
  42. print(g.diff(y))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement