function for_stats, vals, binsize=binsize, minval=minval, maxval=maxval

if n_elements(binsize) eq 0 then binsize = .01
if n_elements(minval) eq 0 then minval = 0.0
if n_elements(maxval) eq 0 then maxval = 1.0
nbins = (maxval-minval)/binsize + 1
half = binsize/2.0
res = fltarr(3,nbins)

j = 0
for i=minval,maxval-binsize,binsize do begin
	w = where( vals ge i and vals lt i+binsize, c )
	if c lt 3 then continue
	res[*,j] = [i+half,mean(vals[w]),stddev(vals[w])]
	j++
endfor

return,res

end