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]