from math import *
from pylab import *

for k in range (0,25):
    x1=8*(sqrt(3)/2)**k*cos(-k*pi/6)
    y1=
    x2=
    y2=
    a=[x1,x2]
    b=[y1,y2]
    plot(a,b,'red')
show()