# -*- coding: utf-8 -*-
"""
Created on Fri Dec 10 14:18:13 2021

@author: Sim
"""
import numpy as np
import matplotlib.pyplot as plt

def normpdf(x, m, s):
    y = 1/(np.sqrt(2*np.pi)*s) * np.exp(-(x-m)**2/(2*s**2))
    return y
    

x = np.linspace(-5, 7, 100)
y1 = normpdf(x, 0, 1)
y2 = normpdf(x, 2, 1)
y3 = normpdf(x, 0, 2)
y4 = normpdf(x, 2, 2)

plt.plot(x, y1, label='N(0,1)')
plt.plot(x, y2, label='N(2,1)')
plt.plot(x, y3, label='N(0,2^2)')
plt.plot(x, y4, label='N(2,2^2)')
plt.axhline(y=0)
plt.legend()