Python SciPy special.it2struve0用法及代码示例

本文简要介绍 python 语言中 scipy.special.it2struve0 的用法。

用法:

scipy.special.it2struve0(x, out=None) = <ufunc 'it2struve0'>#

与 0 阶 Struve 函数相关的积分。

返回积分，

\[\int_x^\infty \frac{H_0(t)}{t}\,dt\]

其中\(H_0\) 是 0 阶 Struve 函数。

参数 ：：

x： array_like

积分的下限。

out： ndarray，可选

函数值的可选输出数组

返回 ：：

I： 标量或 ndarray

积分值。

注意：

Wrapper for a Fortran routine created by Shanjie Zhang and Jianming Jin [1].

参考：

[1]

张善杰和金建明。 “特殊函数的计算”，John Wiley and Sons，1996 年。https://people.sc.fsu.edu/~jburkardt/f_src/special_functions/special_functions.html

例子：

在某一点评估函数。

>>> import numpy as np
>>> from scipy.special import it2struve0
>>> it2struve0(1.)
0.9571973506383524

通过提供 x 的数组，在多个点评估该函数。

>>> points = np.array([1., 2., 3.5])
>>> it2struve0(points)
array([0.95719735, 0.46909296, 0.10366042])

绘制 -10 到 10 之间的函数。

>>> import matplotlib.pyplot as plt
>>> x = np.linspace(-10., 10., 1000)
>>> it2struve0_values = it2struve0(x)
>>> fig, ax = plt.subplots()
>>> ax.plot(x, it2struve0_values)
>>> ax.set_xlabel(r'$x$')
>>> ax.set_ylabel(r'$\int_x^{\infty}\frac{H_0(t)}{t}\,dt$')
>>> plt.show()