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