I am trying to find the best way to draw from a normal distribution. I want to be able to use the normal probability density function (and its cumulative) in Haskell.

To put it simply, I want to use the functions on this page without using the GSL binding... I am struggling to find the right package. Which one would you say is the best?

Thank you and have a great day!

# Best How To :

Here's an example which uses random-fu:

```
import Data.Random -- for randomness
import Text.Printf -- for printf
import Data.Foldable -- for the for_ loop
-- pdf and cdf are basically “Distribution -> Double -> Double”
main = do
-- defining normal distribution with mean = 10 and variation = 2
let normal = Normal (10 :: Double) 2
-- CDF
for_ [0..10] $ \i ->
printf "cdf(%2d): %.4f\n" i (cdf normal (fromInteger i))
-- PDF
putStrLn "---"
for_ [0..10] $ \i ->
printf "pdf(%2d): %.4f\n" i (pdf normal (fromInteger i))
```

Run it and you'll see the following output:

```
cdf( 0): 0.0000
cdf( 1): 0.0000
cdf( 2): 0.0000
cdf( 3): 0.0002
cdf( 4): 0.0013
cdf( 5): 0.0062
cdf( 6): 0.0228
cdf( 7): 0.0668
cdf( 8): 0.1587
cdf( 9): 0.3085
cdf(10): 0.5000
---
pdf( 0): 0.0000
pdf( 1): 0.0000
pdf( 2): 0.0001
pdf( 3): 0.0004
pdf( 4): 0.0022
pdf( 5): 0.0088
pdf( 6): 0.0270
pdf( 7): 0.0648
pdf( 8): 0.1210
pdf( 9): 0.1760
pdf(10): 0.1995
```

There are 2 important points here:

`Normal`

is a constructor defining a distribution. There are other distributions, such as `Uniform`

, etc. Distributions have different types, but they all belong to the `Distribution`

class.

`pdf`

and `cdf`

are class methods, which can operate on many distributions (perhaps not all, but many). The 1st parameter is the distribution and the 2nd – point at which PDF/CDF should be evaluated.