Метод на Якоби (проста итерация) за решаване на СЛАУ

Clear["Global*`"]
A=(({{4, 1, -2}, {1, 9, -3}, {2, 3, -10}}));
b={1,0,-2};
n=Length[b]

съставяне на матрицата C и вектора d от итерационна формула x^(k+1)=C x^k +d

C[[i,i]]=0
C[[i,j]]=-((A[[i,j]])/(A[[i,i]]))
C[[i]]=-((A[[i]])/(A[[i,i]]))
d[[i]]=((b[[i]])/(A[[i,i]]))

CC={{0,0,0},{0,0,0},{0,0,0}};
d={0,0,0};
For[i=1, i≤n, i++,
(*body*)
CC[[i]]=-((A[[i]])/(A[[i,i]]));
CC[[i,i]]=0;
d[[i]]=((b[[i]])/(A[[i,i]]))

]
MatrixForm[CC]
d

проверка на първа норма:

Max[(∑_(j=1)^n)Abs[CC[[1,j]]],(∑_(j=1)^n)Abs[CC[[2,j]]], (∑_(j=1)^n)Abs[CC[[3,j]]]]

Извод: първа норма на CC е по-малка от 1 => методът ще е сходящ

Итерации:

x={0,0,0};
Print["k = ", 0," x = ", x]
eps =(10^(-6));
epsk =3;
For[k=1, epsk >eps, k++,
(*body*)
xnew = CC.x +d;
epsk=Max[Abs[xnew-x]];
Print["k = ", k," x = ", N[xnew], " epsk = ", N[epsk]];
x = xnew
]

Получихме решението: x1=0.382759, x2 = 0.0551726, x3=0.293103 с точност до 6-ия знак след десетичната точка