Plot 在 python 中具有可变的积分限制

Question

我有 function f(x) 并且我想将 function 从 0 积分到区间 (0,2) 中的某个点 x。 我知道如何在数学中解决这个问题，但是我不知道如何在 python 中解决这个问题。

数学

OM=0.3
OL=0.7
f[x_] := ((1 + x)^3 OM + OL)^(-1/2);
Plot[NIntegrate[f[x], {x, 0, z2}], {z2, 0, 2}]

我在python中的代码：

import numpy as np
import matplotlib.pyplot as plt
from scipy import integrate

y = np.linspace(0, 2, 20)
Om=0.3
Ol=1-Om

def H(x):
    return (Om*(1+x)**3 +Ol)**(-1./2)

#Integral from a single point in the inverval
print quad(H,0,2)

我怎样才能在 python 中做同样的事情？

Answer 1

[integrate.quad(H, 0, t)[0] for t in y]

quad的输入是f, a, b ，在这里您对y中的所有t进行从0到t的积分。