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

kT = lambda t: 1./t
ni = lambda x, t: 1. / ( exp((x-1.)/kT(t)) + 1. )
x = linspace(0,5,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,100),'c',label=r'$kT=\mu/100$',lw=2)
ax.plot(x,ni(x,10),'b',label=r'$kT=\mu/10$',lw=2)
ax.plot(x,ni(x,2),'r',label=r'$kT=\mu/2$',lw=2)
ax.plot(x,ni(x,1),'y',label=r'$kT=\mu$',lw=2)

ax.set_ylabel(r'$\bar{n}$',fontsize=15)
ax.set_xlabel(r'$\epsilon/\mu$',fontsize=12)
ax.set_yticks(arange(0,1.2,.2))
ax.grid(True)
ax.minorticks_on()
ax.set_xlim(0,5)
ax.set_ylim(-.1,1.1)
ax.legend(frameon=False, handletextpad=0,loc='upper right',\
	labelspacing=.05)
setp(ax.get_legend().get_texts(),fontsize=12)
fig.savefig("FD_e_mu.svg",bbox_inches="tight",\
	pad_inches=.15)