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])