因此,我正在编写一个评估泰勒级数的脚本。但是,我希望它能够评估所有类型的功能。因此,例如,我尝试使用功能acot(x)
。
x = sy.Symbol('x')
f = acot(x)
...
func = taylor(f,0,3)
taylor_lambda = sy.lambdify(x, func, 'numpy')
上面的代码毫无异常(exception)地运行(例如,除非我使用
acsch
,但它不会运行)。但是,当它到达这一行时:
plt.plot(x1,taylor_lambda(x1),label='taylor approximation')
我得到:
NameError: name 'acot' is not defined
我试图在lambdify调用中将
numpy
替换为sympy
,但这似乎具有象征意义。某些情况(更罕见的功能)正在发生这种情况,而其他情况则没有。谢谢!我的进口如下:
import sympy as sy
import numpy as np
from sympy.functions import *
from sympy import pi, E,acot
import matplotlib.pyplot as plt
import math
最佳答案
我之前确实有过类似的问题,并且已经设法解决了。
plt.plot(x1,taylor_lambda(x1),label='taylor approximation')
看起来还可以。我给我一个较旧的代码,可以正常工作,您可以进行比较。
from sympy.abc import x
from sympy import sin, series
from sympy.utilities.lambdify import lambdify
import numpy as np
import matplotlib.pyplot as plt
func = sin(x)/x
taylor = series(func, n=6).removeO()
evalfunc = lambdify(x, func, modules=['numpy'])
evaltaylor = lambdify(x, taylor, modules=['numpy'])
t = np.linspace(-7.5, 7.5 , 100)
plt.plot(t, evalfunc(t), 'b', label='sin(x)/x')
plt.plot(t, evaltaylor(t), 'r', label='Taylor')
plt.legend(loc='best')
plt.show()
关于python - (某些功能)未使用SymPy Lambdify定义,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/45273827/