# lusolve

疎な線形システムの解を得る

## 呼び出し手順

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

## パラメータ

b

  通常の実数行列

A

  可逆な実数正方疎行列

hand

  疎行列LU分解を計算した際のハンドル(lufactの出力)

x

  通常の実数行列

## 説明

`x=lusolve(hand,b)` は疎な線形システム `A*x = b`を解きます.

`[hand,rk]=lufact(A)` は lufact の出力です.

`x=lusolve(A,b)`は,疎な線形システム `A*x = b`を解きます.

## 例

    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

## 参照

- sparse
- lufact
- slash
- backslash