In this paper, a new computational method for improving the accuracy of numerically computed solutions is introduced. The computational method is based on the one-step method and conserved quantities of holonomic systems are considered as kinematical constraints in this method.