The method described permits the computation of the concentrations of free ions and ion-ligand complexes in a solution containing arbitrary numbers of divalent cations and ligands. It is required that the pH be known, along with appropriate sets of ligand-hydrogen and ligand-divalent cation concentration binding constants. It is assumed that these sets of constants are chosen to be consistent with the ionic strength of the complete solution which contains the divalent cations and ligands. The technique is an iterative one which provides upper and lower bounds for the values of the unknowns. The method does not require initial guesses at the values of the unknowns, and it gives correct answers even when the concentrations involved are many orders of magnitude apart. The present formulation of the problem is restricted to the case where only one cation can bind to a given ligand at any one time. The method is applicable to large molecules with multiple "sub-ligands" provided these sub-ligands are independent in their function as ion-binding sites. These sub-ligands need not all have the same properties. It is also shown that a simple modification of the method permits the determination of the subset of total ion concentrations that are required in order to produce a specified subset of free ion concentrations. The modifications required to include monovalent cation binding are presented in outline form.