python - 在 python 上集成

Question

嗨，我的讲师向我提出了一个问题，要求通过 python 集成一个函数，但他给了我们很少的信息。边界是+无穷大和-无穷大，函数是

(cos a*x) * (e**-x**2)

到目前为止我已经

def gauss_cosine(a, n):
    sum=0.0
    dx = ((math.cosine(a*x)*math.exp(-x**2)))
    return 
    for k in range (0,n):
        x=a+k*dx
        sum=sum+f(x)
    return dx*sum

不确定这是否正确。亲切的问候

Answer 1

我在这个网站上没有看到太多推荐，但你可以尝试 sympy:

In [1]: import sympy as sp
In [2]: x, a = sp.symbols(('x', 'a'))
In [3]: f = sp.cos(a*x) * sp.exp(-x**2)
In [4]: res = sp.integrate(f, (x, -sp.oo, sp.oo))
In [5]: res
Out[5]: sqrt(pi) * exp
In [6]: sp.pprint(res)
Out[6]:
         2 
       -a  
       ────
  ___   4  
╲╱ π ⋅ℯ    

对于数值积分，请尝试使用 scipy 包。