/* nag_1d_cheb_eval(e02aec) Example Program
 *
 * Copyright 1998 Numerical Algorithms Group.
 *
 * Mark 5, 1998.
 *
 */

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

int main(void)
{

#define NMAX 199
#define NP1MAX NMAX+1

  double xcap;
  double a[NP1MAX], p;

  Integer i, m, n;
  Integer r;

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

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

  while ((scanf("%ld",&m) != EOF))
    {
      if (m > 0)
        {
          Vscanf("%ld",&n);
          if (n >= 0 && n <= NMAX)
            {
              for (i = 0; i < n+1; ++i)
                  Vscanf("%lf",&a[i]);
              Vprintf("\n   R       Argument       Value of polynomial \n");
              for (r = 1; r <= m; ++r)
                {
                  xcap = (double) (2*r - m - 1) / (double) (m - 1);
                  e02aec(n+1, a, xcap, &p, NAGERR_DEFAULT);
                  Vprintf(" %3ld%14.4f    %14.4f\n",r,xcap,p);
                }
            }
          else
            {
              Vprintf( "Incorrect input value of n.\n");
              return EXIT_FAILURE;
            }
        }
    }
  return EXIT_SUCCESS;
}