Plot
Needs["PlotLegends`"]
y[x_] := 1/(2 x - 1) Log[x/(1 - x)];
z[x_] := 1/(2 x) + 1/(2 (1 - x));
a = Plot[{y[x], z[x]}, {x, 0, 1}, AxesLabel -> {φ, χN},
PlotLegend -> {Style[
"coexistence", 13],
Style[
"spinodal", 13]}, LegendPosition -> {-0.72, -0.02},
LegendShadow -> False, LegendSize -> {1.1, 0.5}, LegendBorder -> None,
LegendBorderSpace -> Automatic, LegendTextSpace -> 7,
LegendBackground -> White, BaseStyle -> {FontSize -> 18},
PlotRange -> {0, 9}]
Export["Flory.eps", a]