We investigate the algorithmic approximation of ordinary differential equations having a known conservation law, with finite difference schemes which inherit a discrete version of the conservation law. We use the method of moving frames on a multispace due to Olver. We assume that the system of ODEs to be studied has a variational principle and that the conservation law arises from a variational symmetry via Noether's theorem.