Elliptic Integral Third Kind RJ

      x        y        z        p           actual          rj(x,y,z,p)
    0.10     0.10     0.10     0.10       3.162278e+01       3.162278e+01
    2.00     3.00     4.00     5.00       1.429758e-01       1.429758e-01
  100.00   100.00   100.00   100.00       1.000000e-03       1.000000e-03