Méthode numérique pour résoudre $$A u = b$$
où $A \in S_n^{++}(\mathbb{R})$ et on décompose $A$ sous la forme de $D + T + {}^t T$ où $D$ est la diagonale de $A$ et $T$ est la partie triangulaire supérieure stricte. On introduit alors $w >0$, $M = D/w + {}^t T$ et $N = (1-w)/wD - T$ de sorte que $A = M-N$.