# one-tailed z-test
h_0 = 0.8
h_1 = 0.75

n = 100
se = np.sqrt(h_0 * (1 - h_0) / n)

power = compute_power(h_0, h_1, se, tail="left")
print("power: %.3f, beta: %.3f"%(power, 1 - power))

# acceptance region > 0.734
# power: 0.346, beta: 0.654