  The derivation of second-order ordinary differential equations (ODEs) as
continuous-time limits of optimization algorithms has been shown to be an
effective tool for the analysis of these algorithms. Additionally, discretizing
generalizations of these ODEs can lead to new families of optimization methods.
We study discretizations of an Euler-Lagrange equation which generate a large
class of accelerated methods whose convergence rate is $O(\frac{1}{t^p})$ in
continuous-time, where parameter $p$ is the order of the optimization method.
Specifically, we address the question asking why a naive explicit-implicit
Euler discretization of this solution produces an unstable algorithm, even for
a strongly convex objective function. We prove that for a strongly convex
$L$-smooth quadratic objective function and step size $\delta<\frac{1}{L}$, the
naive discretization will exhibit stable behavior when the number of iterations
$k$ satisfies the inequality $k < (\frac{4}{Lp^2 \delta^p})^{\frac{1}{p-2}}$.
Additionally, we extend our analysis to the implicit and explicit Euler
discretization methods to determine end behavior.
