c=.001
a=10.
b=10000. #-10000000.
d=-0.02

kB=8.9*10.**(-5.)  #eV/K


#p "Sun_FeAs2_B.dat" u 1:(-$2) w l lw 4 ti "Sun et al : FeAs2" 
#rep "0.2/d.0005/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lw 2
#rep "0.2/d.0005/response.dat_1" u 1:((-c*$1/(b+a*$1)-kB*10.**3.*$4) / (1./(b+a*$1)+$3)) w l lw 2


set log x
set log y

set xrange[1:1000]
set yrange[0.001:1000000.]

p "0.2/d.0005/response.dat_1" u 1:(10.**(-3.)/$3) w l lw 2

rep "0.2/d.0005/response.dat_1" u 1:(1./( 10.**(3.)*$3 + 1./(b*exp(d*$1)) ) ) w l lw 2