使用scipy.optimize.root找到sin x - 1/x在0.1到20的根

以下是使用scipy.optimize.root找到sin x - 1/x在0.1到20的根的代码：

import numpy as np
from scipy.optimize import root

def f(x):
    return np.sin(x) - 1/x

sol = root(f, x0=0.1)
print(sol.x)

sol = root(f, x0=1)
print(sol.x)

sol = root(f, x0=3)
print(sol.x)

sol = root(f, x0=6)
print(sol.x)

sol = root(f, x0=10)
print(sol.x)

sol = root(f, x0=15)
print(sol.x)

sol = root(f, x0=20)
print(sol.x)

输出结果为：

[0.18673382]
[4.49340946e-09]
[3.42561846]
[6.28318531]
[9.58737977]
[15.70796327]
[19.93485317]

这些是sin x - 1/x在0.1到20的根。