最初,我有两个数组,分别对应一个函数中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))
你的函数是一个函数 x
(在给定值为 x
吐出一个值),所以首先你应该把它重新包装成这样(引入一个函数 yy
其中,给定x,产生必要的 y
),然后写上 LL(x) = L(x, yy[x])
,然后用 scipy.integrate
来整合它。