Gibbs ensemble Monte Carlo simulations have been used to calculate vapour-liquid equilibria of diatomic and dipolar diatomic fluids. Dipolar two-centre Lennard-Jones diatomics with elongations L* = L/sigma between 0.33 and 1.00 have been studied for various values of dipole moments mu(d)* = mu(d)/(epsilonsigma3)1/2. The results for phase equilibrium are in good agreement with recent theoretical predictions based on cluster-expansion perturbation theory.