def u_beta(beta, N):
sum = 0
for m in range(1, N+1):
sum += ex(m*beta)
return sum
def F_a_N(a, N):
sum = 0
for p in prime_list:
if p > N:
break
sum += np.log(p)*ex(a*p)
return sum
