Aqueous and non-aqueous electrolyte solutions are ubiquitous in chemical and biochemical applications, especially in innovative processes, and they play a major role in geochemistry, environmental science, and numerous other scientific fields. Despite the obvious importance of electrolyte systems, research success on electrolyte thermodynamics is still behind all of the advances on non-electrolyte thermodynamics. After decades of research, several issues of thermodynamic models for electrolytes remain the object of discussion in the thermodynamic community. Still today, only a few simulation packages offer a general approach to calculate phase equilibria of electrolyte systems in a broad application range regarding the kind of salts and solvent, the number of phases, and the number of non-ionic or ionic species involved. In this work, the general background and the assumptions behind the equilibrium conditions of multiphase electrolyte systems with distributed ions are reviewed. A general methodology is proposed, which can be used to determine the number of liquid phases and their composition at equilibrium of any electrolyte system independent of the number of components. The algorithm was implemented in a FORTRAN routine using the equation of state “ePC-SAFT advanced” (Bülow, M.; Ascani, M.; Held, C. ePC-SAFT advanced-Part I: Physical meaning of including a concentration-dependent dielectric constant in the born term and in the Debye–Hückel theory. Fluid Phase Equilib. 2021, 535, 112967) to estimate the fugacity of each species. The algorithm was successfully tested against experimental data using case studies including three-phase liquid–liquid–liquid equilibria with two ionic species and two-phase liquid–liquid equilibria with three ionic species.