pro histo,array,x,h,m,s
; compute the histogram of an integer array
h=histogram(array)
x=lindgen(max(array)-min(array)+1)+min(array)
plot,x,h,psym=10
s=stdev(array,m)
return
end