def expected_value(ratio,flips,success,failure):
P_row=mwmath.pascal_row(flips)
denom = 2**flips
expected_val = 0
for ii in range(flips+1):
ratio = P_row[ii]/denom
expected_val+=(ratio*((success)**(flips-ii))*((failure)**(ii)))
return expected_val
def expected_value(ratio, flips, success, failure):
P_row = mwmath.pascal_row(flips)
denom = 2**flips
expected_val = 0
for ii in range(flips + 1):
ratio = P_row[ii] / denom
expected_val += (ratio * ((success)**(flips - ii)) * ((failure)**(ii)))
return expected_val
def count_wins(ratio,flips,success,failure,limit):
P_row=mwmath.pascal_row(flips)
above_limit = 0.0
for ii in range(flips+1):
#s = long(success)
#if (s)**(flips-ii)>=limit:
# above_limit+=P_row[ii]
# continue
new_val=(((success)**(flips-ii))*((failure)**(ii)))
if new_val>=limit:
above_limit+=P_row[ii]
return above_limit
def count_wins(ratio, flips, success, failure, limit):
P_row = mwmath.pascal_row(flips)
above_limit = 0.0
for ii in range(flips + 1):
#s = long(success)
#if (s)**(flips-ii)>=limit:
# above_limit+=P_row[ii]
# continue
new_val = (((success)**(flips - ii)) * ((failure)**(ii)))
if new_val >= limit:
above_limit += P_row[ii]
return above_limit
