As I did in cases zArr and rzArr I want to use function dl(z) to create an array or list dlArr that contains the generated values from dl given z from zArr and rz from rzArr. I tried doing it like this but the syntax is wrong:
dlArr= [(dl(z,rz) for z in zArr & rz in rzArr)]
Relating code for reference:
zArr = np.linspace(0.01, 2.0, 1048)
def dist_integrand(z):
dist_integrand = 1.0 / np.sqrt(Omegam * (1 + z) ** 3 + Omegal)
return dist_integrand
def rz(z):
rz = integrate.quad(dist_integrand, 0, z)
return rz
rzArr = (rz(z) for z in zArr)
def dl(z,rz):
dl = (1+z) * rz
return dl
#!!
dlArr= [(dl(z,rz) for z in zArr & rz in rzArr)]