title: Solving system of linear equations with square roots
body: <p>Say I have a system of linear equations with a square root</p>
  
  <pre><code>  1  1  |  1  
    (1/2 + sqrt(5) / 2)  (1/2 - sqrt(5) / 2)  | 1
  </code></pre>
  
  <p>Using np.linalg.solve to solve this system of equations I would normally do</p>
  
  <pre><code>vars = [[1, 1], [1/2 + sqrt(5)/2, -sqrt(5)/2 + 1/2]]
  outcomes = [1, 1]
  solution = np.linalg.solve(vars, outcomes) 
  #solution has to be only whole numbers, no crazy decimals. Preferably in the following form
  [ sqrt(x), sqrt(y) ]
  </code></pre>
  
  <p>However this returns an error because it doesn't know what to do with the sqrt(). How can I solve this system of equations with square roots and get full numbers back, so no decimals?</p>