Найти ошибку в коде Wolfram Mathematica

[Yuriyy]

Необходимо проверить метод простых итераций. Но код не работает. Объясните, пожалуйста, что с ним не так

SI[Matr_, Vect_, eps_] :=
Module[{q = Norm[Matr, 2], f = Vect, h = Matr, xs, xn = {}, k = 0},
xs = f; xn = h.xs; q = q/(1 - q); k++
While[(q*Max[Abs[xs - xn]] >= eps), xs = xn; xn = h.xs + f; k++];
Print[k "- количество шагов"]; Return[xn]]
S1 = N[SI[H, f, eps]]

[ura_111]

Я только что скачал эту математику...

Так?

45.jpg

p.s.: я не знаю эту программу - у меня методом тыка что-то вывело. Скажи что должно получится.

[ura_111]

Начал осваиваться с этой программой. Поставь запятую:

46.jpg

ЩАС РАЗБЕРЁМСЯ.

[ura_111]

Скажи мне чему равны эти параметры [H, f, eps] ?
И что нужно получить? Число?

70.jpg

p.s.: есть подозрения, что параметры [H, f, eps] не определены.

[Yuriyy]

Сообщение от ura_111

Начал осваиваться с этой программой. Поставь запятую:

Вложение 83592

ЩАС РАЗБЕРЁМСЯ.

Насчет запятой спасибо, не знаю, как можно было тут ошибку допустить :D
программа сама работает, но количество итераций не может быть равно 0.

[Yuriyy]

Сообщение от ura_111

Скажи мне чему равны эти параметры [H, f, eps] ?
И что нужно получить? Число?

Вложение 83593

p.s.: есть подозрения, что параметры [H, f, eps] не определены.

всё определено, H- матрица, f- столбец свободных членов, eps - погрешность, eps = 1/2*10^-9

[ura_111]

А что такое "k++" в строчке:

xs = f; xn = h.xs; q = q/(1 - q); k++

на картинке такого нет.