Данные SN не сходятся с LCDM-моделью

[iLushkersky]

Доброго всем времени суток!

Решил построить график зависимости модуля расстояния и красного смещения с тремя космологическими моделями:

1. \Omega_M = 0.3, \Omega_\Lambda = 0.7
2. \Omega_M = 0.3, \Omega_\Lambda = 0
3. \Omega_M = 1.0, \Omega_\Lambda = 0

Значение постоянной Хаббла сост. 68 км/с/Мпк (также пробовал для 65 и 72, разницы особенной нет).

Вот сама функция (для \Omega_k = 0):

\mu_p=5\log\left(cH^{-1}_0(1+z)\int_{0}^{z}[(1+z)^2(1+\Omega_Mz)-z(2+z) \Omega_{\Lambda}]^{-1/2}\,dz\right)+25

В качестве подопытных использовал свечи из обзоров: SDSS-II, SNLS, DES - каково было мое удивление, когда практически все показатели соответствовали модели с плотностью материи, равной единице; решил взять сверхновые из оригинального исследования Supernova Cosmology Project и High-z Search Team - результат тот же (фотографии ниже):






Делал в MatLab (код):

c = 299792.458;
Omega_m1 = .3;
Omega_lam1 = .7;
Omega_m2 = .3;
Omega_lam2 = 0;
Omega_m3 = 1.0;
Omega_lam3 = 0;
H = 68;

syms z
f = 5 * log10(c * H^(-1) * (1+z) *...
    int((...
        (1+z)^2 * (1+Omega_m1*z) - z*(2+z)*Omega_lam1...
        )^(-1/2), [0, z]) )+ 25;

g = 5 * log10(c * H^(-1) * (1+z) *...
    int((...
        (1+z)^2 * (1+Omega_m2*z) - z*(2+z)*Omega_lam2...
        )^(-1/2), [0, z]) )+ 25;

u = 5 * log10(c * H^(-1) * (1+z) *...
    int((...
        (1+z)^2 * (1+Omega_m3*z) - z*(2+z)*Omega_lam3...
        )^(-1/2), [0, z]) )+ 25;

(далее fplot и вставка данных из xlsx-таблицы)

Сетую на какую-то ошибку в записи, хотя перепроверял - все ок, не заметил. Тут глаз много, может быть кто-нибудь что-нибудь найдет?

[Pphantom]

iLushkersky, а вы не могли бы объяснить цель этих упражнений в построении графиков? Я уж не говорю о том, что вообще-то это возобновление темы из Карантина.