PROGRAM Test
IMPLICIT NONE
REAL, PARAMETER :: Omega_m0 = 0.27, Omega_rel0 = 8.24E-5, Omega_L0 = 0.73, Omega_0 = 1.02
REAL :: z_max, z, dz
REAL :: I, g_z, g_zplus, S
REAL :: d_p0, d_L, d_A, d_nr, d_r
INTEGER :: n = 100000
OPEN(UNIT = 10)
z_max = 0
DO WHILE (z_max <= 4.001)
dz = z_max/n...