Tupper's self-referential formula

tupper.py

import numpy
import matplotlib.pyplot as plt


def tupper(x,y):
    return 0.5 < ((y // 17) // (2**(17 * x + y % 17))) % 2

N = 4858450636189713423582095962494202044581400587983244549483093085061934704708809928450644769865524364849997247024915119110411605739177407856919754326571855442057210445735883681829823754139634338225199452191651284348332905131193199953502413758765239264874613394906870130562295813219481113685339535565290850023875092856892694555974281546386510730049106723058933586052544096664351265349363643957125565695936815184334857605266940161251266951421550539554519153785457525756590740540157929001765967965480064427829131488548259914721248506352686630476300

x, y = numpy.meshgrid(numpy.arange(107), numpy.arange(N, N + 17)[::-1])

plt.imshow(tupper(x, y))
plt.show()