FUNCTION zdist,z dh = 3000.0/0.7 omegak = -0.010 dm = dh*(1/sqrt(abs(omegak)))*sin((((sqrt(abs(omegak)))))*dc(z)/dh) return, dm/(1+z) end