companion コンパニオン行列

## 呼び出し手順

    A=companion(p)

## 引数

p

  多項式または多項式のベクトル

A

  正方行列

## 説明

特性多項式を有する 行列`A`を返します. `p`がモニックな場合,特性多項式は `p`に等しくなります. `p`がモニックでない場合, `A`の特性方程式は `p/c`に等しくなります. ただし,`c`は`p` の最大次数の係数です.

`p` がモニックな多項式のベクトルの場合, `A` はブロック対角となり, i番目の特性多項式は `p(i)`となります.

## 例

    s=poly(0,'s');
    p=poly([1,2,3,4,1],'s','c')
    det(s*eye(4,4)-companion(p))
    roots(p)
    spec(companion(p))

## 参照

- spec
- poly
- randpencil