Efficient Computation of Higher-Order Variational Integrators in Robotic Simulation and Trajectory Optimization.