In classical analysis, one builds the catalog of special functions by repeatedly adjoining solutions of differential equations whose coefficients are previously known functions. Consequently, the properties of special functions depend crucially on the basic properties of ordinary differential equations. This naturally led to the study of formal differential equations, as in the seminal work of Turrittin [167]; this may be viewed retroactively as a theory of differential equations over a trivially valued field. After the introduction of p-adic analysis in the early 20th century, there began to be corresponding interest in solutions of p-adic differential equations; however, aside from some isolated instances (e.g., the proof of the Nagell-Lutz theorem; see Theorem 3.4), a unified theory of p-adic ordinary differential equations did not emerge until the pioneering work of Dwork on the relationship between p-adic special functions and the zeta functions of algebraic varieties over finite fields (e.g., see [57,58]). At that point, serious attention began to be devoted to a serious discrepancy between the p-adic and complex-analytic theories: on an open p-adic disc, a nonsingular differential equation can have a formal solution which does not converge in the entire disc (e.g., the exponential series). One is thus led to quantify the convergence of power series solutions of differential equations involving rational functions over a nonarchimedean field; this was originally done by Dwork in terms of the generic radius of convergence [59]. This and more refined invariants were studied by numerous authors during the half-century following Dwork's initial work, as documented in the author's book [93].At around the time that [93] was published, a new perspective was introduced by Baldassarri [13] (and partly anticipated in prior unpublished work of Baldassarri and Di Vizio [15]) which makes full use of Berkovich's theory of nonarchimedean analytic spaces. Given a differential equation as above, or more generally a connection on a curve over a nonarchimedean field, one can define an invariant called the convergence polygon; this is a function from the underlying Berkovich topological space of the curve into a space of Newton polygons, which measures the convergence of formal horizontal sections and is well-behaved with respect to both the topology and the piecewise linear structure on the Berkovich space. One can translate much of the prior theory of p-adic differential equations into (deceptively) simple statements about the behavior of the convergence polygon; this process was carried out in a series of papers by Poineau and Pulita [141,133,135,136], as supplemented by work of this author [98] and upcoming joint work with Baldassarri [16].