In [1]:

    

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

x=linspace(0,4,1000)
figure(figsize=(12,6))
fs=30
plot(x,1/(exp(10*(x-2))+1),lw=3)
grid()
ylim(0,1.2);
aprop=dict(color='k',width=.003)
yticks(linspace(0,1,3),(list(map(lambda x:r'$'+str(x)+r'$',linspace(0,1,3)))),fontsize=fs)
xticks(linspace(0,4,5),[r'$0$','',r'$\mu$',''],fontsize=fs)
xlabel(r'$E$',fontsize=fs)
ylabel(r'$f(E)$',fontsize=fs)
arrow(2,0.5,3/10,0,**aprop)
arrow(2,0.5,-3/10,0,**aprop)
text(1.8,0.6,r'$k_\mathrm{B}T$',fontsize=fs);


    

In [3]:

    

from scipy.special import *


    

In [4]:

    

x=linspace(-2,4,1000)
figure(figsize=(12,6))
fs=30
plot(x,zetac(x)+1,lw=3)
grid()
ylim(-3,3)
yticks(linspace(-3,3,7),(list(map(lambda x:r'$'+str(x)+r'$',linspace(-3,3,7)))),fontsize=fs)
xticks(linspace(-2,4,7),(list(map(lambda x:r'$'+str(x)+r'$',linspace(-2,4,7)))),fontsize=fs)
xlabel(r'$x$',fontsize=fs)
ylabel(r'$\zeta(x)$',fontsize=fs);