Numerical calculation of strongly coupled plasma composition

The ionization balance equations are considered for multicomponent plasma of arbitrary complex structure. Partial degeneracy of electrons and ionization potential shift, caused by non-ideality of plasma (interaction of charged particles), are taken into account. Quick enough and robust algorithm for solution of these equations is developed. The example with Li on conditions that plasma density verge towards condenced matter and non-ideality parameter reaches 6-7 illustrates efficiency of the method.