In [14]:
import matplotlib.pyplot as plt
from numpy import convolve, arange, linspace, sum
# x = [1/3,1/3,1/3]
x = [0.0,0.4,0.6]
y = [1]
N = 21
for count in arange(N):
y = convolve(y,x)
percent = linspace(0,100,len(y))
plt.plot(percent,y)
print(sum(y * percent))
In [ ]: