Factorisation de matrices LDM LDL(transpose) - Page 1

Bonjour voilà je cherche à réaliser en TI-basic ces deux algorithmes (matlab/octave ,mais je pense que ça se comprend) à moins que quelqu'un ne connaisse déjà des programmes tout fait ça m'arrangerait (mais même sur ticalc : nada ). N'ayant pas vraiment le temps de devenir un expert en TI-basic, je reste un malheureux novice, si des personnes plus expérimentés pouvaient me venir en aide je leur en serai très reconnaissantes.

LDM

d1 = A11
L2:n,1 = A2:n,1/A11
for j = 2 : n − 1 do
solve L1:j,1:j v1:j = A1:j,j (Systeme triangulaire)
Mj,1:j−1 = v1:j−1./d1:j−1
dj = vj
Lj+1:n,j = (Aj+1:n,j − Lj+1:n,1:j−1 v1:j−1)/vj
end for
solve L1:n,1:n v1:n = A1:n,n (Systeme triangulaire)
Mn,1:n−1 = v1:n−1./d1:n−1
dn = vn



LDLt

for j = 1 : n do
for i = 1 : j − 1 do
vi = Lj,i di
end for
vj = Ajj − Lj,1:j−1 v1:j−1
dj = vj
Lj+1:n,j = (Aj+1:n,j − Lj+1:n,1:j−1 v1:j−1)/vj
end for


Merci d'avance

Il manque les parenthèses autour des indices dans ton code Matlab, c'est limite illisible.

Erf désolé j'ai repris l'algo de mon prof qui fait également des cours illisibles ... Plus sérieusement voici deux screenshots à partir de mon LaTeX

Fichier joint : (LDL.PNG)

Fichier joint : (LDM.PNG)