C++实现行最简型矩阵

昨天学习了线性代数的行最简矩阵的计算方法，思路比较简单，但是做起来比较的麻烦，因此，我查阅了一些资料，编写了一个用C++实现的代码；代码本人已经写好了详细的注释，并且在VS2015编译运行成功，大家可以放心阅读；

// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。

#include "stdafx.h"
#include<math.h>
#include<stdlib.h>
//定义两个宏，M代表矩阵的行数，N代表矩阵的列数
#define M 4
#define N 4   
using namespace std;
//flag数组起到标识作用，用来标识某列的数据是否为0，若flag[n]=0,则该列不为0
int flag[M];
int count1[M];//用来存储3行从开头起所具有的数据为0 的元素的个数
//输出矩阵的函数
void print(float(*a)[N])//以下参数都是用float类型，适应范围更广
{
   
	int i, j;
	for (i = 0; i < M; i++) //二层循环打印
	{
   
		for (j = 0; j < N; j++)
		{
   
			if (fabs(a[i][j]) < 0.000001) //先判断是否为零
				a[i][j] = 0;
			cout << a[i][j] << "   ";
		}
		cout << endl;  //一整行打印完，回车，打印下一行
	}
}
//统计矩阵每行从开头起所具有的连续为0的元素的个数并赋给数组count1，例如0 1 2 0 4，则统计出来为1；
 void statistics(float(*a)[N])
{
   
	int i = 0, j = 0;
	for (i = 0; i<M; i++)