def sierpinski(N):
'''Generates the Sierpinski triangle by taking the modulo-2 of each element in Pascal's triangle'''
return [([0] * (N // 2 - a // 2)) + [binomial(a, b) % 2 for b in range(a + 1)] + ([0] * (N // 2 - a // 2)) for a in range(0, N, 2)]
@interact
def _(N=slider([2 ** a for a in range(12)], label='Number of iterations', default=64), size=slider(1, 20, label='Size', step_size=1, default=9)):
M = sierpinski(2 * N)
matrix_plot(M, cmap='binary').show(figsize=[size, size])
