  We propose a numerical scheme for the time-integration of nonholonomic
mechanical systems, both conservative and nonconservative. The scheme is
obtained by simultaneously discretizing the constraint equations and the
Herglotz variational principle. We validate the method using numerical
simulations and contrast them against the results of standard methods from the
literature.
