(require (planet "geometric.ss" ("williams" "science.plt")
                 "random-distributions"))
(require (planet "discrete-histogram-with-graphics.ss"
                 ("williams" "science.plt")))

(let ((h (make-discrete-histogram)))
  (do ((i 0 (+ i 1)))
      ((= i 10000) (void))
    (discrete-histogram-increment! h (random-geometric 0.5)))
  (discrete-histogram-plot h "Histogram of Geometric Distribution"))