Discretization scheme in Volterra integro-differential equations that preserves stability and boundedness

A nonstandard discretization scheme is applied to continuous Volterra integro-differential equations. We will show that under our discretization scheme the stability of the zero solution of the continuous dynamical system is preserved. Also, under the same discretization, using a combination of Lyapunov functionals, Laplace transforms and z-transforms, we show that the boundedness of solutions of the continuous dynamical system is preserved.