1 2 3 4
def central_limit_theorem(mu, sigma, len): mu *= 2 f = lambda: mu * (1.0/int(mu**2) * (sigma * np.random.random(int(mu**2))).sum()) - sigma/2.0*mu + mu/2.0 return np.array([f() for _ in range(len)])