【Java例题】4.3 3. 使用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个元的系数
",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