代码中自变量名称错误,应将x改为r:
f=1/(sy.log(0.84/r))-120*r
f_root=sy.solve([f],[r])
改写成如下:
import sympy as sy
r=sy.symbols('r')
f=1/(sy.log(0.84/r))-120*r
f_root=sy.solve([f],[r])
print(f_root)
sy.plot(f,(r,-1,1))
fdiff1=f.diff(r,1)
fdiff2=f.diff(r,2)
# print(fdiff1)
# print(fdiff2)
phi=abs(f*fdiff2/fdiff1**2)
res=phi.subs({r:-1})
print(res)