set xlabel 'E (MeV)' set ylabel 'Events' set log y set style fill solid filter(x,y)=floor(x/y)*y plot 'fort.93' u (filter($1,0.01)):(1) smooth frequency with boxes lt 3 title ''