The release of two catalogs for the distance modulus of Supernova (SN) of type Ia, namely, the Union 2.1 compilation, see  , and the joint light-curve analysis (JLA), see  , allows matching the observed distance modulus with the theoretical distance modulus of various cosmologies. In this fitting procedure, the cosmological parameters are derived in a scientific and reproducible way.
We now focus our attention on the flat Friedmann-Lemaître-Robertson-Walker (flat-FLRW) cosmology. A first fitting formula has been derived by  and an approximate solution in terms of Padé approximant has been introduced by  . The presence of the elliptical integrals of the first kind in the integral for the luminosity distance in flat-FLRW cosmology has been noted by    . As a practical example the luminosity distance can be expanded into a series of orthonormal functions and the two cosmological parameters turn out to be and , see  . This paper first introduces in Section 2 a framework useful to build a new solution for the luminosity distance in flat-FLRW cosmology, which will be derived in Section 3.
This section reviews the adopted statistical framework, the ΛCDM cosmology, and an existing solution for the luminosity distance in flat-FLRW cosmology.
2.1. The Adopted Statistics
In the case of the distance modulus, the merit function is
where N is the number of SNs, is the observed distance modulus evaluated at redshift , is the error in the observed distance modulus evaluated at , and is the theoretical distance modulus evaluated at , see formula (15.5.5) in  . The reduced merit function is
where is the number of degrees of freedom, N is the number of SNs, and k is the number of parameters. Another useful statistical parameter is the associated Q-value, which has to be understood as the maximum probability of obtaining a better fitting, see formula (15.2.12) in  :
where GAMMQ is a subroutine for the incomplete gamma function.
The goodness of the approximation in evaluating a physical variable p is evaluated by the percentage error
where is an approximation of p.
2.2. The Standard Cosmology
We follow  , where the Hubble distance is defined as
The first parameter is
where G is the Newtonian gravitational constant and is the mass density at the present time. The second parameter is
where is the cosmological constant, see  . These two parameters are connected with the curvature by
The comoving distance, , is
where is the “Hubble function”
The above integral does not have an analytical solution but a solution in terms of Padé approximant has been found, see  .
2.3. A First Formula for a Flat-FLRW Universe
The first model starts from Equation (2.1) in  for the luminosity distance, ,
where is the Hubble constant expressed in km・s−1・Mpc−1, c is the speed of light expressed in km・s−1, z is the redshift and a is the scale-factor. The indefinite integral, , is
The solution is in terms of F, the Legendre integral or incomplete elliptic integral of the first kind, and is given in  .
The luminosity distance is
where means the real part. The distance modulus is
3. A New Formula for a Flat-FLRW Universe
The second model for the flat cosmology starts from Equation (1) for the luminosity distance in 
The above formula can be obtained from formula (9) for the comoving distance inserting and the variable of integration, t, denotes the redshift.
A first change in the parameter introduces
and the luminosity distance becomes
The following change of variable, , is performed for the luminosity distance which becomes
Up to now we have followed  which continues introducing a new function ; conversely we work directly on the resulting integral for the luminosity distance: which is
where s is given by Equation (16) and is Legendre’s incomplete elliptic integral of the first kind,
see  . The distance modulus is
with s as defined by Equation (16).
In recent years, the extraction of the cosmological parameters from the distance modulus of SNs has become a common practice, see among others    . The best fit to the distance modulus of SNs is here obtained by implementing the Levenberg-Marquardt method (subroutine MRQMIN in  ). This method requires the fitting function, in our case Equation (22), as well the first derivative
, which has a simple expression, and the first derivative ,
which has a complicated expression. A simplification can be introduced by imposing a fiducial value for the Hubble constant, namely , see   . We call this model “flat-FLRW-1”, where the “1” stands for there being one parameter. Table 1 presents and for the Union 2.1 compilation of SNs and Figure 1 displays the best fit. The reading of this table allows to evaluate the goodness of the approximation, see (4), in the derivation of the Hubble constant in going from the supposed true value ( ) to the deduced value ( ), which is . The
Table 1. Numerical values from the Union 2.1 compilation of , and Q, where k stands for the number of parameters.
Figure 1. Hubble diagram for the Union 2.1 compilation. The solid line represents the best fit for the exact distance modulus in flat-FLRW cosmology as represented by Equation (22). Parameters as in first line of Table 1; Union 2.1 compilation.
JLA compilation is available at the Strasbourg Astronomical Data Centre (CDS), and consists of 740 type I-a SNs for which we have the heliocentric redshift, z, the
apparent magnitude in the B band, the error in , , the parameter X1,
the error in X1, , the parameter C, the error in C, , and . The observed distance modulus is defined by Equation (4) in  :
The adopted parameters are , and
see line 1 in Table 10 of  . The uncertainty in the observed distance modulus, , is found by implementing the error propagation equation (often called the law of errors of Gauss) when the covariant terms are neglected, see Equation (3.14) in  ,
Table 2. Numerical values from the JLA compilation of , and Q, where k stands for the number of parameters.
Figure 2. Hubble diagram for the JLA compilation. The solid line represents the best fit for the exact distance modulus in flat-FLRW cosmology as represented by Equation (22). Parameters as in first line of Table 2.
As an example the luminosity distance for the Union 2.1 compilation with data as in the first line of Table 1 is
and the distance modulus is
We now derive some approximate results without Legendre integral for the flat-FLRW case and Union 2.1 compilation with data as in Table 1, first line. A Taylor expansion of order 6 around z = 0 of the luminosity distance as given by Equation (19) for the flat-FLRW case and Union 2.1 compilation gives
The upper limit in redshift, 0.197, is the value for which the percentage error, see Equation (4), is . The asymptotic expansion of the luminosity distance with respect to the variable z to order 5 for the flat-FLRW case and Union 2.1 compilation gives
At the lower limit of the percentage error is . The two above approximations at low and high redshift have a limited range of existence but does not contain the Legendre integral as solutions (28) and (29) which cover the overall range .
A Taylor expansion of order 6 of the distance modulus as given by Equation (22) around for the flat-FLRW case and Union 2.1 compilation gives
The upper limit in redshift, 0.197, is the value at which the percentage error is . Figure 3 reports both the numerical and the Taylor expansion of distance modulus in the above range.
Figure 3. Distance modulus in flat-FLRW cosmology as represented by Equation (22) with parameters as in first line of Table 1 (full red line) and Taylor solution (dash-dot-dash line) (blue line).
The asymptotic expansion of the distance modulus with respect to the variable z to order 5 for the flat-FLRW case and Union 2.1 compilation gives
The lower limit in redshift, 1.27, is the value at which the percentage error is . The ranges of existence in z for the analytical approximations here derived have the percentage error <2%, see Equation (4).
We now introduce the best minimax rational approximation, see    , of degree (2, 1), for the distance modulus ,
In the case in which the distance modulus is represented by Equation (29) and given the interval , the coefficients of the best minimax rational approximation are presented in Table 3; the maximum error for the fit is . Figure 4 displays the data and the fit.
We have presented an analytical approximation for the luminosity distance in terms of elliptical integrals with real argument. The fit of the distance modulus
Table 3. Maximum error and coefficients of the distance modulus for the best minimax rational approximation for the flat-FLRW case and Union 2.1 compilation. Interval of existence .
Figure 4. Distance modulus in flat-FLRW cosmology as represented by Equation (22) with parameters as in first line of Table 1 (full line) and minimax rational approximation (empty stars); Union 2.1 compilation.
of SNs of type Ia allows finding the parameters and for the two compilations in flat-FLRW cosmology
A first comparison with  in the case of the Union 2.1 compilation gives a percentage error for the derivation of and for the derivation of . A second comparison can be done with Equation (13) in 
Figure 5. Distance modulus for ΛCDM cosmology (full line), flat-FLRW-1 (dot-dash-dot- dash line) and flat-FLRW cosmology (dashed line). Parameters as in Table 1 and interval of existence .
In the case of the Union 2.1 compilation, the percentage error for the derivation of and for . A Taylor expansion at low redshift and an asymptotic expansion are presented both for the luminosity distance and the distance modulus. A simple version of the distance modulus is determined through the best minimax rational approximation. Adopting the cosmological parameters found here, the cosmological constant turns out to be, for the Union 2.1 compilation,
flat-FLRW Union 2.1,
or introducing and the Planck time, ,
flat-FLRW Union 2.1.
The statistical parameters of the fits are given in Table 1 and Table 2 where the other two models are presented. The values of the in the above table say that for the Union 2.1 compilation the flat cosmology produces a better fit than the ΛCDM does, but the situation is the reverse for the JLA compilation. As a concluding remark we point out that, thanks to the calibration on the distance modulus of SNs, the differences between the solutions here analyzed are minimum. Therefore a restricted range in redshift should be adopted in order to visualize the diverseness, see Figure 5.
 Suzuki, N., Rubin, D., Lidman, C., Aldering, G., Amanullah, R., Barbary, K. and Barrientos, L.F. (2012) The Hubble Space Telescope Cluster Supernova Survey. V. Improving the Dark-Energy Constraints above Z Greater than 1 and Building an Early-Type-Hosted Supernova Sample. ApJ, 746, 85.
 Adachi, M. and Kasai, M. (2012) An Analytical Approximation of the Luminosity Distance in Flat Cosmologies with a Cosmological Constant. Progress of Theoretical Physics, 127, 145.
 Benitez-Herrera, S., Ishida, E.E.O., Maturi, M., Hillebrandt, W., Bartelmann, M. and Röpke, F. (2013) Cosmological Parameter Estimation from SN Ia Data: A Model-Independent Approach. MNRAS, 436, 854.
 Zaninetti, L. (2016) An Analytical Solution in the Complex Plane for the Luminosity Distance in Flat Cosmology. Journal of High Energy Physics, Gravitation and Cosmology, 2, 581.
 Jones, D.O., Scolnic, D.M., Riess, A.G., Rest, A., Kirshner, R.P., Berger, E., Kessler, R., Pan, Y.C., Foley, R.J., Chornock, R., Ortega, C.A., Challis, P.J., Burgett, W.S., Chambers, K.C., Draper, P.W., Flewelling, H., Huber, M.E., Kaiser, N., Kudritzki, R.P., Metcalfe, N., Tonry, J., Wainscoat, R.J., Waters, C., Gall, E.E.E., Kotak, R., McCrum, M., Smartt, S.J. and Smith, K.W. (2018) Measuring Dark Energy Properties with Photometrically Classified Pan-STARRS Supernovae. II. Cosmological Parameters. The Astrophysical Journal, 857, 51.
 Ben-Dayan, I., Gasperini, M., Marozzi, G., Nugier, F. and Veneziano, G. (2013) Average and Dispersion of the Luminosity-Redshift Relation in the Concordance Model. Journal of Cosmology and Astroparticle Physics, 6, 2.