In [26]:

    

from sdm.utils import calculate_probabilities


    

In [29]:

    

bits = 256
cdf = calculate_probabilities(bits)


    

In [30]:

    

for idx, p in enumerate(cdf):
    if p > 0.001:
        break
        
print 'Radius: {} (p={:.6f})'.format(idx, cdf[idx])
print ''

for i in range(idx-5, idx+5):
    print i, cdf[i]


    

    




Radius: 103 (p=0.001067)

98 0.000106566561594
99 0.000174241902136
100 0.000280492186788
101 0.000444601537339
102 0.000693983393569
103 0.00106684558638
104 0.00161538323541
105 0.0024094567845
106 0.00354063702894
107 0.00512640372676

In [ ]: