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

Question

我有这个等式：

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

应该是这样的：

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

标签： pythonsympy

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]