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("}"); } }