A method for computing a solution speciation in a closed system or at equilibrium with respect to mineral or gas phases is presented. Each species is expressed as a function of a basis of components. An optimum choice of these components in relation to the prescribed conditions makes it possible on the one hand, to clearly pose the unknowns and the equations, and on the other hand, to reduce computer time and numerical problems.