python - 有没有一种简单的方法可以通过 sympy 获取某个域中函数的所有不连续性？

Question

给定一个 sympy 表达式，有没有办法找到给定区间内的所有不连续点？例如，给定从 -2 到 2 的 1/(x^2-1)，它将返回 -1 和 1。它不必是符号的。数值解实际上可能更适合我的目的。

Answer 1

为此，您可以使用 singularities 模块。

In [ ]: from sympy import *  
In [ ]: init_printing()
In [ ]: x = symbols('x')
In [ ]: singularities(1/(x**2 - 1), x)
Out[ ]: (-1, 1) # A tuple of SymPy objects

引用:http://docs.sympy.org/latest/modules/calculus/index.html#sympy.calculus.singularities.singularities