scipy.special.ellipkm1#
- scipy.special.ellipkm1(p, out=None) = <ufunc 'ellipkm1'>#
在 m = 1 附近的椭圆积分第一类
此函数被定义为
\[K(p) = \int_0^{\pi/2} [1 - m \sin(t)^2]^{-1/2} dt\]其中 m = 1 - p。
- 参数:
- p类似数组
将椭圆积分的参数定义为 m = 1 - p。
- outndarray,可选
函数值的可选输出数组
- 返回值:
- K标量或 ndarray
椭圆积分的值。
备注
Cephes [1] 例程 ellpk 的包装。
对于 p <= 1,计算使用近似,
\[K(p) \approx P(p) - \log(p) Q(p),\]其中 \(P\) 和 \(Q\) 是十阶多项式。实参 p 在内部使用,而不是 m,以便将 m = 1 处的对数奇点转移到原点;这会保留最大的精度。对于 p > 1,标识
\[K(p) = K(1/p)/\sqrt(p)\]被使用。
参考资料
[1]Cephes 数学函数库,http://www.netlib.org/cephes/