We study the acceptance probability in the hybrid Monte Carlo method applied to lattice QCD with dynamical fermions We investigate in detail its dependence upon the step-size and the parameters [beta], m q and V, and suggest a simple phenomenological formula for this dependence