import numpy as np
#import scipy.stats as stats
import seaborn as sns # seaborn package를 이용
import matplotlib.pyplot as plt # matplotlib.pyplot package를 이용
from numpy import random

# set the random seed:

r=1000    

np.random.seed(12)
z_1 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(34)
z_2 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(56)
z_3 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(78)
z_4 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(910)
z_5 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(21)
z_6 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(43)
z_7 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(65)
z_8 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(87)
z_9 = random.normal(loc=0.0, scale=1.0, size=r)

np.random.seed(109)
z_10 = random.normal(loc=0.0, scale=1.0, size=r)

chi_5_1 = z_1**2+z_2**2+z_3**2+z_4**2+z_5**2
chi_5_2 = z_6**2+z_7**2+z_8**2+z_9**2+z_10**2

f_5_5 = (chi_5_1/5)/(chi_5_2/5)

sns.histplot(data=f_5_5, x=None, kde=True).set(title='Histogram of F w/ df=(5,5)')
plt.show()