We have simulated the behaviour of a model, binary mixture of nematogens composed of cylindrically symmetric particles using the Monte Carlo technique. The characteristics of the model mixture were chosen to be in accord with most of the assumptions made in the Humphries–James–Luckhurst theory of liquid crystalline mixtures. The results of the simulation experiments allow us to test, for the first time, the validity of the molecular field approximation in this theory. In addition to the second rank long range orientational order parameters for both components of the mixture we have also determined certain orientational pair correlation functions. These enable us to investigate the ability of one component to enhance the order of the other component or the same component in its vicinity.