Numerical verification of the order of the asymptotic solutions of a nonlinear differential equation (Q2455816)

Consider the initial value problem \[ \ddot x+ u-{\varepsilon\over c} u^2= c,\quad u(0)= a,\quad \dot u(0)= 0,\tag{\(*\)} \] where \(\varepsilon\) is a small parameter, \(c\) and \(a\) are given constants. The author applies the Lindstedt-Poincaré method to get the asymptotic expansion of the solution \(u(t,\varepsilon)\) of \((*)\) in the form \[ u(t,\varepsilon)= \sum_{n= 0} \varepsilon^n u_n(\tau),\quad \tau=\omega t, \] \[ \omega(\varepsilon)= \sum_{n=0} \varepsilon^n \omega_n. \] Finally, he presents a numerical procedure to verify the order of an asymptotic expansion.