var = alpha_np * (alpha_0 - alpha_np) / (alpha_0 **2 * (alpha_0+1))