'''
lecture_18_gaussian_mult.py
Gaussian distribution via repeated multiplication of "growth factors" [1, 1.1]
CS251: Data Analysis and Visualization
Oliver W. Layton
Spring 2019
'''

import numpy as np
import pymc3 as pm
import scipy.stats as stats
import matplotlib.pyplot as plt

plt.style.use(['seaborn-colorblind', 'seaborn-darkgrid'])

# Simulate picking sequence of 5 growth factors [1.0, 1.1] 50,000 times
mult = np.prod(np.random.uniform(low=1, high=1.1, size=(5, 50000)), axis=0)

pm.kdeplot(mult)
plt.xlabel('x', fontsize=14)
plt.ylabel('Probability', fontsize=14)

plt.show()