/* nag_cos_integral (s13acc) Example Program.

 *

 * Copyright 1990 Numerical Algorithms Group.

 *

 * Mark 2 revised, 1992.

 */



#include <nag.h>

#include <stdio.h>

#include <nag_stdlib.h>

#include <nags.h>



int main(void)

{

  double x, y;



  /* Skip heading in data file */

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

  Vprintf("nag_cos_integral (s13acc) Example Program Results\n");

  Vprintf("     x           y\n");

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

    {

      /* nag_cos_integral (s13acc).

       * Cosine integral Ci~(x)

       */

      y = nag_cos_integral(x, NAGERR_DEFAULT);

      Vprintf("%12.3e%12.3e\n", x, y);

    }

  return EXIT_SUCCESS;

}