glever inverso do feixe de matrizes

## Seqüência de Chamamento

    [Bfs,Bis,chis]=glever(E,A [,s])

## Parâmetros

E, A

  duas matrizes de reais quadradas de igual dimensão

s

  string (o valor padrão é '`s`')

Bfs,Bis

  duas matrizes de polinômios

chis

  polinômio

## Descrição

Computação de

`(s*E-A)^-1`

pelo algoritmo generalizado de Leverrier para um feixe de matrizes.

    (s*E-A)^-1 = (Bfs/chis) - Bis.

`chis` = polinômio característico (até uma constante multiplicativa).

`Bfs` = nmatriz de polinômios como numerador.

`Bis` = matriz de polinômios ( - expansão de `(s*E-A)^-1` ao infinito).

Note o sinal - antes de `Bis`.

## Cuidado

Esta função usa `cleanp` para simplificar `Bfs,Bis` e `chis`.

## Exemplos

    s=%s;F=[-1,s,0,0;0,-1,0,0;0,0,s-2,0;0,0,0,s-1];
    [Bfs,Bis,chis]=glever(F)
    inv(F)-((Bfs/chis) - Bis)

## Ver Também

- rowshuff
- det
- invr
- coffg
- pencan
- penlaur