/* nag_real_symm_eigenvalues(f02aac) Example Program
 *
 * Copyright 1990 Numerical Algorithms Group.
 *
 * Mark 1, 1990.
 */

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

#define NMAX 8
#define TDA NMAX


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

  Vprintf("f02aac 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]);
  f02aac(n, &a[0][0], (Integer)TDA, r, NAGERR_DEFAULT);
  Vprintf("Eigenvalues\n");
  for (i=0; i<n; i++)
    Vprintf("%9.4f%s",r[i],(i%8==7 || i==n-1) ? "\n": " ");
  return EXIT_SUCCESS;
}