С.П. Шарый - Институт вычислительных технологий СО РАН

3.6. Прямые методы решения линейных систем 295⎛C =⎜⎝⎞↓0↓ ↓ . ↓ ↓ ..,. ... . .. ⎟ ↓ ⎠ ··· ×Рис. 3.15. Схема определения элементов в методе Холесского.Более точно,при j = 1{c11 = √ a 11 ,c i1 = a i1 /c 11 ,i = 2,3,...,n,при j = 2{c22 = √ a 22 −c 2 21 ,c i2 = ( a i2 −c i1 c 21)/c22 ,i = 3,4,...,n,при j = 3{c33 = √ a 33 −c 2 31 −c2 32 ,c i3 = ( a i3 −c i1 c 31 −c i2 c 32)/c33 ,i = 4,5,...,n,и так далее для остальных j. Псевдокод этого процесса выглядит следующимобразом:DO FOR j = 1 TO n∑j−1c jj ← √ ajj −k=1c 2 jkDO FOR i = j +1 TO n( )∑j−1/c ij ← a ij − c ik c jk cjjEND DOEND DOk=1(3.72)