# -*- coding: utf-8 -*-
"""
Created on Wed Jun 23 12:59:33 2021

@author: Sim
"""
from scipy.stats import norm
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

x = np.linspace(-5,5, num=100)
pdf1 = norm.pdf(x, 0, 1)
pdf2 = norm.pdf(x, 2, 1)
pdf3 = norm.pdf(x, 0, 2)

plt.plot(x, pdf1, 'k-', label='N(0,1)')
plt.plot(x, pdf2, 'r-', label='N(2,1)')
plt.plot(x, pdf3, 'b-', label='N(0,2^2)')
plt.legend()
plt.show()


quantile = norm.ppf([0.005, 0.01, 0.025, 0.5, .95, .975, .99, .995])
moments = norm.stats(0, 1,moments='mvsk')
median = norm.median(0,1)
stdev = norm.std(0,1)


xx = norm.rvs(0,1, size=1000)
xx2 = pd.Series(xx)
xx2.hist(density=True)
plt.plot(x,pdf1, 'r-')
plt.show()