Cixi

βi ci−1 , i = 2, . . . , n. (5.33) α1 = a1 , βi = αi−1 Avec (5.33), il est facile de résoudre les deux systèmes bidiagonaux Ly = b et Ux = y, pour obtenir les formules suivantes y1 = b1 , yi = bi − βi yi−1 , i = 2, . . . , n, (5.34) yn , xi = (yi − cixi+1 )/αi, i = n − 1, . . . , 1. (5.35) (Ux = y) xn = αn (Ly = b) 154 5 Systèmes linéaires Cette technique est connue sous le nom d’algorithme de Thomas. Son coût est de l’ordre de n opérations. La commande spdiags de MATLAB permet de construire