/* nag_real_jacobian_elliptic (s21cac) Example Program
 *
 * Copyright 2002 Numerical Algorithms Group.
 *
 * Mark 7, 2002.
 */

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

int main(void)
{
  double u, m, sn, cn, dn;

  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vprintf("s21cac Example Program Results\n");
  Vprintf("      u            m            sn           cn           dn\n");
  while (scanf("%lf %lf", &u, &m) != EOF)
    {
      s21cac(u, m, &sn, &cn, &dn, NAGERR_DEFAULT);
      Vprintf("%12.3e %12.3e %12.3e %12.3e %12.3e\n", u, m, sn, cn, dn);
    }
  return EXIT_SUCCESS;
}