As in our prior work, we neglect reactions occurring in the vapor phase and the effects of the very slow film growth chemistry on the composition of the bulk liquid in each tray of the distillation column. We solve this tightly coupled vapor-liquid system of ordinary differential equations (ODEs) using a variant of Speth-Strang operator splitting iteratively solving,