Lyapunov stability of the Euler method

We extend the Lyapunov stability criterion to Euler discretizations of differential inclusions. It relies on a pair of Lyapunov functions, one in continuous time and one in discrete time. In the context of optimization, this yields sufficient conditions for the stability of nonisolated local minima when using the Bouligand subgradient method.