﻿功能:采用改进的平方根算法求解AX=B这类线性方程组

格式:SolveEquationLDL(a,b)
a:为对称正定的矩阵 
b:为n*m的矩阵变量,其中n表示a矩阵的阶数,m表示列数,其中m不限制为1,即本函数求解的方程不限定B只为1列的情况

参考://王能超.数值分析简明教程[M].高等教育出版社,北京,2009:186-187

例子:

a =
 [45.6   32.56  21
  32.56  40     17.8
  21     17.8   52.55]
b =
 [2062.7408
  1968.1584
  2498.552 ]

x=SolveEquationLDL(a,b)//回车后得到如下结果
x =
[ 12.3400000000000
  23.7800000000000
  34.5600000000000 ]

//上面的解可以使用{Mul<矩阵运算\Mul>}(a,x)进行验证,看是否等于b

//修改于2012/3/22