  We present a simulation algorithm for dynamical fermions that combines the
multiboson technique with the Hybrid Monte Carlo algorithm. We find that the
algorithm gives a substantial gain over the standard methods in practical
simulations. We point out the ability of the algorithm to treat fermion
zeromodes in a clean and controllable manner.
