A model for sampling inelastic hadron/nucleon interactions at energies below 5 GeV is described. The model describes the total hadronic cross section as a sum of quasi-two-body reaction channels. Inelastic events are sampled by first selecting one of these reaction channels and the subsequent decay of all resonances produced. Comparison with exclusive and inclusive data indicates very good agreement. The model is very well suited for hadron cascade calculations, detector simulation, and as input for models of hadron/nucleus interactions.