数值分析课题二 矩阵求逆

模拟求逆过程写就可以了

e是单位矩阵，只通过行变换，把原矩阵变成单位矩阵。

消元的时候不能选主元。

#include <iostream>
#include <cstdio>
#include <math.h>
using namespace std;
const int MAXN = 1000;
int n;
double a[MAXN][MAXN];
double e[MAXN][MAXN];
void init()
{
    for(int i = 0;i<n;i++){
        for(int j = 0;j<n;j++){
            e[i][j] = i == j?1:0;
        }
    }


}
void low_mat()
{
    for(int k = 0;k<n;k++){///往下扫
        for(int i = k+1;i<n;i++){
            double l = -a[i][k]/a[k][k];
            for(int j = k;j<n;j++){
                a[i][j] += l*a[k][j];
                e[i][j] += l*e[k][j];
            }
        }
    }
}
void up_mat()
{
    for(int k = n-1;k>=0;k--){///再往上扫
            for(int i = k-1;i>=0;i--){///千万别选主元，会跪
                doubl