功能: 计算Bulirsch第1类不完全椭圆积分
格式: EllipticBulirschEl1(x, kc)
x  : 数值或者矩阵变量
kc : 数值或者矩阵变量


说明:
1、函数定义如下
$$
\mathbf{EllipticBulirschEl1}(x, k_c) = \int_0^x \dfrac{1}{\sqrt{(1+t^2)(1+k_c^2t^2)}}dt
$$
参考:
[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 = EllipticBulirschEl1(1,2)//执行得到如下结果
a =
[ 0.65264594921422 ]
