/* nag_elliptic_integral_rf (s21bbc) Example Program.
 *
 * Copyright 1990 Numerical Algorithms Group.
 *
 * Mark 2 revised, 1992.
 *
 * Mark 3 revised, 1994.
 */

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

int main(void)
{
  double rf, x, y, z;
  Integer ix;

  Vprintf("nag_elliptic_integral_rf (s21bbc) Example Program Results\n");
  Vprintf("     x      y      z       s21bbc\n");
  for (ix=1; ix<=3; ix++)
    {
      x = ix*0.5;
      y = (ix+1)*0.5;
      z = (ix+2)*0.5;
      /* nag_elliptic_integral_rf (s21bbc).
       * Symmetrised elliptic integral of 1st kind R_F(xyz)
       */
      rf = nag_elliptic_integral_rf(x, y, z, NAGERR_DEFAULT);
      Vprintf("%7.2f%7.2f%7.2f%12.4f\n", x, y, z, rf);
    }
  return EXIT_SUCCESS;
}