  We discuss a simulation algorithm for dynamical fermions, which combines the
multiboson technique with the Hybrid Monte Carlo algorithm. The algorithm turns
out to give a substantial gain over standard methods in practical simulations
and to be suitable for dealing with fermion zero modes in a clean and
controllable way.
