from ppn import *
from math import *

def hline(x,y1,y2,Color):
    if(y1<y2):
        r=range(y1,y2+1)
    else:
        r=range(y2,y1+1)
    for i in r:
        set_pixel(x, i, Color)


B=(0,)*3
W=(255,)*3
G=(0,255,0)
DG=(0,128,0)

f = lambda x: -(sqrt(cos(x))*cos((ppn_w-84)*x)+sqrt(fabs(x))-0.7)*(4-x**2)**0.01*1.3*ppn_w/8
y4, y2 = f(pi/4), f(pi/2)
l=[]
for i in range(ppn_w):
    x = ((i-ppn_w//2)/ppn_w)*5
    l.append(ceil((cos(x)<0 and y2 or f(x))+y4-y2))

imax = int(input("N="))*ppn_w

clear_screen()
for i in range(ppn_w):
    for j in range(ppn_h):
        set_pixel(i, j, ppn_has_color and B or W)

show_screen()
i = 0
while imax > 0:
    hline(i%ppn_w,l[i],l[i-1],ppn_has_color and W or B)
    hline((i-5)%ppn_w,l[i-5],l[i-6],ppn_has_color and G or B)
    hline((i-(ppn_w-29))%ppn_w,l[i-(ppn_w-29)],l[i-(ppn_w-28)],ppn_has_color and DG or B)
    hline((i-(ppn_w-24))%ppn_w,l[i-(ppn_w-24)],l[i-(ppn_w-23)],ppn_has_color and B or W)
    show_screen()
    i = (i + 1) % ppn_w
    imax -= 1