(require (planet "bivariate-gaussian.ss" ("williams" "science.plt")
                 "random-distributions"))
(require (planet "histogram-2d-with-graphics.ss"
                 ("williams" "science.plt")))

(let ((h (make-histogram-2d-with-ranges-uniform
          20 20 -3.0 3.0 -3.0 3.0)))
  (do ((i 0 (+ i 1)))
      ((= i 10000) (void))
    (let-values (((x y) (random-bivariate-gaussian 1.0 1.0 0.0)))
      (histogram-2d-increment! h x y)))
  (histogram-2d-plot
   h "Histogram of Bivariate Gaussian Distribution"))