主元消去法（N元一次方程组）



代码是我自己写的，发现错误或者bug的请评论告诉我，谢谢

学过线代的话，就会知道算法很简单，把方程组的系数和常数放在一个矩阵里，进行

矩阵变换，最后得到的单位矩阵就是方程组的解

#include<stdio.h>
#include<math.h>

#define N 3

float A[N][N+1];

void findmain (int n)    //寻找主元函数，将方程组系数的最大放在矩阵前面，并前除以自己变成一
{
    float a,b;
 a = fabs (A[n][n]);
    int i;
 int j = n;

 if (n != N - 1)
 {
     for (i = n; i < N; i++)
  {
      if (fabs(A[i][n]) > a)
   {
       a = A[i][n];
                j = i;
   }
  }
 }
 if (j != n)
 {
     for (i = n; i < N + 1; i++)
  {
      b = A[n][i];
   A[n][i] = A[j][i];
   A[j][i] =