因此,我正在编写一个评估泰勒级数的脚本.但是,我希望它能够评估所有类型的功能.因此,我尝试使用例如acot(x)函数.

x = sy.Symbol('x')
f = acot(x)
...
func = taylor(f,0,3)
taylor_lambda = sy.lambdify(x, func, 'numpy')

上面的代码毫无例外地运行(例如,除非我使用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()