lyap equação de Lyapunov

## Seqüência de Chamamento

    [X]=lyap(A,C,'c')
                [X]=lyap(A,C,'d')

## Parâmetros

A, C

  matrizes quadradas de reais, `C` deve ser simétrica

## Descrição

`X= lyap(A,C,flag)` resolve as equações matriciais de tempo contínuo ou de tempo discreto de Lyapunov:

    A'*X + X*A = C       ( flag='c' )
    A'*X*A - X = C       ( flag='d' )

Perceba que existe uma única solução se e só se um autovalor de `A` não é um autovalor de `-A` (`flag='c'`) ou 1 sobre um autovalor de `A` (`flag='d'`).

## Exemplos

    A=rand(4,4);C=rand(A);C=C+C';
    X=lyap(A,C,'c');
    A'*X + X*A -C
    X=lyap(A,C,'d');
    A'*X*A - X -C

## Ver Também

- sylv
- ctr_gram
- obs_gram