/* nag_elliptic_integral_rd (s21bcc) 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 rd, x, y, z;
  Integer ix, iy;

  Vprintf("s21bcc Example Program Results\n");
  Vprintf("      x      y      z       s21bcc  \n");
  for (ix=1; ix<=3; ix++)
    {
      x = ix*0.5;
      for (iy=ix; iy<=3; iy++)
        {
          y = iy*0.5;
          z = 1.0;
          rd = s21bcc(x, y, z, NAGERR_DEFAULT);
          Vprintf(" %7.2f%7.2f%7.2f%12.4f\n", x, y, z, rd);
        }
    }
  return EXIT_SUCCESS;
}