python - 集成依赖于两个数组的函数

Question

最初，我有两个数组对应于一个函数中的 x 和 y 的值，但我不知道那个函数，我只知道 y 的值取决于 x。然后，我计算一个依赖于两个数组的函数。

我需要在 Python 中计算最后一个函数的积分，以获得 x 的第一个值和最后一个值之间曲线下的总面积。知道如何做到这一点吗？

x = [array]
y(x) = [array]

a = 2.839*10**25
b = 4*math.pi
alpha = 0.5
z = 0.003642
def L(x,y,a,b,alpha,z):
    return x*((y*b*a)/(1+z)**(1+alpha))

Answer 1

你的函数是 x 的函数(给定 x 的值，它会吐出一个值)，所以首先你应该这样重新打包它(引入一个函数 yy，给定 x，生成必需的 y)，然后编写 LL(x) = L(x, yy[x]) , 然后使用 scipy.integrate 来整合它。