我使用比率测试来查找元素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的正确限制?