2004.sagews

T = RealDistribution('uniform', [0,1])

def f():
    return tan(pi*(T.get_random_element ()-1/2))
    
a = -10
b = 10

def f():
    return T.get_random_element ()
    
a = 0
b = 1

nbins = 50
n = 2


histdata = {}

for k in range (0, nbins):
    histdata[k] = 0
    
for _ in range (0, 100):
    x = 0
    for _ in range(0, n):
        x = x + f()
    x = x / n
    l = int((x-a)/(b-a)*nbins)
    if l >= 0 and l < nbins:
        histdata[l] += 1
    
mdata = []
for k in range (0, nbins):
    mdata.append (histdata[k])
    
print mdata

bar_chart(mdata)

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 3, 1, 2, 1, 2, 2, 0, 3, 8, 5, 5, 12, 13, 11, 9, 6, 4, 4, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1]