One of the continuing challenges in cosmology has been to determine the large-scale space-time metric from observations with a minimum of assumptions -- without, for instance, assuming that the universe is almost Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW). If we are lucky enough this would be a way of demonstrating that our universe is FLRW, instead of presupposing it or simply showing that the observations are consistent with FLRW. Showing how to do this within the more general spherically symmetric, inhomogeneous space-time framework takes us a long way towards fulfilling this goal. In recent work researchers have shown how this can be done both in the traditional Lema\^{i}tre-Tolman-Bondi (LTB) 3 + 1 coordinate framework, and in the observational coordinate (OC) framework. In this paper we investigate the stability of solutions, and the use of data in the OC field equations including their time evolution and compare both approaches with respect to the singularity problem at the maximum of the angular-diameter distance, the stability of solutions, and the use of data in the field equations. This allows a more detailed account and assessment of the OC integration procedure, and enables a comparison of the relative advantages of the two equivalent solution frameworks. Both formulations and integration procedures should, in principle, lead to the same results. However, as we show in this paper, the OC procedure manifests certain advantages, particularly in the avoidance of coordinate singularities at the maximum of the angular-diameter distance, and in the stability of the solutions obtained. This particular feature is what allows us to do the best fitting of the data to smooth data functions and the possibility of constructing analytic solutions to the field equations.

A rigorous demonstration that given appropriate data on our past light cone leads to the determination of the metric functions and all their time derivatives on our past light cone is presented, thus showing how to evolve the solution we obtain from data on the light cone off it in a well-defined and straightforward way. It also automatically gives a procedure for constructing the solution for all spherically symmetric, inhomogeneous cosmological Lema\^itre-Tolman-Bondi models in observational coordinates as a Taylor series in time of however many terms we need. Our procedure takes into account the essential data giving the maximum of the observer area (angular-diameter) distance, and the redshift $z_{max}$ at which that occurs. This enables the determination of the vacuum-energy density $\mu_{\Lambda}$, which would otherwise remain undetermined.

Extending the approach developed by Ara\'ujo and Stoeger [1] and improved in Ara\'ujo {\it et al} [2], we have shown how to construct dust-filled $\Lambda \neq 0$ Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) cosmological models from FLRW cosmological data on our past light cone. Apart from being of interest in its own right -- demonstrating how such data fully determines the models -- it is also illustrated in the flat case how the more general spherically symmetric (SS) Einstein field equations can be integrated in observational coordinates with data fit to FLRW forms arrayed on our past light cone, thus showing how such data determines a FLRW universe -- which is not {\it a priori} obvious. It is also shown how to integrate these exact SS equations, in cases where the data are not FLRW, and the space-time is not known to be flat. It is essential for both flat and non-flat cases to have data giving the maximum of the observer area (angular-diameter) distance, and the redshift $z_{max}$ at which that occurs. This enables the determination of the vacuum-energy density $\mu_{\Lambda}$, which would otherwise remain undetermined.

Assuming that the cosmological principle holds, Maartens, Ellis and Stoeger (MES) recently constructed a detailed scheme linking anisotropies in the cosmic background radiation (CMB) with anisotropies and inhomogeneities in the large scale structure of the universe and showed how to place limits on those anisotropies and inhomogeneities simply by using CMB quadrupole and octupole limits. First we indicate and discuss the connection between the covariant multipole moments of the temperature anisotropy used in the MES scheme and the quadrupole and octupole results from COBE. Then we introduce those results into the MES limit equations to obtain definite quantitative limits on the complete set of cosmological measures of anisotropy and inhomogeneity. We find that all the anisotropy measures are less than 10^{-4} in the case of those not affected by the expansion rate H, and less than 10^{-6} Mpc^{-1} in the case of those which are. These results quantitatively demonstrate that the observable universe is indeed close to Friedmann-Lemaitre-Robertson-Walker (FLRW) on the largest scales, and can be adequately modelled by an almost-FLRW model -- that is, the anisotropies and inhomogeneities characterizing the observable universe on the largest scales are not too large to be considered perturbations to FLRW.