In [3]:

    

from aleato import Conmix,aleatvalue
import numpy as np
import matplotlib.pylab as plt

x = np.arange(2,22)
def pobrax():
    z=[]
    for i in range(1,8):
        z.append(i)
    for i in range(6):
        z.append(8)
    for j in range(7,0,-1):
        z.append(j)
    pj=[]
    for i in range(len(z)):
        pj.append(float(z[i])/float(104))

    return pj

def ExpectedValue(x,pobrax):
    e_x = []
    for index in range( len(x)):
        e_x.append(x[index]*pobrax[index]) 
    ex = np.sum(e_x)
    
    return ex   

    
x = np.arange(2,22)
pobrax = pobrax()

N = 1000
print ExpectedValue(x,pobrax)

g= aleatvalue(N,x,pobrax) 
%matplotlib inline
    
count, bins, ignored = plt.hist(g, 100, normed=True)
plt.plot(bins, np.ones_like(bins), linewidth=2, color='green')
plt.show()


    

    




11.5
------






    

In [ ]: