nlev Leverrierのアルゴリズム

## 呼び出し手順

    [num,den]=nlev(A,z [,rmax])

## 引数

A

  実数正方行列

z

  文字列

rmax

  オプションのパラメータ (`bdiag`参照)

## 説明

`[num,den]=nlev(A,z [,rmax])` は, `(z*eye()-A)^(-1)`を計算します.

計算の際には,Aのブロック対角化の後に Leverrierアルゴリズムを各ブロックに適用します.

このアルゴリズムは通常の Leverrierアルゴリズム より優れていますが,まだ完全ではありません!

## 例

    A=rand(3,3);x=poly(0,'x');
    [NUM,den]=nlev(A,'x')
    clean(den-poly(A,'x'))
    clean(NUM/den-inv(x*eye()-A))

## 参照

- coff
- coffg
- glever
- ss2tf