数值方法高斯若当消元法解线性方程组
import java.util.Scanner;
public class Gauss {
static final int MAXN= 20;
static double a[][]=new double[MAXN][MAXN];
static int num;
public static void main(String[] args) {
System.out.println("输入未知数个数:");
Scanner sc=new Scanner(System.in);
num=sc.nextInt();
System.out.println("输入用矩阵表示的线性方程组:");
for (int i = 0; i < num; i++) {//输入方程组
for (int j = 0; j <= num; j++) {
a[i][j]=sc.nextDouble();
}
}
for (int j = 0; j < num-1; j++)
{
guiyi(j);
dengjia(j);
guiling(j+1);
}
guiyi(num-1);
for(int i=0;i<num;i++)
{
System.out.println("x"+i+"="+a[i][num]);
}
}
static void dengjia(int jj) {//把原矩阵转化成三角形的
for (int i = jj + 1; i < num; i++) {
double k = a[i][jj]/ a[jj][jj];
for (int j = jj; j =0;i--){
a[jj][i]=a[jj][i]/a[jj][jj];
}
}
static void guiling(int jj){//使其他值变为0
for(int i=0;i<jj;i++){
double k=(-a[i][jj])/a[jj][jj];
for(int j=jj;j<=num;j++){
a[i][j]+=a[jj][j]*k;
}
}
}
}
作者：菜虚困爱菜籽
相关文章
Brigitte
2021-08-06
Serwa
2020-03-20
Opal
2020-03-09
Iolanthe
2020-02-24
Rhoda
2023-07-22
Hester
2023-07-22
Grace
2023-07-22
Vanna
2023-07-22
Peony
2023-07-22
Dorothy
2023-07-22
Dulcea
2023-07-22
Zandra
2023-07-22
Serafina
2023-07-24
Kathy
2023-08-08
Olivia
2023-08-08
Elina
2023-08-08
Jacinthe
2023-08-08
Viridis
2023-08-08
Hana
2023-08-08