import numpy as np
import matplotlib.pyplot as plt


print('----- 3.3 ----------------')
print("Testons quelques opérations matricielles.")
A = np.array([[2, 3], [4, 5]], np.float64)
print("A=\n", A)
B = np.array([[3, 5], [7, 9]], np.float64)
print("B=\n", B)
C = A * B # produit terme à terme
print("C=\n", C)
D = np.dot(A, B) # produit matriciel
print("D=\n", D)
print("type(D[0,0])=", type(D[0,0]))
vec = np.array([4, 6], np.float64)
print("A*vec=", np.dot(A, vec)) # produit Matrice fois un vecteur.

print('----- 3.4 ----------------')
Ainv = np.linalg.inv(A) # Matrice inverse de la matrice A
print("Ainv=\n", Ainv)
print("A*Ainv=\n", np.dot(A, Ainv)) # Matrice fois son inverse donne l'identité