记录下泰勒展开式和极限的代码实现,多看官方文档,幸福长久远。
以如下函数为例:
import sympy
"""自定义函数,求其在x=0附近的泰勒展开(4项)"""
x = sympy.symbols('x')
fx = (x - sympy.sin(x)) / x ** 3
"""级数展开:series(函数表达式, x0, n),使用.removeO()去除皮亚诺余项,remove后面的o大写"""
taylor_expansion_1 = fx.series(x, 0, 5)
taylor_expansion_2 = fx.series(x, 0, 5).removeO()
print(taylor_expansion_1)
print(taylor_expansion_2)
"""求极限"""
print(sympy.limit(fx, x, sympy.oo)) # 趋于无穷
print(sympy.limit(fx, x, 0)) # 趋于0,默认值 dir='+',也就是趋于+0
print(sympy.limit(fx, x, 0, dir="-")) # 趋于0,默认值调整为dir='-',也就是趋于-0