kBJ=1.38*10.**(-23.) #J/K
kB=8.9*10.**(-5.)  #eV/K
e=1.602*10.**(-19.)
r=-0.5
m=9.81*10.**(-31.)   #kg
#h=1.   #Js this is hbar.
hbar=(6.62607*10.**(-34.))#/(2.0*3.141592653589793238)  #hbar in Js

set xlabel "T [eV]"
set ylabel "n [1/m^3]

# S=pi^2 / 3 kB/e kBT/eF (r+3/2)

# eF=h^2/(2m) (3n/(8pi))^(2/3)
#or
# eF=pi^2/3 kB/e kBT/S (r+3/2)

# i.e.
# n = (2me_F/h^2)^3/2 *8pi/3




#n=
#( 2.*m/hbar^2*((pi**(2.))/3.*kB/e *kB*$1/$2 * (r+3./2.)))**(3./2.) *8.*pi/3.

#    S was in mV/K., #n in 1/cm^3
units=1000. # / (100.)**3.

p "Sun_FeSb2_model.dat" u 1:( (2.*m/(hbar**2.)*(pi**2.)/3.*kBJ/e *kBJ*$1/$2 * (r+3./2.))**(3./2.) *8.*pi/3.) *units w l lw 2 ti ""

# Js*Js/Kg  = J m2   S in mV/K
#J=kg m2/s2 , n=1/m3