The resulting functions \(t(\mu )^{a_j}\)  parametrize the entropic path
for \( \mu \in (0, \epsilon ]\) .
The starting
point \(t(\epsilon )\)  is found by solving the
binomial equations \(x^*(\epsilon )_j = t(\epsilon )^{a_j} ,\, j = 1, \ldots n\) . This can be done by a Smith normal form computation. The tracking for \(\mu \rightarrow 0^+\)  is carried out with standard predictor-corrector techniques from numerical homotopy continuation [1]}.
