import matplotlib.pyplot as plt # matplotlib 모듈안에 pyplot라는 함수를 import 하는데 별칭은 plt로 지정
from scipy.stats import norm # scipy의 stats패키지로부터 norm이라는 모듈을 import 해라
import numpy as np
# 모듈을 다 패키징한 것이 패키지
avg = 148.5
std = 30
N = 1000000
population = np.random.randn(N) * std + avg
result = []
for i in range(1, 10001):
sample = np.random.choice( population, 100 ).mean()
result.append(sample)
x = np.arange(140, 160, 0.001) # 140~160까지 0.001 간격으로 숫자를 생성
y = norm.pdf(x, np.mean(result), np.std(result))
# y = norm.pdf(x, 표본 평균들의 평균, 표본평균값들의 표준편차)
# 초등학생 키의 표본 평균값들에 대한 확률 밀도함수 값이 출력 됨
plt.plot(x,y,color='blue')
plt.fill_between( x, y, interpolate = True, color = 'lime', alpha=0.4)
# plt 모듈안의 fill_between 함수를 이용해서 확률밀도 함수 그래프의 아래 영역을 색깔로 채움
# interpolate = True 를 이용해서 아래의 영역이 색깔로 채워지게 됨
# alpha는 색깔 투명도인데 0.0~1.0사이로 기술 할 수 있음
plt.show()