All pastes #676108 Raw Edit

Someone

public text v1 · immutable
#676108 ·published 2007-08-30 19:55 UTC
rendered paste body
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]