# Written by Robert Richardson for Stat 624

import numpy as np              
import pandas as pd             
import matplotlib.pyplot as plt 

n_samp = 5

# Loop
tic = time.perf_counter()
B = 100000
mus = [0]*B
meds = [0]*B
for i in range(B):
  x = np.random.normal(0,1,n.samp)
  mus[i] = np.mean(x)
  meds[i] = np.median(x)

toc = time.per_counter()
loop_time = toc - tic

# Vectorize
tic = time.perf_counter()
B = 100000
x = np.random.normal(0,1,[B,n.samp])
mus = x.mean(axis=1)
meds = np.median(x,axis=1)
toc = time.per_counter()
vec_time = toc-tic