To implement many-qubit gates for use in quantum simulations on quantum computers efficiently, we develop and present methods reexpressing exp[-i(H-1+H-2+...)Delta t] as a product of factors exp[-iH(1)Delta t], exp[-iH(2)Delta t],..., which is accurate to third or fourth order in Delta t. The methods we derive are an extended form of the symplectic method, and can also be used for an integration of classical Hamiltonians on classical computers. We derive both integral and irrational methods, and find the most efficient methods in both cases. [S1050-2947(99)07209-1].