import numpy as np

print('------ 3.8 ---------------')
print("Système de 4 équations à 4 inconnues.")
A = np.array([[4, -2, 4, 4], [2, -5, 7, -9], [3, -2, 3, 4], [7, 3, -3, 8]], np.float64)
b = np.array([12, 13, 8, 6], np.float64)
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