counts=intarr(n_intervals)

max_dist=max(distances)
interval=ceil(max_dist)/float(n_intervals)
radii=interval*(indgen(n_intervals)+1)
	FOR j=0, n_intervals-1 DO BEGIN
		counts[j]=n_elements(where(distances LE radii[j]))						
	ENDFOR

a=linfit(alog10(radii), alog10(counts), /double)
	dims=a[1]