如何在python中使用积分方程的双积分?

时间:2015-07-06 23:45:37

标签: python scipy integration gaussian quad

我试图在Python中集成此功能:

I =积分(积分(g(y)* ln [f(x,y)/ g(y)] dx)dy) 其中g(y)=积分(f(x,y)dx)从0到1

eq = lambda x: x #NOTE: want to be able to be able to input any equation
f = lambda x,y: math.exp(-1.0*(eq(x)-y)**2)
g = lambda y: integrate.quad(f, 0, 1, args=(1,))[0]
combined = lambda x, y: g(y) * math.log(f(x,y) / g(y))
combInt = lambda y: integrate.quad(combined, 0,1, args=(1,))[0]
I = integrate.quad(combInt, -np.inf, np.inf)[0]

然而,当我运行我的代码时,我发现" g"是不变的,"我"是否定的。我是否错误地使用了四合一功能?怎么能不变?

0 个答案:

没有答案