杜利特尔分解(LU分解）的JAVA代码

杜利特尔分解(LU分解）的JAVA代码

数学基础：设A为n阶方阵，满足其各阶主子式：
det (Subscript[A, p]) != 0
A=LU

其代码如下：
public static double[][] LUDecomposition(double a[][],String X)//X代表返回矩阵的种类，若为U，返回u矩阵；若为l,返回l矩阵。若为其他，返回a;

    {
    int row=a.length;//row为矩阵行数
    int line=a[0].length;//line为矩阵列数
    double temp=0;
    if(row!=line) {
        System.out.println("矩阵行列数不等！！");
        return a;
        }//如果矩阵行列数不等，返回原矩阵。
    double[][] l=new double[row][line];//定义L矩阵
    double[][] u=new double[row][line];//定义U矩阵
    int k=0;
    for(int j=0;j

下面在主函数里调用这个函数：
public static void main(String[] args) {
double[][] a= {{6, 10, 7}, {1, 7, 8}, {10, 6, 3}};double[] b= {1,0,-3};
double[][] l=LUDecomposition(a,”l”);
double[][] u=LUDecomposition(a,”u”);
double[] y=SolveLy(l,b);
double[] x=SolveUx(u,y);
}