#importation du module numpy
import numpy as np

#écriture des coefficients de la matrice M dans un tableau M, ici pour un système 5x5,
#d'autres dimensions possibles
M=np.array([[1,-1,0,0,-1],[0,1,1,-1,0],[0,0,-20,-30,0],[-10,-50,20,0,0],[0,50,0,30,-40]])

#écriture du vecteur des coefficients constants b dans un tableau b, système 5x5
b=np.array([0,0,-10,0,0])

#appel de la méthode permettant de résoudre le système (M,b)
sol=np.linalg.solve(M,b)

#affichage des résultats, noms des réponses selon donnée problème
for i in range(0,len(b)):    
    print('I'+str(i+1)+'=',sol[i])