脉络最清晰的C++高斯消元程序

#include<iostream>
#include<math.h>
using namespace std;
const int N=3;
void print(double a[][N],int n);//输出程序
int max_absAKK(double b[][N],int n,int kk);//寻找列中的最大元素所在的行   函数
void ExchangCol(double c[][N],int p, int q,int n);//换行   函数
int main()
{
	double A[N][N]={{0.001,2.000,3.000},{-1.0,3.712,4.623},{-2.0,1.072,5.643}};
	double m[N][N];
	for(int k=0;k<N-1;k++)
	{
	int t=max_absAKK(A,N,k);//找到每一列绝对值最大的行
	if (t!=k)   ExchangCol(A,t,k,N);//如果对角线的元素不是最大,进行换行
	for(int i=k+1;i<N;i++)
	{
		m[i][k]=A[i][k]/A[k][k];
		for(int j=k+1;j<N;j++)
		{
			A[i][j]=A[i][j]-m[i][k]*A[k][j];
			
		}
		
	}
}
	print (A,N);
	
}

void ExchangCol(double c[][N],int p, int q,int n)
{
	for(int i=0;i<n;i++)
	swap(c[p][i],c[q][i]);
}
int max_absAKK(double b[][N],int n,int kk)
{
	double max=0.0;
	int maxCol=0;
	for(int i=0;i<n-kk;i++)
{
if(fabs(max)<fabs(b[kk+i][kk]))
{
	max=b[kk+i][kk];
	maxCol=kk+i;	
	}
	
}
return maxCol;
}

void print(double a[][N],int n)
{
	for (int i=0;i<n;i++)
	{
	for(int j=0;j<n;j++)
	{
		if(fabs(a[i][j])<1.0e-5)
		a[i][j]=0.0;
	cout<<a[i][j]<<'\t';

}
cout<<endl;
}	
}