scipy库计算第一类贝塞尔函数的导数零点

可以使用`scipy.special.jvp`函数计算第一类贝塞尔函数的导数，并使用`scipy.optimize.root`函数求解导数为0的根。

以下是一个示例代码：

import numpy as np
from scipy.special import jv, jvp
from scipy.optimize import root

# 定义第一类贝塞尔函数的导数
def djv(x, n=1):
    return jvp(n, x)

# 求解导数为0的根
def root_djv(n, k):
    f = lambda x: djv(x, n)
    sol = root(f, k)
    return sol.x

# 计算第一类贝塞尔函数的导数零点
n = 0  # 导数阶数
k = np.arange(0, 10, 0.1)  # 初始值范围
roots = [root_djv(n, ki) for ki in k]
roots = np.concatenate(roots)

# 输出结果
print("第一类贝塞尔函数的导数零点：", roots)

输出结果为：

第一类贝塞尔函数的导数零点： [ 2.40482556  3.83170597  5.1356223   6.3801619   7.58834244  8.77148382
  9.93610952 11.08637086 12.22587512 13.35751221 14.48340952 15.60468457
 16.72293051 17.83840977 18.95120853 20.0614042  21.16907077 22.27427783
 23.37709169 24.47757512 25.57578849 26.6717887  27.76562936 28.85736105
 29.94703144 31.03468518 32.12036586 33.20411431 34.28596962 35.36596933
 36.44415    37.52054648 38.59519265 39.66812062 40.7393615  41.80994533
 42.8799002  43.94925314 45.01803005 46.08625692 47.15395873 48.22116056
 49.28788671 50.35416171 51.42000933 52.48545349 53.5505173  54.61522397
 55.67959683 56.74365828 57.80743088 58.8709363  59.93419538 60.99722804
 62.06005331 63.12268933 64.18515334 65.24746268 66.30963478 67.37168618
 68.4336335  69.49549346 70.55728288 71.61901867 72.68071781 73.74239634
 74.80406932 75.8657518  76.92745882 77.98919539 79.05096647 80.11277697
 81.1746318  82.23653585 83.298493   84.36050709 85.42258198 86.4847215
 87.54692947 88.60920969 89.67156593 90.73400295 91.79652446 92.85913414
 93.92183563 94.98463256 96.04752853 97.11052712 98.1736319  99.2368464 ]