We study general gauge-dependent dynamical equations describing homogeneous isotropic cosmologies coupled to a scalar field ψ (scalaron). For flat cosmologies (k = 0), we analyze in detail the previously proposed gauge-independent equation describing the differential, χ(α) ≡ ψ ′ (α), of the map of the metric α to the scalaron field ψ, which is the main mathematical characteristic locally defining a 'portrait' of a cosmology in the so-called 'α-version'. In a more habitual 'ψ-version', the similar equation for the differential of the inverse map,χ(ψ) ≡ χ −1 (α), can be solved asymptotically or for some 'integrable' scalaron potentials v(ψ). In the flat case,χ(ψ) and χ(α) satisfy the first-order differential equations depending only on the logarithmic derivative of the potential, l(ψ) ≡ v ′ (ψ)/v(ψ). Once we know a general analytic solution for one of these χ-functions, we can explicitly derive all characteristics of the cosmological model.In the α-version, the whole dynamical system is integrable for k = 0 and with any 'α-potential',v(α) ≡ v[ψ(α)], replacing v(ψ). There is no a priori relation between the two potentials before deriving χ(α) orχ(ψ), which implicitly depend on the potential itself, but relations between the two pictures can be found by asymptotic expansions or by inflationary perturbation theory. We also consider alternative proposals -to specify a particular cosmology by guessing one of its portraits and then finding (reconstructing) the corresponding potential from the solutions of the dynamical equations.The main subject of this paper is the mathematical structure of isotropic cosmologies, but some explicit applications of the results to a more rigorous treatment of the chaotic inflation models and to their comparison with the ekpyrotic-bouncing ones are outlined in the frame of our 'α-formulation' of isotropic scalaron cosmologies. In particular, we establish an inflationary perturbation expansion for χ(α). When all the conditions for inflation are satisfied, which are: v > 0, k = 0, χ 2 (α) < 6, and χ(α) obeys a certain boundary (initial) condition at α → −∞, the expansion is invariant under scaling of v and its first terms give the standard inflationary parameters, with higher-order corrections. When v < 0 and 6χ 2 < 1 our general approach can be applied to studies of more complex ekpyrotic solutions alternative to inflationary ones. * Alexandre.Filippov@jinr.ru 1