Mini_Sample_Monte_Carlo_Simulation_for_Ingrid_2017.py

import numpy as np
from numpy import random as rng

samples = rng.random(500000)
prices = (np.array([10, 20, 30, 40]) * 1000).reshape(4, 1)
f = np.vectorize(lambda x: x * 10)
final = f(samples) * prices
means = final.mean(1)

print(means, np.argmax(means))
print(final.std(1))