Le problème se trouvait autre part!
En effet le problème venait de la formule 9 page 11 du complément. Dans cette formule on utilise l'inverse de la matrice Ia. Le déterminant de cette dernière est de l'ordre de 10^-7 alors qu'il nous a été demandé de retourner la matrice nulle dans la méthode d'invertion si le déterminant de la matrice était inférieur à 10^-5.
En tout cas merci de la réponse!