# lusolve

solucionador de sistemas lineares esparsos

## Seqüência de Chamamento

    x=lusolve(hand,b)
                x=lusolve(A,b)

## Parâmetros

b

  matriz de reais completa

A

  matriz quadrada de reais esparsa e invertível

hand

  manipulador para fatores de lu esparsos previamente computados (saída de lufact)

x

  matriz de reais completa

## Descrição

`x=lusolve(hand,b)` resolve o sistema linear esparso `A*x = b`.

`[hand,rk]=lufact(A)` é a saída de lufact.

`x=lusolve(A,b)` resolve o sistema linear esparso `A*x = b`

## Exemplos

    non_zeros=[1,2,3,4];rows_cols=[1,1;2,2;3,3;4,4];
    sp=sparse(rows_cols,non_zeros);
    [h,rk]=lufact(sp);x=lusolve(h,[1;1;1;1]);ludel(h)
    rk,sp*x

    non_zeros=[1,2,3,4];rows_cols=[1,1;2,2;3,3;4,4];
    sp=sparse(rows_cols,non_zeros);
    x=lusolve(sp,-ones(4,1));
    sp*x

## Ver Também

- sparse
- lufact
- slash
- backslash