      (mvnormal mean variance))))
 (defn normal
   [mean variance]
-  (new NormalDistribution mean (Math/sqrt variance)))
+  (new NormalDistribution mean (sqrt variance)))
 (defn pareto
   [scale shape]
   (new ParetoDistribution scale shape))

        p-double (double-array p)]
    (new EnumeratedIntegerDistribution x-int p-double)))
(defn exponential
  [mean]
  (new ExponentialDistribution mean))
(defn f-distribution
  [df1 df2]
  (new FDistribution df1 df2))
(defn gamma
  [shape rate]
  (new GammaDistribution shape (/ 1 rate)))
(defn geometric
  [p]
  (new GeometricDistribution p))

    ([d x] (/ (pdf (new TDistribution (:df d)) (/ (- x (:location d)) (:scale d))) (:scale d)))
    ([d] (fn [x] (pdf d x))))
  (log-pdf
    ([d x] (- (log-pdf (new TDistribution (:df d)) (/ (- x (:location d)) (:scale d))) (log (:scale d)) ))
    ([d] (fn [x] (log-pdf d x))))
  MultivariateNormalDistribution
  (log-pdf
