SymPy渐近线怪癖及其解决方法

Question

今天，我偶然发现了SymPy的一个奇怪功能。 假设我们尝试在渐近线处评估一个函数：

>>from sympy import *
>>
>>x = Symbol("x")
>>N(cot(x),subs={x: 0})
cot(x)

它既不返回NaN，也不返回错误或异常，而是函数本身！
我的问题：
1.如何自动检测到我偶然发现了一个渐近线（除了在结果上调用type（）之外）？
2.此设计决定背后的动机是什么？

Answer 1

在表达式上使用subs而不是N的参数：

In [3]: N(cot(x).subs(x,0))
Out[3]: zoo

这可能是sympy中的错误

Answer 2

要回答问题1，您可以自动尝试将结果转换为float() ，如下所示：

try:
    result=float(N(cot(x),subs={x: 0}))
except Exception as e:
    plus=limit(cot(x), x, 0, '+')
    minus=limit(cot(x), x, 0, '-')
    func=str(inspect.trace()[-1][0].f_locals['result'])
    if plus==minus:
        print("Function %s has a point discontinuity at the point specified. An approximation of the function taking the limit of both sides is %s" % (func,plus))
    else:
        print("Function %s has a jump discontinuity at the point specified" % func)

（感谢有关检查当地人的问题 ）