2

我无法将它集成到 python (1/1+t))+(-1/exp(t)) [0,np.inf]

import numpy as np
from math import exp
from scipy.integrate import quad

print(quad(lamba t: (1/1+t))+(-1/exp(t)),0,np.inf)

它必须显示欧拉常数

这就是我想要整合的

4

2 回答 2

1

顺便说一句,有了这个积分,你永远不会有欧拉常数。

欧拉常数定义为(抱歉我不能发布图片):

http://imgur.com/W5brbW4

因此,您必须将lambda功能更改为以下内容:

import numpy as np
from math import exp
from scipy.integrate import quad
f = lambda t: 1/t * (1/(1+t) - exp(-t))

In [21]: quad(f, 0.0, np.inf)
Out[21]: (0.5772156649015537, 3.613579096292482e-10)

结果就是0.5772156649015537欧拉常数的值。

于 2016-12-29T11:18:37.177 回答
0

你写lamba的不是lambda,这应该引发一个SyntaxError.

此外,阅读您的方程式,您传递的函数应该是:

lambda t: 1/(1+t) - exp(-t)

尝试将最后一行更改为:

print(quad(lambda t: (1/(1+t) - exp(-t)), 0, np.inf))
于 2016-12-29T11:06:05.967 回答