set xlabel "Energy (Ry)"
set ylabel "Bloch spectrum function (arb. unit)"

set xrange [-0.1:0.1]
#set xrange [0.04:0.06]
#set xrange [-0.04:-0.02]
set yrange [0:1000]
unset ytics

plot "GaAsGamma.txt" u 1:2 w p pt 1 notitle

f1(x)=a11/(1+a13*(x-a12)*(x-a12))
a11 = 500 ; a12= 0.05 ; a13 = 1000000
fit [0.04:0.06] f1(x) "GaAsGamma.txt" using 1:2 via a11,a12,a13
rep f1(x)

f2(x)=a21/(1+a23*(x-a22)*(x-a22))
a21 = 1000 ; a22= -0.03 ; a23 = 1000000
fit [-0.04:-0.02] f2(x) "GaAsGamma.txt" using 1:2 via a21,a22,a23
rep f2(x)

set terminal pngcairo size 520,390
set output "GaAsGamma.png"
rep