Почему scipy поиск корня бесселя не возвращает корни в нуле?

Я пытаюсь использовать код, который использует нули функции Бесселя для других вычислений. Я заметил, что следующий фрагмент кода дает результаты, которые я считаю неожиданными.

    import scipy
    from scipy import special

    scipy.special.jn_zeros(1,2)

Я ожидаю, что результат этого вызова будет

    array([0., 3.83170597])

вместо

    array([3.83170597, 7.01558667])

Есть ли причина, по которой корень при x = 0,0 не возвращается?

Из того, что я вижу, корни симметричны вдоль оси x, за исключением тех, которые находятся в начале координат, но я не думаю, что это было бы достаточной причиной, чтобы полностью отказаться от корня.

На компьютере, который я использую, установлена ​​​​версия python 2.7.10 и используется версия 0.19.0 scipy.

P.S. следующая функция - это то, что я пытаюсь найти нули

    scipy.special.j1

python scipy bessel-functions
person WolcottR    schedule 02.05.2018    source источник

Ответы (1)


arrow_upward
2
arrow_downward

Похоже, принято не считать ноль при нуле, см., например, здесь. Может быть, это считается излишним?

person Paul Panzer    schedule 03.05.2018
comment
Это мое предположение, но, может быть, только для j0 нет корня в 0. Таким образом, я мог видеть там избыточность. В этом случае, оставив пользователю возможность обрабатывать один особый случай, это нормально. Благодарю вас! - person WolcottR; 03.05.2018