def print_flat_cosmo(): """Print cosmological functions for a flat cosmology.""" H0 = 70. Om0 = 0.3 cosmo = FlatLambdaCDM(H0=H0, Om0=Om0, Tcmb0=0.) print('===============') print('Flat cosmology:') print('===============') print('H0={} [km/s/Mpc], Om0={}, Ode0={}, Orel0={}, Ok0={}\n'.format( H0, Om0, cosmo.Ode(0.), cosmo.Ogamma(0.) + cosmo.Onu(0.), cosmo.Ok(0.))) redshifts = (0., 0.1, 10., 1000.) print_curves(redshifts, cosmo)