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