You can drop the {z,M} points into a spread sheet like Excel and take Dm = (10^(M/5))/(1+z) to get {z,Dm} pairs, which are well-behaved (Dm approx az - bz²). You'll need some integrating software to compare the LCDM model which predicts that "a" is proportional to sinh(qt) where q is a constant that depends on the parameters. Then Dm is proportional to the time integral of (1/a)dt. The limits are from the time tz where a(tz) = 1/(1+z) to the present time t0 (which is fixed).