Calculate and plot Hubble diagrams of Supernovae Ia (i.e their distance moduli as a function of their cosmological redshifts) for the following three flat (Omega_k = 0) cosmological models with H_0=71 km/s/Mpc: realistic (Omega_M = 0.3), matter dominated (Omega_M = 1) and dark energy dominated (Omega_M = 0), and compare them with the observed cosmological sample of SN Ia from the Table 5 from Riess et al. 2004, ApJ, 607, 665. Which cosmological model gives the best agreement?