  A geometric derivation of numerical integrators for optimal control problems
is proposed. It is based in the classical technique of generating functions
adapted to the special features of optimal control problems.
