Pablo Francisco Pérez Hidalgo
27,044
reputation
8
36
62
http://careers.stackoverflow.com/pablofph
from numpy import *
n = 20
a = transpose(tril(ones((n+1,n+1))))
b = copy(a)
for i in xrange(n):
b = dot(b,a)
print(b[0,n])