Log10[x/(1 - x)] == A + y
z == B w
z + v == C y
Ee == w + y
v == D x

eqns = {x4, x5, x6, x8, x9}

vars = {x, y, z, w, v};

(* Makes a table of solutions for values of Ee that vary from 0.3 to 0.7 *)
sols[evalue_] := NSolve[eqns /. Ee -> evalue, vars, Reals]
sollist = Table[sols[Ee], {Ee, 0.3, 0.7, 0.05}]

(* Compiles a list of log10[x/(1-x)] *)
masterlist = {}
Do[logx1x[sollist[[i, 1]]], {i, 9}]

logx1x = Function[singlesol,
frac = singlesol[[1, 2]];
frac = N[Log10[frac/(1 - frac)]];
Print[frac];
AppendTo[masterlist, frac]
]