使用Gauss高斯消元法求解n元一次方程組的根,

3. 使用Gauss消元法求解n元一次方程組的根,
舉例，三元一次方程組：
0.729x1+0.81x2+0.9x3=0.6867
x1+x2+x3=0.8338
1.331x1+1.21x2+1.1x3=1 

package chapter4;

import java.util.Scanner;

public class demo3 {
    public static void main(String[] args) {
        Scanner sc=new Scanner(System.in);
        System.out.println("輸入未知數的個數");
        int a=sc.nextInt();
        double [][]jz = new double[a][a+1];
        for(int m=0;m<a;m++) {
            for(int n=0;n<a+1;n++) {
                if(n==a) {
                    System.out.println("請輸入常係數");
                    jz[m][a]=sc.nextDouble();
                    continue;
                }
                System.out.printf("請輸入%d行第%d個元的係數\n",m+1,n+1);
                jz[m][n]=sc.nextDouble();
            }
        }
        for(int i=0;i<a;i++) {
            for(int m=i;m<a;m++) {
                double jzmi=jz[m][i];
                for(int n=i;n<a+1;n++) {
                    jz[m][n]=jz[m][n]/jzmi;
                }
            }
            for(int m=i+1;m<a;m++) {
                for(int n=i;n<a+1;n++) {
                    jz[m][n]=jz[m][n]-jz[i][n];
                }
            }
        }
        for(int i=a-1;i>0;i--) {
            for(int m=i-1;m>=0;m--) {
                jz[m][a]=jz[m][a]-jz[i][a]*jz[m][i];
                jz[m][i]=0;
            }
        }
        System.out.print("解向量爲：{");
        for(int m=0;m<a;m++) {
            System.out.print("  "+jz[m][a]);
        }
        System.out.print("}");
    }
}

轉載於：https://www.cnblogs.com/LPworld/p/10723978.html