;+
; NAME:
; bluenoise
;
;
; PURPOSE:
; return a blue (f) noise spectra
;
;
; CATEGORY:
; stats, noise
;
;
; OPTIONAL INPUTS:
; n: number of samples to generate, default=1000
;
;
; KEYWORD PARAMETERS:
; RANGE: range over which to generate samples, default=[5,100]
; SEED: random number seed
;
;
; OUTPUTS:
; array of size N distributed by a blue-noise spectra (f power law)
;
;
; EXAMPLE:
; histoplot, /xlog, /ylog, /fill, yrange=[1,100000], bluenoise(10000, range=[5,100]) 
;
;
; MODIFICATION HISTORY:
;
;       Tue Mar 24 12:21:57 2009, Brian Larsen
;		written and tested
;
;-
FUNCTION bluenoise, n, $
                   RANGE = range, $
                   SEED = seed

  SetDefaultValue, n, 1000
  SetDefaultValue, range, [5, 100]
  randomp, ans, 2, n, range_x = range, SEED = seed


  RETURN, ans
END