# two-tailed t-test
h_0 = 300
h_1 = 290

n = 50
se = 20 / np.sqrt(n)

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

# acceptance region [294.456, 305.544]
# power: 0.942, beta: 0.058