Elliptic-integral-rd-func
calculates an approximate value for the integral
where x, y ≥ 0, at most one of x and y is zero, and z > 0.
For more information please review the s21bcc function in the NAG document.