Leverrier's algorithm
[num,den]=nlev(A,z [,rmax])
real square matrix
character string
optional parameter (see bdiag
)
[num,den]=nlev(A,z [,rmax])
computes
(z*eye()-A)^(-1)
by block diagonalization of A followed by Leverrier's algorithm on each block.
This algorithm is better than the usual leverrier algorithm but still not perfect!