title: How to find the eigenvalues and eigenvectors of a matrix with SymPy?
body: <p>I want to calculate the eigenvectors <strong>x</strong> from a system <em>A</em> by using this: <em>A</em> <strong>x</strong> = <em>λ</em> <strong>x</strong></p>
  
  <p>The problem is that I don't know how to solve the eigenvalues by using SymPy.
  Here is my code. I want to get some values for x1 and x2 from matrix <code>A</code></p>
  
  <pre><code>from sympy import *
  x1, x2, Lambda = symbols('x1 x2 Lambda')
  I = eye(2)
  A = Matrix([[0, 2], [1, -3]])
  equation = Eq(det(Lambda*I-A), 0)
  D = solve(equation)
  print([N(element, 4) for element in D]) # Eigenvalus in decimal form
  print(pretty(D)) # Eigenvalues in exact form
  
  X = Matrix([[x1], [x2]]) # Eigenvectors
  T = A*X - D[0]*X # The Ax = %Lambda X with the first %Lambda = D[0]
  print(pretty(solve(T, x1, x2)))
  </code></pre>