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

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

用法:

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

第一類 0 階貝塞爾函數的積分。

計算積分

\[\begin{split}\int_0^x J_0(t) dt \\ \int_0^x Y_0(t) dt.\end{split}\]

有關 \(J_0\) 和 \(Y_0\) 的更多信息，請參閱 j0 和 y0 。

參數 ：：

x： array_like

評估積分的值。

out： ndarrays 元組，可選

函數結果的可選輸出數組。

返回 ：：

ij0： 標量或 ndarray

j0 的積分

iy0： 標量或 ndarray

y0 的積分

參考：

[1]

S. 張和 J.M. Jin，“特殊函數的計算”，Wiley 1996

例子：

在某一點評估函數。

>>> from scipy.special import itj0y0
>>> int_j, int_y = itj0y0(1.)
>>> int_j, int_y
(0.9197304100897596, -0.637069376607422)

在幾個點上評估函數。

>>> import numpy as np
>>> points = np.array([0., 1.5, 3.])
>>> int_j, int_y = itj0y0(points)
>>> int_j, int_y
(array([0.        , 1.24144951, 1.38756725]),
 array([ 0.        , -0.51175903,  0.19765826]))

繪製從 0 到 10 的函數。

>>> import matplotlib.pyplot as plt
>>> fig, ax = plt.subplots()
>>> x = np.linspace(0., 10., 1000)
>>> int_j, int_y = itj0y0(x)
>>> ax.plot(x, int_j, label="$\int_0^x J_0(t)\,dt$")
>>> ax.plot(x, int_y, label="$\int_0^x Y_0(t)\,dt$")
>>> ax.legend()
>>> plt.show()