#include <math.h>
#include <stdio.h>
#include <nagmk19_t3e.h>

main()
{
  double eps, eta, x;
  int ifail, k, nfmax,nsig,out;
  double f(double *);
  FTCHAR title={"COMPLEX NUMBER:",15};
  COMPLEX cp;

  printf("C05AJF/E Example Program Results\n\n");
  for (k=1; k<=2; k++)
    {
      eps = k==1 ? 0.1e-3 : 0.1e-4;
      x = 1.0;
      eta = 0.0;
      nfmax = 200;
      ifail = 1;
 
      C05AJE(&x,&eps,&eta,f,&nfmax,&ifail);

      if (ifail==0)
        printf("With eps = %e   root = %f\n",eps,x);
      else
        {
          printf("ifail = %d\n",ifail);
          if (ifail==3 | ifail==4)
            printf("With eps = %e final value = %f\n",eps,x);
        }
     } 
  
  printf("\nTest of using complex numbers\n");
  cp=CMPLX(1,3); /* C version */
  k=1;nsig=3;out=6;
  Y07CBE(&k,&cp,&k,title,&nsig,&nsig,&k,&out);
}

double f(double *x)
{
  return exp(-*x) - *x;
}