Tupper's Function.py

from numpy import *
import matplotlib.pyplot as plt

# The inequality itself
k1 = 960939379918958884971672962127852754715004339660129306651505519271702802395266424689642842174350718121267153782770623355993237280874144307891325963941337723487857735749823926629715517173716995165232890538221612403238855866184013235585136048828693337902491454229288667081096184496091705183454067827731551705405381627380967602565625016981482083418783163849115590225610003652351370343874461848378737238198224849863465033159410054974700593138339226497249461751545728366702369745461014655997933798537483143786841806593422227898388722980000748404719
#k1 =4858450636189713423582095962494202044581400587983244549483093085061934704708809928450644769865524364849997247024915119110411605739177407856919754326571855442057210445735883681829823754139634338225199452191651284348332905131193199953502413758765239264874613394906870130562295813219481113685339535565290850023875092856892694555974281546386510730049106723058933586052544096664351265349363643957125565695936815184334857605266940161251266951421550539554519153785457525756590740540157929001765967965480064427829131488548259914721248506352686630476300

k2 = 45030838265404060658944391191720259862457481496834578271716006604162860114935889160588042689225264741351281192225004397013144833542264011622642892391575541437982998169767842389604809784878768371173619879697151174419054415191249346402577205279555771578501597358327448229878362931179478117452221136354272159568803849335470107002997400660220829047068044877606900768208365380176018080797601859084760068277158236598680389585763050536690782294968591350894328896785672781881705234432
k3 = 261272402331062461046756967766527679136125974100196816407595482011141783214027352332456136934671403365341421901634358290042403718036062783757050612374028077540573730360018887434753734699599548547323149757936323277710371752724473077367353115113372067922822562267788090289008006442830715699738122849484810499358799019299035579162473469191258178589029433829821060447640940672317988837477034296306937932267065745239580247686646439838626370116865683548097733174562354170117841595930726467690399474447226913000344596010562188879069184
K = [k1,k2,k3, k2+k3]


def tupper(x, y):
    """Return True if point (x, y) (x and y both start at 0) is to be drawn black, False otherwise
    """
    return ((k + y)//17//2**(17*int(x) + int(y)%17))%2 > 0.5

# plot the function as x*y grid
for k in K:   
    z = zeros((110,17))
    for x in range(z.shape[0]):
        for y in range(z.shape[1]):
            if tupper(x,y) != 0:
                z[x][y] =1
    plt.subplot(1,5,K.index(k)+1)
    plt.axis('off')
    plt.matshow(z,fignum=False)  
    
plt.show()