/* nag_prod_limit_surviv_fn(g12aac) Example Program.
 *
 * Copyright 1996 Numerical Algorithms Group.
 *
 * Mark 4, 1996.
 *
 */

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

#define NMAX 18

int main(void)
{

  double psig[NMAX];
  double p[NMAX];
  double t[NMAX];
  double tp[NMAX];

  Integer i, n;
  Integer ifreq[NMAX], ic[NMAX], nd;

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

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

  Vscanf("%ld ", &n);
  if (n <= NMAX)
    {
      for (i = 0; i < n; ++i)
        Vscanf("%lf %ld %ld ", &t[i], &ic[i], &ifreq[i]);

      g12aac(n, t, ic, ifreq, &nd, tp, p, psig, NAGERR_DEFAULT);

      Vprintf("\n   Time   Survival    Standard\n");
      Vprintf("         probability  deviation\n\n");
      for (i = 0; i < nd; ++i)
        Vprintf(" %6.1f%10.3f  %10.3f\n", tp[i], p[i], psig[i]);
    }
  return EXIT_SUCCESS;

}