Gnuplot scripts

hf.gnu

set key c
set grid y
set grid x
set xtics nomirror
set x2tics
set xlabel 'Cutoff Energy Ry"
set x2label 'Number of RL Vectors in Exchange' 
p "hf.dat" u 1:($4-$3)  w lp lw 2 lt -1 t "HF gap",  using 1:($4-$3):x2tic(2) lt 7 t