/* nag_ref_vec_multi_normal(g05eac) Example Program
 *
 * Copyright 1991 Numerical Algorithms Group.
 *
 * Mark 2, 1991.
 *
 * Mark 3 revised, 1994.
 */

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

#define N 2
#define TDC N

int main(void)
{
  Integer i, j;
  double a[N], c[N][TDC], z[N];
  double *r = (double *)0;
  double eps = 0.01;

  Vprintf("g05eac Example Program Results\n");
  a[0] = 1.0;
  a[1] = 2.0;
  c[0][0] = 2.0;
  c[1][1] = 3.0;
  c[0][1] = 1.0;
  c[1][0] = 1.0;
  g05cbc((Integer)0);
  g05eac(a, (Integer)N, &c[0][0], (Integer)TDC,
         eps, &r, NAGERR_DEFAULT);
  for (i=1; i<=5; i++)
    {
      g05ezc(z, r);
      for (j=0; j<2; j++)
        Vprintf("%10.4f",z[j]);
      Vprintf("\n");
    }
  NAG_FREE(r);
  return EXIT_SUCCESS;
}