import numpy as np

print('------ 3.5 ---------------')
print("Système de 2 équations à 2 inconnues.")
A = np.array([[4, -2], [9, -5]], np.float64)
Ainv = np.linalg.inv(A) # Matrice inverse de la matrice A
b = np.array([11, 23], 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('----- 3.6 ---------------')
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