function y=zD(z) % z-distance relation, distance in units of Mpc h=0.72; DH=3000./h; % in Mpc OmegaM=0.274; OmegaLam=0.726; %fun=strcat('',num2str(OmegaM),'*(1+x)^3+',num2str(OmegaLam),'') DC=DH.*comdis(OmegaM,OmegaLam,z); DM=DC; % for zero curvature y=DM.*(1+z);