Metod Danilevskog i Mathematica 6
Trebam uraditi model racunanja koeficijenata karakteristicnog polinoma matrice metodom Danilevskog. Model treba biti uradjen u okruzenju Mathematica. Uspio sam sklepati neki kod medjutim negdje mi se potkrada neki propust. Ima li neka dobra dusa da me prosvijetli ili mozda da drugi prijedlog kako na efikasniji nacin ovo uraditi. Moram napomenuti da nisam zavrsio finalno "vadjenje" koeficijenata iz zadnje iteracije.Unaprijed zahvalan svima.
Evo kod :
A[0] = {
{3, 4, 5, 6},
{9, 7, 8, 6},
{8, 7, 6, 5},
{1, 2, 3, 4}
}
n = Dimensions[A][[1]]
For[ i = 1, i <= n, i++,
For [j = 1, j <= n, j++,
F [ i ] = IdentityMatrix[n]
F[ i ] [[n - i, j]] = -(A[i - 1][[n - i + 1, j]]/A[i - 1][[n - i + 1, n - i]])
F[ i ] [[n - i, n - i]] = 1/A[i - 1][[n - i + 1, n - i]]
B[ i ] = Inverse [F[ i ]]
A[ i ] = B[ i ].A[i - 1].F[ i ]
Print[A[ i ]]
]
]
 |  |
 |
