Неправильно считается ранг матрицы

#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"<<std::endl;
std::cin>>n;
for(i=0;i<n;i++)
{for(j=0;j<n;j++)
{
std::cout<<"\n["<<i<<"]["<<j<<"]=";
std::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<<std::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<<std::endl;
std::cout<<std::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<<std::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;
rang++;}
else
{rang--;}
 
}
 
 
std::cout <<"rang = "<< rang-1 << "\n";
std::cout <<"\n";
    /*
int rang = 1;
int k1=n-1;
int rang1=1;
//key = true;
/*
for(int i=0;i<k1;i++)
for(int j=0;j<n-1;i++)
if (a[i+1][j+1]!=0)
{
rang++;
if (rang>rang1)
rang1=rang;}
else
{
//k1--;
i++;
if (rang>rang1)
rang1=rang;}
 
std::cout <<"rang = "<< rang1 << "\n";
std::cout <<"\n";  */
getchar();
getchar();
}

#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"<<std::endl;
std::cin>>n;
for(i=0;i<n;i++)
{for(j=0;j<n;j++)
{
std::cout<<"\n["<<i<<"]["<<j<<"]=";
std::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<<std::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<<std::endl;
std::cout<<std::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<<std::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 ) 
        rang++;
}
 
 
std::cout <<"rang = "<< rang << "\n";
std::cout <<"\n";
    /*
int rang = 1;
int k1=n-1;
int rang1=1;
//key = true;
/*
for(int i=0;i<k1;i++)
for(int j=0;j<n-1;i++)
if (a[i+1][j+1]!=0)
{
rang++;
if (rang>rang1)
rang1=rang;}
else
{
//k1--;
i++;
if (rang>rang1)
rang1=rang;}
 
std::cout <<"rang = "<< rang1 << "\n";
std::cout <<"\n";  */
getchar();
getchar();
}