Cette résolution suive la même démarche que celle de

Gauss a la seule différence que la matrice principale est diagonalisé au lieu d être triangularaisé.