algoritmo de Leverrier
[num,den]=nlev(A,z [,rmax])
matriz quadrada de reais
string
parâmetro opcional (ver bdiag
)
[num,den]=nlev(A,z [,rmax])
computa
(z*eye()-A)^(-1)
por diagonalização por blocos de A seguido pelo algoritmo de Leverrier em cada bloco.
Este algoritmo é melhor que o algoritmo usual de Leverrier, mas ainda não está perfeito!
F. Delebecque., S. Steer INRIA;