[ W.A.Mozart @ 31.01.2009. 13:54 ] @
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 ]]

]
]
[ holononi @ 04.06.2009. 09:37 ] @
Prvo: Polazna matrica nije zadata u skladu kako se koristi u programu. Potrebno je

A = {
{3, 4, 5, 6},
{9, 7, 8, 6},
{8, 7, 6, 5},
{1, 2, 3, 4}
}

Drugo: Dimenzije matrice za funkciju IdentityMatrix[n] treba odrediti na sledeći način

n = Length[A]

Treće: Potupak izračunavanja treba korigovati.