bend.gp

set xlabel "q{/Symbol s}" offset 0,-0.75 font "Helvetica,30"
set ylabel "<|h(q{/Symbol s})|^{2}>/{/Symbol s}^{4}" offset -4.5,0 font "Helvetica,30"
#set xtics 0,30000,6000000
#set ytics 0,0.2,1.0
unset key
set logscale xy
set tics scale 3.0, 1.5
set tics font "Helvetica, 30"
plot "q_vs_ha.txt" lw 6 
k=4.79766
f(x) = 0.25 / (k*x**4)
replot f(x) lw 6 lc rgb "blue"
set terminal postscript eps enhanced color solid
set output "bend.eps"
replot