scipy.special.it2j0y0#

scipy.special.it2j0y0(x, выход=None) = 'it2j0y0'>#

Интегралы, связанные с функциями Бесселя первого рода порядка 0.

Вычисляет интегралы

\[\begin{split}\int_0^x \frac{1 - J_0(t)}{t} dt \\ \int_x^\infty \frac{Y_0(t)}{t} dt.\end{split}\]

Подробнее о \(J_0\) и \(Y_0\) см. j0 и y0.

Параметры:

xarray_like: Значения, в которых вычисляются интегралы.
выходкортеж ndarrays, опционально: Необязательные выходные массивы для результатов функции.

Возвращает:

ij0скаляр или ndarray: Интеграл для j0
iy0скаляр или ndarray: Интеграл для y0

Ссылки

[1]

S. Zhang и J.M. Jin, «Computation of Special Functions», Wiley 1996

Примеры

Вычислить функции в одной точке.

>>> from scipy.special import it2j0y0
>>> int_j, int_y = it2j0y0(1.)
>>> int_j, int_y
(0.12116524699506871, 0.39527290169929336)

Вычислить функции в нескольких точках.

>>> import numpy as np
>>> points = np.array([0.5, 1.5, 3.])
>>> int_j, int_y = it2j0y0(points)
>>> int_j, int_y
(array([0.03100699, 0.26227724, 0.85614669]),
 array([ 0.26968854,  0.29769696, -0.02987272]))

Построить графики функций от 0 до 10.

>>> import matplotlib.pyplot as plt
>>> fig, ax = plt.subplots()
>>> x = np.linspace(0., 10., 1000)
>>> int_j, int_y = it2j0y0(x)
>>> ax.plot(x, int_j, label=r"$\int_0^x \frac{1-J_0(t)}{t}\,dt$")
>>> ax.plot(x, int_y, label=r"$\int_x^{\infty} \frac{Y_0(t)}{t}\,dt$")
>>> ax.legend()
>>> ax.set_ylim(-2.5, 2.5)
>>> plt.show()

../../_images/scipy-special-it2j0y0-1.png