/* nag_real_symm_eigensystem(f02abc) 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 8
#define TDA NMAX
#define TDV NMAX



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

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