[求助]LU分解怎么得出的结果和计算的不一致啊？

hoarywolf

使用lu(A)函数作LU分解，可是分解出来的L,U矩阵却是这样的，请问怎么回事啊？<BR>A =
<p>
<>2 4 2 6<BR>4 9 6 15<BR>2 6 9 18<BR>6 15 18 40</P>
<>L =</P>
<>0.3333 1.0000 -0.6667 1.0000<BR>0.6667 1.0000 0 0<BR>0.3333 -1.0000 1.0000 0<BR>1.0000 0 0 0<BR>U =</P>
<P>6.0000 15.0000 18.0000 40.0000<BR>0 -1.0000 -6.0000 -11.6667<BR>0 0 -3.0000 -7.0000<BR>0 0 0 -0.3333<BR>可是我用杜立特算法手工算出的却是<BR>L=<BR>1<BR>2 1<BR>1 2 1<BR>3 3 2 1</P>
<P>U=<BR>2 4 2 6<BR>1 2 3<BR>3 6<BR>1</P>
<P>怎么会出现这种结果呢？难道是lu(A)函数的求LU的算法和我的算法不一致吗？<BR></P>