/* nag_real_eigenvalues(f02afc) Example Program
 *
 * Copyright 1990 Numerical Algorithms Group.
 *
 * Mark 2 revised, 1992.
 */

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

#define NMAX 4
#define TDA NMAX
#define COMPLEX(A) A.re, A.im

int main(void)
{
  Integer  i, j, n;
  double a[NMAX][TDA];
  Complex r[NMAX];
  Integer iter[NMAX];

  Vprintf("f02afc Example Program Results\n");
  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vscanf("%ld",&n);
  if (n<1 || n>NMAX)
    {
      Vfprintf(stderr, "N is out of range: N = %5ld\n", n);
      return EXIT_FAILURE;
    }
  for (i=0; i<n; i++)
    for (j=0; j<n; j++)
      Vscanf("%lf", &a[i][j]);
  f02afc(n, &a[0][0], (Integer)TDA, r, iter, NAGERR_DEFAULT);
  Vprintf("Eigenvalues\n");
  for (i=0; i<n; i++)
    Vprintf("( %7.3f  , %7.3f ) \n", COMPLEX(r[i]));
  return EXIT_SUCCESS;
}