REAL eps, eta, x;
INTEGER ifail, k, nfmax;
REAL, EXTERNAL :: f

WRITE(*,*)"C05AJF/E Example Program Results"

  DO k=1,2

      IF (K == 1 ) THEN
         eps =  0.1e-3 
      ELSE
         eps = 0.1e-4;
      END IF
      x = 1.0;
      eta = 0.0;
      nfmax = 200;
      ifail = 1;
!      /* Remove underscore as appropriate */
      CALL c05aje(x,eps,eta,f,nfmax,ifail);
      IF (ifail==0) THEN
        WRITE(*,*) "With eps = ", eps, " root = ", x
      ELSE
        WRITE(*,*) "ifail = ", ifail
          IF (ifail==3 .OR. ifail==4)  THEN
            WRITE(*,*) "With eps = ", eps, " final value = ", x
          END IF
      END IF
  END DO
END PROGRAM

FUNCTION  f(x)
F = exp(-x) - x;
RETURN
END FUNCTION f