sylv シルベスタ方程式.

## 呼び出し手順

    sylv(A, B, C, flag)

## 引数

A,B,C

  適当な次元の実数行列.

flag

  文字列 (`'c'` または `'d'`)

## 説明

`X= sylv(A,B,C,'c')` は以下の "連続時間" シルベスタ方程式の解`X`を計算します,

    A*X+X*B=C

`X=sylv(A,B,C,'d')` は以下の "離散時間" シルベスタ方程式の解`X`を計算します,

    A*X*B-X=C

## 例

    A = rand(4, 4); C = rand(4, 3); B = rand(3, 3);
    X = sylv(A, B, C, 'c');
    norm(A*X+X*B-C)
    X = sylv(A, B, C, 'd');
    norm(A*X*B-X-C)

## 参考

- lyap