Не работает функция решения СЛАУ методом Гаусса, переведенная из рабочего кода C++

// √аусса.cpp: определ¤ет точку входа дл¤ консольного приложени¤.
//
 
#include "stdafx.h"
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
 
#define N    5    
#define N1    N+1  
 
float matrix[N][N1]={{ 5,    1.2,    -2,   2.4,   1,   9  },
                     { 1.8,    5,     1,   2.2,   4,   7  },
                     { 3,    1.8,     1,   1.6,   1,   2.5},
                     { 4.5,    2,     4,   2.6,   2,   9  },
                     { 7,      1,   2.4,     0,   9,  -8  }}; 
 
 
void ShowMatrix(void)
{ 
    printf("systema:    \n");
    int i, j;
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
            printf("%+3.1f*x%d ", matrix[i][j],j + 1);
        printf(" = %3.1f\n", matrix[i][N]);
    }
}
 
 
int _tmain(int argc, _TCHAR* argv[])
{
      float tmp,xx[N1];
      short int i,j,k;
 
      ShowMatrix();
  
      /*ћетод √аусса*/
  
     for (i = 0; i < N; i++)
     {
        tmp = matrix[i][i];
        for (j = N; j >= i; j--)
            matrix[i][j] /= tmp;
 
        for (j = i + 1; j < N; j++)
        {
          tmp = matrix[j][i];
          for (k = N; k >= i; k--)
            matrix[j][k] -= tmp * matrix[i][k];
        }
     }
 
     xx[N-1] = matrix[N-1][N];
     for (i = N-2; i >= 0; i--)
     {
         xx[i] = matrix[i][N];
         for (j = i + 1; j < N; j++) 
             xx[i] -= matrix[i][j] * xx[j];
     }
 
     printf("\n       Metod Gaussa:    \n");
      for (i = 0; i < N; i++)
        printf("x%d = %3.3f\n", i + 1, xx[i]);
     return 0;  
}