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

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

#define NMAX 8
#define TDA NMAX

main()
{
  
  double   a[NMAX][TDA], b[NMAX], x[NMAX];
  Integer  i, j, n;
  
  Vprintf("f04arc 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);
      exit(EXIT_FAILURE);
    }
  for (i=0; i<n; i++)
    for (j=0; j<n; j++)
      Vscanf("%lf",&a[i][j]);
  for (i=0; i<n; i++) 
    Vscanf("%lf",&b[i]);
  f04arc(n,(double *)a, (Integer)TDA, b, x, NAGERR_DEFAULT);
  Vprintf("Solution\n");
  for (i=0; i<n; i++) 
    Vprintf("%9.4f\n",x[i]);
  exit(EXIT_SUCCESS);
}