Elliptic Integral Second Kind RD x y z actual rd(x,y,z) 0.10 0.10 0.10 3.162278e+01 3.162277e+01 0.00 2.00 1.00 1.797210e+00 1.797210e+00 100.00 100.00 100.00 1.000000e-03 1.000000e-03