A simple mathematical model describing the separation of a series of aromatic bases by electrokinetic chromatography using beta-cyclodextrin (beta-CD) as a pseudostationary phase is described. The model takes into account changes in electrolyte pH and the different formation constants between the neutral and charged forms of the analytes with the CD. Constants in the model were obtained within the two-dimensional experimental space defined by pH and [beta-CD] with nonlinear regression using only five experimental points. These constants agreed with expected trends in analyte-CD interactions and predicted much higher formation constants for the neutral analyte-CD complex than for the charged analyte-CD complex. Correlation between predicted and observed mobilities using additional 20 points within the experimental space gave r(2) = 0.995. Optimisation of the pH and [beta-CD] was performed using both the normalised resolution product and minimum resolution product criteria and provided two optimum separations which exhibited different selectivities. Differences between predicted and observed migration times at these optima were less than 2.5 and 5% for the normalised resolution product and the minimum resolution criteria, respectively. In both cases the correct migration order was predicted. The model was also applied successfully to the optimisation of conditions for the separation of a specific mixture of analytes or for conditions under which particular analytes migrated in a desired order.