The mathematical model outlined in Part I is recast in a form suitable for numerical computation. The spatial derivatives are replaced by finite-difference expressions, which leads to a set of ordinary differential equations coupled to a set of nonlinear algebraic relations. This system is solved using existing integration techniques. The resulting algorithm simulates the characteristic behavior of the classical modes of electrophoresis, which is shown by examples involving moving boundary electrophoresis and isoelectric focusing. In the first example two different integration schemes are used and their accuracy and stability investigated. The second example illustrates the versatility of the methodology.