请教高手

从数学上可以证明：对称非奇异矩阵A可以分解为下三角矩阵L与其转置LT之积，即A=LLT。
，做一个函数(函数名为LLT)，参数为对称非奇异矩阵，返回值为下三角矩阵，即，
在matrix.h中申明原型：   
CMatrix LLT(const CMatrix& A);//不是类的成员函数，是一般的函数
在matrix.cpp实现：
CMatrix LLT(const CMatrix& A)
{
  …………………………………
}
要在主函数中测试。
提示：先推导公式，再编写程序