(require (planet "pareto.ss" ("williams" "science.plt" 2 0)
                 "random-distributions"))
(require (planet "histogram-with-graphics.ss"
                 ("williams" "science.plt" 2 0)))

(let ((h (make-histogram-with-ranges-uniform 40 1.0 21.0)))
  (do ((i 0 (+ i 1)))
      ((= i 10000) (void))
    (histogram-increment! h (random-pareto 1.0 1.0)))
  (histogram-plot h "Histogram of Pareto Distribution"))