import numpy as np ws=.5 wb=0 wr=.5 y5=[] for i in range(10000): y5.append((np.exp(np.random.normal(.06,.55))*ws+wb*np.exp(np.random.normal(0,.4))+wr*1.03)*(np.exp(np.random.normal(.06,.55))*ws+wb*np.exp(np.random.normal(0,.4))+wr*1.03)*(np.exp(np.random.normal(.06,.55))*ws+wb*np.exp(np.random.normal(0,.4))+wr*1.03)*(np.exp(np.random.normal(.06,.55))*ws+wb*np.exp(np.random.normal(0,.4))+wr*1.03)*(np.exp(np.random.normal(.06,.55))*ws+wb*np.exp(np.random.normal(0,.4))+wr*1.03)) print(np.mean(y5))