>>> for u in [1,2,2000]:
      print('Les valeurs propres de la Jacobienne en ('+str(u)+','+str(u)+') sont : ')
      print(sol.solve(carac.subs({x:u,y:u})))
Les valeurs propres de la Jacobienne en (1,1) sont : 
[2, 2000]
Les valeurs propres de la Jacobienne en (2,2) sont : 
[1/2, 1000]
Les valeurs propres de la Jacobienne en (2000,2000) sont : 
[1/2000, 1/1000]