An addendum regarding integrators - notice all of the expressions, right before the call to dSolveLCP, containing "stepsize1" (in dInternalStepIsland_x2). These appear to be the dependencies that the time-stepping LCP formulation has on the (currently Euler) integration method; probably, it is these which would need to be changed to reflect a non-Euler integrator. -Norman