from matplotlib.pyplot import *
from numpy import *
from mpl_toolkits.axisartist import *

ni = lambda x: 1. / ( exp(1./x) + 1. )
x = linspace(0,15,100)

fig = figure(figsize=(4,4))
ax = Subplot(fig,"111",axisbg='none')
fig.add_subplot(ax)
ax.axis["right"].set_visible(False)
ax.axis["top"].set_visible(False)
ax.plot(x,ni(x),'k',lw=2)
ax.plot([0,15],[.5,.5],'k--',lw=2.)
ax.set_ylabel(r'$\bar{n}$',fontsize=17)
ax.set_xlabel(r'$\frac{kT}{\epsilon-\mu}$',fontsize=17)
ax.set_xticks(arange(0,16,3))
ax.grid(True)
ax.minorticks_on()
ax.set_xlim(0,15)
ax.set_ylim(0,.6)
fig.savefig("FD_kT_e.svg",bbox_inches="tight",\
	pad_inches=.15)