A new Monte Carlo method which permits the direct calculation of incoherent phase diagrams is presented. The method involves maintaining each lattice quite separately during the calculation and determining the phase boundaries by equating both the chemical potential and the pressure between phases. The results of several calculations are presented as examples. This new simulation procedure may be applied to any type of incoherent phase diagram without the complication of metastability.