invalid floating point operation

#include <iostream.h>
#include <math.h>
#include <stdlib.h>
 
int main(int argc, char* argv[])
{
   int N;
   int i, j;
   int kol = 0;
   double sum = 0;
 
   // задание в качестве первого числа последовательности текущего времени
   srand(time(NULL));
 
   cout << "Enter dimension of a array" << endl;
   cin >> N;
 
   double **Array = new double*[N];
   for( i = 0; i < N; i++ )
   {
      Array[i] = new double[N];
   };
 
   double min;
   cout << "Threshold value is " << endl;
   cin >> min;
 
 
 
   for(i = 0; i < N; i++)
   {
      for(j = 0; j < N; j++)
      {
        // Генерация элементов в диапазоне от -250 до 250
        Array[j][i] = (rand()%31 - 15);
        cout << "a[" << j+1 <<"][" << i+1<< "]= " <<Array[j][i]<<'\t';
        sum += Array[j][i]/N;
      }
      cout << endl;
      if (sum <= min)
         kol++;
      sum = 0;
   }
 
   for(i = 0; i < N-1; i++)
   {
      for( j = i; j < N; j++)
      {
         if(Array[i][i] != 0)
         {
            for(int k = 0; k < N; k++)
            {
[B][COLOR="Black"]               Array[j][k] = (long double)Array[j][k] *(1-(Array[j][i])*1.0/(Array[i][i]));[/COLOR][/B]
            }
         }
      }
   }
   for(i = 0; i < N; i++)
   {
      for(j = 0; j < N; j++)
      {
        cout << "a[" << j+1 <<"][" << i+1<< "]= " <<Array[j][i]<<'\t';
      }
      cout << endl;
   }
 
   for(i = 0; i < N; i++ ) {
        delete [] Array[i];
   }
   delete [] Array;
 
   getchar();
   getchar();
 
   return 0;
}
 
 
//---------------------------------------------------------------------------

   for(i = 0; i < N-1; i++)
   {
      for( j = i+1; j < N; j++)
      {
         if(Array[i][i] != 0)
         {
            long double koef = Array[i][j]/Array[i][i];
            for(int k = 0; k < N; k++)
            {
               Array[k][j] -= Array[k][i]*koef;
            }
         }
      }
   }