/* nag_elliptic_integral_rj (s21bdc) 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 r, rj, x, y, z;
  Integer ix, iy, iz;

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