我尝试使用Sympy包计算f(x)= x ** 2 / /(x-1)的二阶导数,但垂直渐近线似乎好像x逼近1。如何消除呢?这是我的代码:
import sympy as sp
import numpy as np
import pylab as pyl
t = sp.symbols('t')
def fp(t):
return t**2 / (t -1)
f = sp.lambdify(t, fp(t), modules=['numpy'])
dp = sp.diff(fp(t))
d2p = sp.diff(dp)
d2 = sp.lambdify(t, d2p, modules=['numpy'])
x = np.linspace(-7.5, 10, 10**3)
pyl.plot(x, d2(x), 'b-')
pyl.ylim(-6, 10)
pyl.grid()
pyl.show()
答案 0 :(得分:0)
我会在同一通话中尝试二阶导数:
d2p = sp.diff(t**2 / (t-1), t, 2)
,看看是否有区别。您还可以评估d2p并查看d2p(1)
的值是什么。