我知道这里有一些类似的问题,但他们都没有真正解决我的问题。Python:绘制一个整体
我的代码如下所示:
import numpy
import matplotlib.pyplot as plt
from scipy import integrate as integrate
def H(z , omega_m , H_0 = 70):
omega_lambda=1-omega_m
z_prime=((1+z)**3)
wurzel=numpy.sqrt(omega_m*z_prime + omega_lambda)
return H_0*wurzel
def H_inv(z, omega_m , H_0=70):
return 1/(H(z, omega_m, H_0=70))
def integral(z, omega_m , H_0=70):
I=integrate.quad(H_inv,0,z,args=(omega_m,))
return I
def d_L(z, omega_m , H_0=70):
distance=(2.99*(10**8))*(1+z)*integral(z, omega_m, H_0=70)
return distance
的功能做的工作,我的问题:我如何可以绘制D_L与Z'就像这显然是一个问题,我在我的d_L的定义中有这个整数函数,它依赖于z和一些args =(omega_m,)。
你可以举一个'z'的例子,你如何在'd_L'中使用它? –
最后,z应该有足够的步数(500+)从0运行到2。 – user7248647