  We report on large-scale Monte Carlo simulations of the three-dimensional
4-state Potts ferromagnet subject to quenched, random bond dilution. For small
dilutions the rather strong first-order phase transition of the pure system is
found to persist, whereas for larger dilutions the theoretically expected
softening to a continuous transition is confirmed. The properties of the
underlying disorder distributions of thermal observables are discussed and
illustrated with a few selected examples.
