试推导矩阵A的Crout分解A=LU的计算公式，其中L为下三角阵，U为单位上三角阵。