np=40000 # number of points
n=100 #number of intervals
max=3. #max value
min=-3. #min value
width=(max-min)/n #interval width
#function used to map a value to the intervals
hist(x,width)=width*floor(x/width)+width/2.0
set boxwidth width*0.9

set term png
set o "gaussrng.png"

plot "res" u (hist($1,width)):(1.0) smooth freq w boxes lc rgb"green" notitle, width*np*exp(-x*x/2)/sqrt(2*pi)