# pencan

forma canônica de feixe de matrizes

## Seqüência de Chamamento

    [Q,M,i1]=pencan(Fs)
                [Q,M,i1]=pencan(E,A)

## Parâmetros

Fs

  um feixe regular `s*E-A`

E,A

  duas matrizes quadradas de reais

Q,M

  duas matrizes não-singulares de reais

i1

  inteiro

## Descrição

Dado o feixe regular `Fs=s*E-A`, `pencan` retorna as matrizes `Q` e `M` tais que `M*(s*E-A)*Q` está na forma "canônica".

`M*E*Q` é uma matriz de blocos

    [I,0;
     0,N]

com `N` nilpotente e `i1` = tamanho da matriz acima `I`.

`M*A*Q` é uma matriz de blocos:

    [Ar,0;
     0,I]

## Exemplos

    F=randpencil([],[1,2],[1,2,3],[]);
    F=rand(6,6)*F*rand(6,6);
    [Q,M,i1]=pencan(F);
    W=clean(M*F*Q)
    roots(det(W(1:i1,1:i1)))
    det(W($-2:$,$-2:$))

## Ver Também

- glever
- penlaur
- rowshuff