normal distribution

# N(0, 1), N - количество элементов
def NStandart(N):
    xk = np.array([])
    # R(0, 1)
    a1 = RStandart(N)
    time.sleep(3) # pause
    a2 = RStandart(N)

    for i in range(N):
        xk = np.append(xk, np.sqrt(-2 * np.log(a1[i])) * np.cos(2 * np.pi * a2[i]))
    return xk

N = NStandart(5000)
DrawHistogram(-5, 5, 0, 0.05, N)