Ранг матрицы

#include <algorithm>
#include <iostream>
#include<stdio.h>
#include<math.h>
#include<conio.h>
#include<stdlib.h>
#include <time.h>
#define NN 6
float rndup(float n)//round up a float type and show one decimal place
{
      float t;
      t=n-floor(n);
      if (t>=0.5)    
      {
              n*=10;//where n is the multi-decimal float
              ceil(n);
              n/=10;
              }
      else 
      {
              n*=10;//where n is the multi-decimal float
              floor(n);
              n/=10;
              }
      return n;
}     
int main()
{
 
    srand(time(NULL));
 float a[NN][NN],kst;
 bool key = true;
int i,j,h,n,p,h1,k=0,t;
std::cout<<"Vvedite razmernost"<<endl;
std::cin>>n;
    for(i=0;i<n;i++)
for(j=0;j<n;j++)
//cin>>a[i][j];
a[i][j]=1+rand()%9;
 
 
//--------Вывод сгенерированной матрицы-----//
for(i=0;i<n;i++)
{
for(j=0;j<n;j++) 
std::cout<<a[i][j]<<" ";
std::cout<<endl;
}
 
//---<Подсчет ранга>----/
for (i = 0; i < n; ++i)
for (k=i+1;k<n; k++)
{
  double multi = a[k][i] / a[i][i];
for(j=0; j<n; j++)
{
     double i2j = a[k][j];
                     i2j = rndup(i2j);
 
                     double   mulxmj = (multi * a[i][j]);
                    mulxmj = rndup(mulxmj);
                       a[k][j] = i2j - mulxmj;
}
}
std::cout<<endl;
std::cout<<endl;
//-----Вывод полученной матрицы---//
for(i=0;i<n;i++)
{
    for(j=0;j<n;j++)
        printf("%4.0f\t",a[i][j]);
printf("\n\n");
    }
std::cout<<endl;
 
//---подсчет ранга--- //
int rang = 0;
key = true;
 
for ( i=0; i<n; ++i)
{
        key = false;
    for ( j=0; j<n; ++j)
                if (a[i][j] != 0.0)
                        key = true;
                if (!key)
                        h1++;
}
rang = n - h1;
 
std::cout <<"rang = "<< rang << "\n";
std::cout <<"\n";
 
return 0;
}