# NumPy와 matplotlib 라이브러리를 호출합니다
import numpy as np
import matplotlib.pyplot as plt

# 평균 및 표준편차를 정의합니다
mu, sigma = 0, 0.1 

# np.random.normal 함수를 사용해서 평균 0, 표준편차 0.1인 샘플 1000개를 추출합니다
s = np.random.normal(mu, sigma, 1000)

# 샘플들의 histogram을 출력합니다
# (s: 배열 혹은 배열들로 구성된 시퀀스, 
# 30: 해당 막대의 영역(bin)을 얼마나 채우는지 결정하는 변수)
count, bins, ignored = plt.hist(s, 30, normed=True)

# 샘플들을 이용해서 정규분포의 모양으로 출력합니다
# (plot(x축 데이터, y축 데이터)꼴로 사용)
plt.plot(bins, 1/(sigma * np.sqrt(2 * np.pi)) * 
    np.exp( - (bins - mu)**2 / (2 * sigma**2) ), linewidth=2, color='r')
plt.show()