/* nag_tsa_auto_corr(g13abc) Example Program
 *
 * Copyright 1991 Numerical Algorithms Group.
 *
 * Mark 2, 1991.
 */

#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nagg13.h>

#define NXMAX 50
#define NKMAX 10

int main(void)
{

  double mean, xv, stat;
  Integer i, nk, nx;
  double r[NKMAX], x[NXMAX];

  Vprintf("g13abc Example Program Results\n");
  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vscanf("%ld %ld", &nx, &nk);

  if (nk>0 && nk<=NKMAX && nx>0 && nx<=NXMAX)
    {
      for (i=0; i<nx; ++i)
        Vscanf("%lf", &x[i]);
      Vprintf("\nThe first %2ld coefficients are required\n", nk);

      g13abc(x, nx, nk, &mean, &xv, r, &stat, NAGERR_DEFAULT);

      Vprintf("The input array has sample mean %12.4f\n", mean);
      Vprintf("The input array has sample variance %12.4f\n", xv);
      Vprintf("The sample autocorrelation coefficients are\n\n");
      Vprintf("   Lag    Coeff\n");
      for (i=0; i<10; ++i)
        Vprintf("%6ld%10.4f\n", i+1, r[i]);
      Vprintf("\nThe value of stat is %12.4f\n", stat);
    }
  else
    {
      Vfprintf(stderr, "One or both of nk and nx are out of range:"
               " nk = %-3ld while  nx = %-3ld\n", nk, nx);

      return EXIT_FAILURE;
    }
  return EXIT_SUCCESS;
}