如何用乔列斯基分解求下三角矩阵？

楚香饭

程序说明：
CHODCM(A,N,D,T)
N     整型变量，输入参数，矩阵A的阶数；
A      NxN个元素的二维实型数组，输入、输出参数，输入时存放A的上三角部分，严格下三角部分可置为任意数，输出时，严格下三角部分，存放分解L的严格下三角 ，而上三角部分不变；
D   N个元素的一维实型数组，输出参数，存放分解中对角均值D的对角元素；
T    N个元数的一维实型数组，工作单元。