Scilab Function

lusolve - solveur de système linéaire creux

### Sequence d'appel

  `lusolve(hand,b)`
  `lusolve(A,b)`

### Parametres

- `b` : matrice réelle pleine
- `A` : matrice réelle carrée creuse inversible
- `hand` : pointeur vers des facteurs L,U creux déjà calculés (sortie de lufact)

### Description

résout le système linéaire `A*x = b`, où `hand` est obtenu par un appel préliminaire à `lufact` :

`[hand,rk]=lufact(A)`.

### Exemples

    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

### Voir aussi

`sparse` `lufact` `slash` `backslash`