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)
