import random as randomnumber
from scitools.std import *

N = [10**1, 10**2, 10**3, 10**6]
interval = range(0., 1., 10.)

for j in N:
    M = 0
    for i in range(N[j]):
        number = randomnumber.randint(interval) # random int... i want float?
        if 0.5 <= number <= 0.6:
            M += 1
    p = float(M)/N
    print "Probability with", N, "attempts:", p, "\n --- "