:ԽԷ

ʽ:SolveSparseDiagonalLU(A,B)

˵:A*x=B,AΪm*nľ,BΪm*1ľ,nʾԽ,ر,n3ʱʱԽǵĴ洢ʽ,nȡ,5Խǡ7Խǡ

ԭ:ȶԷ̽LUֽ,Ȼ

:

//֪
A =
[ 5    2    3    0    0    0    0
  4    8    2    4    0    0    0
  3    2    9    2    1    0    0
  0    2    3    6    2    1    0
  0    0    2    3    7    4    1
  0    0    0    3    1    5    2
  0    0    0    0    2    3    9 ]
y =
[ 7
  6
  4
  2
  0
  8
  7 ]

 A*x=y x

//ɵ


B = [ 0    0    5    2    3
        0    4    8    2    4
        3    2    9    2    1
        2    3    6    2    1
        2    3    7    4    1
        3    1    5    2    0
        2    3    9    0    0 ];
Y = [ 7
        6
        4
        2
        0
        8
        7 ];
x=SolveSparseDiagonalLU(b,y)//سõ½
x =
[ 1.52759453283318
 -0.29997197430259
 -0.01267623852024
  0.57868753503212
 -1.02612857327642
  1.21810460052602
  0.59977148277497 ]