#!/usr/bin/env python

#-------------
# Load modules
#-------------

import random, scipy, numpy

print random.uniform(0.0, 1.0) # a <= N < b
print random.gauss(50.0, 10.0) # gauss(mu, sigma)
print random.gammavariate(1.0, 1.0) # gammavariate(alpha, beta)
    #Gamma distribution. (Not the gamma function!) Conditions on the parameters are alpha > 0 and beta > 0.
# http://docs.python.org/lib/module-random.html
a = scipy.array([numpy.random.exponential(50.0) for i in range(100)])
print "exponential(50.0)"
print scipy.mean(a)
print scipy.std(a)
a = scipy.array([random.gauss(50.0, 10.0) for i in range(100)])
print "normal(50.0, 10.0)"
print scipy.mean(a)
print scipy.std(a)