Python SciPy special.itmodstruve0用法及代碼示例

本文簡要介紹 python 語言中 scipy.special.itmodstruve0 的用法。

用法:

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

0 階修正 Struve 函數的積分。

\[I = \int_0^x L_0(t)\,dt\]

參數 ：：

x： array_like

積分上限(浮點數)。

out： ndarray，可選

函數值的可選輸出數組

返回 ：：

I： 標量或 ndarray

積分為\(L_0\)從 0 到x.

注意：

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 itmodstruve0
>>> itmodstruve0(1.)
0.3364726286440384

通過提供 x 的數組，在多個點評估該函數。

>>> points = np.array([1., 2., 3.5])
>>> itmodstruve0(points)
array([0.33647263, 1.588285  , 7.60382578])

繪製 -10 到 10 之間的函數。

>>> import matplotlib.pyplot as plt
>>> x = np.linspace(-10., 10., 1000)
>>> itmodstruve0_values = itmodstruve0(x)
>>> fig, ax = plt.subplots()
>>> ax.plot(x, itmodstruve0_values)
>>> ax.set_xlabel(r'$x$')
>>> ax.set_ylabel(r'$\int_0^xL_0(t)\,dt$')
>>> plt.show()