import numpy as np

print('------ 3.7 ---------------')
print("Système de 3 équations à 3 inconnues.")
A = np.array([[1, 1, 1], [1, 2, 3], [3, 5, -6]], np.float64)
Ainv = np.linalg.inv(A) # Matrice inverse de la matrice A
b = np.array([-3, -4, 15], np.float64)
v = np.dot(Ainv, b) # produit Ainv*b, donne la solution du sytème.
print("v=", v)  # Affichage de la solution
print("A*v=", np.dot(A, v)) # Vérification que l'on retrouve b
print("Meilleure manière de résoudre.")
v_sol = np.linalg.solve(A, b) # Résolution du système d'équations.
print("v_sol=", v_sol)  # Affichage de la solution
print("A*v_sol=", np.dot(A, v_sol)) # Vérification que l'on retrouve b