“…Theoretical advances can be found in [1,4,5,7,9,13,23], including the derivation of expansions for solutions in the vicinity of a cosmological singularity. Of course, in most cases the solutions are not known in a close form (with the exception of, for instance, [6,[17][18][19]) and therefore we must appeal to numerical approximation in order to analyze the qualitative properties of solutions. Since the solutions of interest blow-up on the singularity t = 0, it it necessary to develop an adapted methodology in order to reach reliable conclusions on, for instance, the stability of solutions with respect to their initial data.…”