/* nag_shapiro_wilk_test(g01ddc) Example Program.
 *
 * Copyright 1996 Numerical Algorithms Group.
 *
 * Mark 4, 1996.
 *
 */

#include <nag.h>
#include <nag_stdlib.h>
#include <stdio.h>
#include <nagg01.h>
#include <nagm01.h>

int main(void)
{

#define NMAX 20

  /* Local variables */
  double a[NMAX];
  Integer i, j, n;
  double w, x[NMAX], pw;
  Boolean calwts;

  Vprintf("g01ddc Example Program Results\n");

  /* Skip heading in data file */
  Vscanf("%*[^\n] ");

  calwts = TRUE;
  Vscanf("%ld ", &n);
  if (n > 0 && n <= NMAX)
    {
      for (j = 1; j <= 2; ++j)
        {
          for (i = 1; i <= n; ++i)
            Vscanf("%lf ", &x[i - 1]);

          m01cac(x, (size_t)n, Nag_Ascending,NAGERR_DEFAULT);
          g01ddc(n, x, calwts, a, &w, &pw, NAGERR_DEFAULT);

          Vprintf("\n For sample number %2ld, value of W statistic = %7.4f\n",
                  j,w);
          Vprintf("                       Significance level is %8.4f\n", pw);
          calwts = FALSE;
        }
    }
  return EXIT_SUCCESS;
}