static const uint32_t N_RUNS = 5;

  // We want the number of observations in each bin to be > 5

  static const uint32_t N_BINS = 50;

  // The following values should yield many more than 5 per bin

  static const uint32_t N_BINS = 100;

  double sum = 0.;

  double confidence = 0.99;

  double maxStatistic = gsl_cdf_chisq_Qinv (0.05, N_BINS);

  double maxStatistic = gsl_cdf_chisq_Pinv (confidence, (N_BINS-1));

  NS_LOG_DEBUG ("Chi square required at " << confidence << " confidence for " << N_BINS << " bins is " << maxStatistic);

  double result = maxStatistic;

  // If chi-squared test fails, re-try it up to N_RUNS times

      double result = ChiSquaredTest (u);

      result = ChiSquaredTest (u);

      sum += result;

      NS_LOG_DEBUG ("Chi square result is " << result);

      if (result < maxStatistic)

        {

          break;

        }

  sum /= (double)N_RUNS;

  NS_TEST_ASSERT_MSG_LT (result, maxStatistic, "Chi-squared statistic out of range");

  NS_TEST_ASSERT_MSG_LT (sum, maxStatistic, "Chi-squared statistic out of range");