module computing cosmological functions
angular diameter distance in units of h^-1 Mpc
volume element in units of h^-3 Mpc^3
E=H(z)/H0
coordinate distance , in units of h^-1 Mpc
integrand for coordinate distance