Nonholonomic Newmark method

Using the nonholonomic exponential map, we generalize the well-known family of Newmark methods for nonholonomic systems. We give numerical examples including a test problem where the structure of reversible integrability responsible for good energy behaviour as described in [16] is lost. We observe that the composition of two Newmark methods is able to produce good energy behaviour on this test problem.