/* nag_incomplete_gamma (s14bac) Example Program.

 *

 * Copyright 1990 Numerical Algorithms Group.

 *

 * Mark 2 revised, 1992.

 */



#include <nag.h>

#include <stdio.h>

#include <nag_stdlib.h>

#include <nags.h>

#include <nagx02.h>



int main(void)

{

  double a, p, q, tol, x;



  /* Skip heading in data file */

  Vscanf("%*[^\n]");

  Vprintf("nag_incomplete_gamma (s14bac) Example Program Results\n");

  /* nag_machine_precision (x02ajc).

   * The machine precision

   */

  tol = X02AJC;

  Vprintf("        a           x           p           q\n");

  while (scanf("%lf %lf", &a, &x) != EOF)

    {

      /* nag_incomplete_gamma (s14bac).

       * Incomplete Gamma functions P(ax) and Q(ax)

       */

      nag_incomplete_gamma(a, x, tol, &p, &q, NAGERR_DEFAULT);

      Vprintf("%12.4f%12.4f%12.4f%12.4f\n", a, x, p, q);

    }

  return EXIT_SUCCESS;

}