The two equations, then, are
[tex]\epsilon'= \epsilon_\infty+ \frac{\epsilon_s- \epsilon_\infty}{1+ (\omega \tau)^2}[/tex]
and
[tex]\epsilon"= \frac{(\epsilon_s- \epsilon_\infty)\omega\tau}{1+ (\omega \tau)^2}[/tex]}
Of course, since there are only two equations in 3 unknown values, you cannot solve for all three. What you could do is solve for two of them in terms of the third.
You might, for example, multiply the first equation by [itex]\omega\tau[/itex], so that the two fractions are the same, and then subtract one equation from the other, eliminating the fractions:
[tex]\epsilon'- \epsilon"= \epsilon_\infty (\omega\tau)[/tex]
Then
[tex]\epsilon_\infty= \frac{\epsilon'-\epsilon"}{\omega\tau}[/tex]
so you have solved for [itex]\epsilon_\infty[/itex] in terms of [itex]\tau[/itex].
Replace [itex]\epsilon_\infty[/itex] by that in either of the two equations, and you can then also solve for [itex]\epsilon_\infty[/itex] in terms of [itex]\tau[/itex].