现在LU部分解决。

    存在问题是要求H = [A B]之后的B矩阵的顺序主子式不为0，这样才能保证之后B = LU符合要求，否则在S=(U\L\(A·S'))会出现NaN。

    此外，又发现一个问题，有时会存在矩阵奇异，导致不收敛，结果可能出错......不太清楚咋解决。

　

　2012.11.19：仅仅是简单的让H矩阵后半部分B的对角线上元素不为0，可以不出现NaN情况，但有的时候会出现由于收敛问题导致的可能出错情况，有些时候不会出现。