Abstract
A general and numerical algorithm using Mathematica v. 2.0 is presented for calculating the stability ratios and interaction potential energy between two spherical colloidal particles according to the DLVO theory. The algorithm is applicable to electrolytes with any number of ionic species having any valence state and superficial potentials on the particles. The algorithm is versatile to allow change the values of the parameters that characterizes the colloidal system and is suitable to plot the stability curves for a given set of radii, Hamaker constant, superficial potential, electrolyte concentration range, and electrolyte valence. The algorithm was implemented on two models of DLVO interaction potential and is appropriate for analysis of coagulation tendencies of spherical colloidal particles.