Bonjour,

Dans un exercice de maths je dois retranscrire l'algorithme de résolution, sur Python, de l'équation
A*x=B avec A(n,n) et B(n,1) où A est déjà une matrice triangulaire supérieure.
Par exemple :
A = 1 1 1 1    B=4
         0 1 0 1          2
         0 0 1 2          4
         0 0 0 3          2
J'ai donc suivi l'algorithme qui nous est donné pour l'écrire en Python :

def Remontee(A,B):
    n=A.shape[0]    Je récupère la dimension de A
    x=np.zeros([n]) Je créer ma liste de X vide
    x[n-1]=b[n-1]/a[n-1,n-1] J'initialise mon dernier X
    for k in range (2,n):
       for j in range(n-k,n-1):
            b[n-k]=b[n-k]-a[n-k,j]*x[j]
        x[n-k]=b[n-k]/a[n-k,n-k]
    return x

Mais ça ne fonctionne par car en comparant le résultat avec la fonction solve_triangular je trouve un résultat différent.
Je pense que je fais une erreur dans la partie en rouge, mais j'ai beau cherché je ne trouve pas mon erreur.

Pourriez vous m'aider ?