功能: 计算Bulirsch完全椭圆积分
格式: EllipticBulirschCel(kc, p, a, b)
kc : 数值或者矩阵变量
p  : 数值或者矩阵变量
a  : 数值或者矩阵变量
b  : 数值或者矩阵变量

说明:
1、函数定义如下
$$
\begin{aligned}
\mathbf{EllipticBulirschCel}(k_c, p, a, b) &= \int_0^{\pi/2}\dfrac{a\cos^2(t)+b\sin^2(t)}{(\cos^2(t)+p\sin^2(t))\sqrt{\cos^2(t)+k_c^2\sin^2(t)}}dt
\\
\\
\mathbf{s}\cdot\mathbf{t}\cdot&:p\neq 0, k_c \neq 0
\end{aligned}
$$
参考:
[1] Bulirsch R. An extension of the Bartky-transformation to incomplete elliptic integrals of the third kind[J]. Numerische Mathematik, 1969, 13(3): 266-284.
[2] Carlson B C. Numerical computation of real or complex elliptic integrals[J]. Numerical Algorithms, 1995, 10(1): 13-26.
[3] NIST handbook of mathematical functions hardback and CD-ROM[M]. Cambridge university press, 2010.

例子:
a = EllipticBulirschCel(1,2,3,4)//执行得到如下结果
a =
[ 3.79223779587408 ]