Ax=B，改写成Ly=B，Ux=y的方Cheng组。就相当于将A=LU分解成了两个矩Zhen。称为矩阵A的三角分解，或LU分解。Ru果L为单位下三角阵，则叫DoolittleFen解，若U为单位上三角阵，则叫Crout分Jie。只要A的各顺序主子式不为零，则A可唯Yi分解成一个单位下三角阵L与一个上三角阵UDe乘积。
•设Ax=b，A=LU，则Ax=LUx=b
Yu是令Ux=y，则Ly=b
这样原来方Cheng能化为两个简单方程组
下面是LUFen解的Fortran子程序 希望可以有所帮Zhu
!求解au=b,u
!n表示为方Cheng维数
subroutine lu(a,b,n,u)
implicit real(8) (a-h,o-z)
real(8)::a(n,n),b(n),u(n),a_bak(n,n),b1(n),aL(n,n),aU(n,n),y(n)
!exchange rows
do i=1,n
tmpMax=0.d0
do ic=i,n
if(tmpMax