![$\left[ \begin{array}{ccc} 1.65 & 1.9 & -15 \\ 12 & 0 & 0.1 \\ 2.85 & -23 & -3.15 \end{array} \right] \mathbf{x} = \left[ \begin{array}{c} 90.9\\ -3.7 \\ -42.6 \end{array} \right]$ $\left[ \begin{array}{ccc} 1.65 & 1.9 & -15 \\ 12 & 0 & 0.1 \\ 2.85 & -23 & -3.15 \end{array} \right] \mathbf{x} = \left[ \begin{array}{c} 90.9\\ -3.7 \\ -42.6 \end{array} \right]$](https://dxdy-01.korotkov.co.uk/f/c/d/0/cd0ed6e38947355c401ce8128aa3884382.png)
Что бы найти обратную диагональную матрицу поменяем местами 1-й и 2-й столбцы (в векторе-столбце b тоже не забуду потом поменять).
![$\mathbf{D}^{-1} = \left[ \begin{array}{ccc} 12 & 0 & 0\\ 0 & 1.9 & 0 \\ 0 & 0 & -3.15 \end{array} \right]^{-1}=\left[ \begin{array}{ccc} 0.083 & 0 & 0\\ 0 & 0.5263 & 0 \\ 0 & 0 & -3.0175 \end{array} \right]$ $\mathbf{D}^{-1} = \left[ \begin{array}{ccc} 12 & 0 & 0\\ 0 & 1.9 & 0 \\ 0 & 0 & -3.15 \end{array} \right]^{-1}=\left[ \begin{array}{ccc} 0.083 & 0 & 0\\ 0 & 0.5263 & 0 \\ 0 & 0 & -3.0175 \end{array} \right]$](https://dxdy-02.korotkov.co.uk/f/5/4/3/5436cc489ae8ab57c81ebfaa588ea8e682.png)
![$\mathbf{x}^{k+1}=-\mathbf{D}^{-1}(\mathbf{L}+\mathbf{U})\mathbf{x}^{k}+\mathbf{D}^{-1}\mathbf{b}=-\left[ \begin{array}{ccc} 0.083 & 0 & 0\\ 0 & 0.5263 & 0 \\ 0 & 0 & -3.0175 \end{array} \right]\left[ \begin{array}{ccc} 0 & 0 & 0.1\\ 1.65 & 0 & -15 \\ 2.85 & -23 & 0 \end{array} \right]\mathbf{x}^{k}+\left[ \begin{array}{ccc} 0.083 & 0 & 0\\ 0 & 0.5263 & 0 \\ 0 & 0 & -3.0175 \end{array} \right]\left[ \begin{array}{c} -3.7\\ 90.9 \\ -42.6 \end{array} \right]=\left[ \begin{array}{ccc} 0 & 0 & -0.0083\\ -0.8684 & 0 & 7.8945 \\ 8.5999 & -69.4025 & 0 \end{array} \right]\mathbf{x}^{k}+\left[ \begin{array}{c} -0.3071\\ 47.8407 \\ 128.5455 \end{array} \right]$ $\mathbf{x}^{k+1}=-\mathbf{D}^{-1}(\mathbf{L}+\mathbf{U})\mathbf{x}^{k}+\mathbf{D}^{-1}\mathbf{b}=-\left[ \begin{array}{ccc} 0.083 & 0 & 0\\ 0 & 0.5263 & 0 \\ 0 & 0 & -3.0175 \end{array} \right]\left[ \begin{array}{ccc} 0 & 0 & 0.1\\ 1.65 & 0 & -15 \\ 2.85 & -23 & 0 \end{array} \right]\mathbf{x}^{k}+\left[ \begin{array}{ccc} 0.083 & 0 & 0\\ 0 & 0.5263 & 0 \\ 0 & 0 & -3.0175 \end{array} \right]\left[ \begin{array}{c} -3.7\\ 90.9 \\ -42.6 \end{array} \right]=\left[ \begin{array}{ccc} 0 & 0 & -0.0083\\ -0.8684 & 0 & 7.8945 \\ 8.5999 & -69.4025 & 0 \end{array} \right]\mathbf{x}^{k}+\left[ \begin{array}{c} -0.3071\\ 47.8407 \\ 128.5455 \end{array} \right]$](https://dxdy-04.korotkov.co.uk/f/7/8/e/78e32bf16c2b8997dc4f4b202395befa82.png)
Посмотрите пожалуйста, правильно ли я делаю.
Дело в том, что написал алгоритм, взял другой пример который решен и прогнал - все работает. Но для вот этого примера ничего не работает -элементы столбца x стремятся к бесконечности, а не к корням! Помогите, я уже совсем измаялся.