set samples 2000
set dummy EN
T = 10000.
e = 1.602e-019
k = 1.380662E-23
k = k / e
#
# Fermi-Energie von Cu
#
E_F = 7.04
#
f_Boltzmann(EN) = 1./(k*T) * exp(-EN/(k*T))
f_Fermi(EN)  = 1./( E_F+k*T*log( 1.+exp( -E_F/(k*T) ) ) ) * 1./( exp( (EN-E_F) / (k*T) ) + 1.)
#
# willk. Normierungskonst. für BE
#
const = 0.5
f_Bose(EN)   = const / ( exp( EN/(k*T) ) - 1. )
#
set xrange [0:12.]
set xlabel "Energie E /eV"
set yrange [0:1.2]
set ytics 0,0.1
set ylabel "Verteilungsfkt. f(E)"
set grid
set title "Boltzmann-, Bose-Einstein- und Fermi-Dirac-Verteilung"
set label 1 " T = 10000 K , Kupfer" at 7.,0.25
plot f_Boltzmann(EN), f_Fermi(EN), f_Bose(EN)
#    EOF