  Traditional numerical discretizations of conservative systems generically
yield an artificial secular drift of any nonlinear invariants. In this work we
present an explicit nontraditional algorithm that exactly conserves these
invariants. We illustrate the general method by applying it to the three-wave
truncation of the Euler equations, the Lotka--Volterra predator--prey model,
and the Kepler problem. This method is discussed in the context of symplectic
(phase space conserving) integration methods as well as nonsymplectic
conservative methods. We comment on the application of our method to general
conservative systems.
