High order symmetric algorithms for nonlinear dynamical systems with non-holonomic constraints