bin.
Спасибо.
Про перестановочную матрицу я знаю, видел ее во многих статьях, в том числе и в Вашей.
Всем спасибо за плодотворную дискуссию.
Свою программу по перестановке строк и матриц я исправил. Теперь она правильно работает на всех матрицах.
(Оффтоп)
Для интересующихся привожу код новой программы:
Код:
Perm2Col[a_, j1_, j2_] :=
Module[{n, am, oldStr1, oldColn1},
am = a; n = Length[am];
oldStr1 = Table[0, {k, 1, n}]; oldColn1 = Table[0, {k, 1, n}];
Do[oldStr1[[j]] = am[[j1, j]]; am[[j1, j]] = am[[j2, j]];
am[[j2, j]] = oldStr1[[j]], {j, 1, n}];
Do[oldColn1[[i]] = am[[i, j1]]; am[[i,j1]] = am[[i, j2]];
am[[i, j2]] = oldColn1[[i]], {i,1, n}];
am]
Таким образом, можно считать, что тема закрыта.
