﻿功能:第三类Carlson对称椭圆积分

格式:EllipticRJ(x, y, z, r)

x: 一个数值不小于0的数值或者矩阵变量
y: 一个数值不小于0的数值或者矩阵变量
z: 一个数值不小于0的数值或者矩阵变量
r: 数值或者矩阵变量

说明:
$$
\begin{aligned}
\mathbf{EllipticRJ}(x, y, z, r) &= \dfrac{3}{2}\int_0^\infty \dfrac{1}{(t + r)\sqrt{(t + x)(t + y)(t + z)}}dt
\\
\\
\mathbf{s}\cdot\mathbf{t}\cdot&:x\geq 0,y\geq 0,z\geq 0
\end{aligned}
$$
参考: Carlson B C . Computing elliptic integrals by duplication[J]. Numerische Mathematik, 1979, 33(1):1-16.
参考:{Carlson Integral<https://mathworld.wolfram.com/CarlsonEllipticIntegrals.html>}

例子:

y = EllipticRJ(1.58,2.36,0.78,1.78) //回车得到如下结果
y =
[ 0.49866476551760 ]
