python - 如何求解积分的极限？

Question

我有这个等式:

integrate(exp(-x**2), (x, c, -c))=1/2

应该是这样的:

如何找到 c 的解决方案，例如使用 sympy 等库？

Answer 1

使用 Sympy，您可以执行以下操作:

from sympy import *
from sympy.abc import x, c

sol = solve(Eq(integrate(exp(-x**2), (x,c,-c)), S(1)/2))
print("symbolic solution:", sol)
print("numeric solution:", [s.evalf() for s in sol])

输出:

symbolic solution: [erfinv(-1/(2*sqrt(pi)))]
numeric solution: [-0.255449286541000]