# -*- coding: utf-8 -*-
"""
Created on Mon Jun 28 10:41:05 2021

@author: Sim
"""
from scipy.stats import f
from scipy.stats import chi2
import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(0,20, 100)

pdf1 = chi2.pdf(x, 1)
pdf2 = chi2.pdf(x, 5)
pdf3 = chi2.pdf(x, 10)
plt.plot(x, pdf1, label='chi2(1)')
plt.plot(x, pdf2, label='chi2(5)')
plt.plot(x, pdf3, label='chi2(10)')
plt.legend()
plt.show()

x2 = np.linspace(0,10, 100)
pdf1 = f.pdf(x2, 1, 10)
pdf2 = f.pdf(x2, 5, 10)
pdf3 = f.pdf(x2, 10, 1)

plt.plot(x2, pdf1, label='f(1,10)')
plt.plot(x2, pdf2, label='f(5,10)')
plt.plot(x2, pdf3, label='f(10,1)')
plt.legend()