Algorithme de Leverrier
[num,den]=nlev(A,z [,rmax])
matrice réelle carrée
chaîne de caractères
paramètre optionnel (voir bdiag)
[num,den]=nlev(A,z [,rmax]) calcule
(z*eye()-A)^(-1) par bloc diagonalisation de
A suivie de l'algorithme de Leverrier sur chaque
bloc.
Cet algorithme est plus efficace que l'algorithme de Leverrier seul, mais il n'est toujours pas parfait !