Elliptic-integral-rf-func
calculates an approximation to the integral
where x, y, z ≥ 0 and at most one is zero.
For more information please review the s21bbc function in the NAG document.