/* nag_regress_confid_interval(g02cbc) Example Program
 *
 * Copyright 1992 Numerical Algorithms Group.
 *
 * Mark 3, 1992.
 */

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

#define NMAX 10

int main(void)
{
  Nag_SumSquare mean;
  Integer n;
  double x[NMAX], y[NMAX], wt[NMAX];
  double clm, clp;
  double yhat[NMAX], yml[NMAX], ymu[NMAX], yl[NMAX], yu[NMAX], h[NMAX],
  res[NMAX], rms;
  Integer i;
  char m, w;

  Vprintf("g02cbc Example Program Results\n");
  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vscanf("%ld\n",&n);
  if (n>=1 && n<= NMAX)
    {
      Vscanf("%lf%lf\n",&clm,&clp);
      Vscanf(" %c %c\n",&m,&w);
      if (m == 'm' || m == 'M')
        mean = Nag_AboutMean;
      else if (m == 'z'|| m == 'Z')
        mean = Nag_AboutZero;
      if (w == 'w' || w == 'W')
        for (i=0; i<n; i++)
          Vscanf("%lf%lf%lf\n",&x[i],&y[i],&wt[i]);
      else
        for (i = 0; i < n; ++i)
          Vscanf("%lf%lf\n",&x[i],&y[i]);

      g02cbc(mean, n, x, y, wt, clm, clp, yhat, yml, ymu, yl, yu, h, res,
             &rms, NAGERR_DEFAULT);

      Vprintf("\ni       yhat[i]    yml[i]    ymu[i]      yl[i]      yu[i]"
              " \n\n");
      for (i=0; i < n; ++i) {
        Vprintf("%ld %10.2f %10.2f", i, yhat[i], yml[i]);
        Vprintf("%10.2f  %10.2f %10.2f\n",ymu[i], yl[i], yu[i]);
      }
    }
  else
    {
      Vfprintf(stderr, "n is out of range:"
               " n = %-3ld\n",n);
      return EXIT_FAILURE;
    }
  return EXIT_SUCCESS;
}