SymPy未通过限制测试

时间:2016-02-16 09:30:29

标签: sympy

我使用比率测试来查找元素a = n**2/3**n的系列是否收敛。在纸上,很容易看出比率a(n + 1)/ a(n) - > 1/3,但SymPy似乎在这个测试中失败了,

>>> n = symbols('n', integer=True)
>>> a = n**2/3**n
>>> a1 = a.subs(n, n+1)
>>> r = a1/a
>>> r
3**n*3**(-n - 1)*(n + 1)**2/n**2
>>> simplify(r)
(n + 1)**2/(3*n**2)
>>> expand(r)
1/3 + 2/(3*n) + 1/(3*n**2)
>>> limit(r, n, oo)
1
>>> limit(expand(r), n, oo)
1/3

为什么SymPy不能评估比率a1 / a的正确限制?

0 个答案:

没有答案