montecarlo, Python 3

montecarlo

#python 3.5.2

from random import random

lista = [-515, 723, -444, 850, -662, -271, 335, -321, 982, -510, -737, 1012, 496, -456, 621, -682, 780, -200, -311, 1168]
rdo =[]
n=5000
ntram = 40
for x in range(0, n, 1):
    matriz=[]

for elem in lista:
        matriz.append([random(), elem])
    
    matrizorden = sorted(matriz)

sdo = 0
    for i in range(0,len(matrizorden),1):
        sdo = sdo + float(matrizorden[i][1])
        if i == 0:
            sdomax = float(matrizorden[i][1])
            mDDA = 0
        else:    
            if sdo >= sdomax:
                sdomax = sdo
            mDD = sdomax - sdo
            if mDD > mDDA:
                mDDA = mDD
    rdo.append(mDDA)

minimo = min(rdo)
maximo = max(rdo)
experimentos = len(rdo)
print('Valor\t\t', 'Frecuencia\t', 'Frec.Acum')
tramo = float(maximo-minimo)/ntram
tramanter =0
for x in range(1, ntram+1): 
    nelem = 0
    
    for valor in rdo:
        if valor<=(minimo+x*tramo):
            nelem = nelem + 1
        frecuencia=(nelem-tramanter)/experimentos
        frecacum = nelem/experimentos
        intervalo = (minimo+x*tramo)
    print('%.3f' %intervalo, '\t', '%.4f' %frecuencia, '\t', '%.4f' %frecacum)
    tramanter=nelem


run \| edit \| history \| help	0