Elliptic-integral-rj-func
calculates an approximation to the integral
Where x,y,z ≥ 0 , ≠0 and at most one of x, y and z is zero.
If <0, the result computed is the Cauchy principal value of the integral.
For more information please review the s21bdc function in the NAG document.