cvs/fract/koch.py

#!/usr/bin/python

base = [(142,142), (426,142), (285,389), (142,142)]

def normal(kord1,kord2):
    x1,y1= kord1
    x2,y2= kord2
    x3= 0.5*(x2-x1+(y2-y1)*pow(3,0.5))+x1
    y3= 0.5*(y2-y1-(x2-x1)*pow(3,0.5))+y1
    return x3, y3

def transform(inbase):
    outbase= []
    for it in range(len(inbase)-1):
        (x1,y1)= inbase[it]
        (x2,y2)= inbase[it+1]
        xplus, yplus= (x2-x1)/3, (y2-y1)/3
        outbase.append(inbase[it])
        outbase.append((x1+xplus,y1+yplus))
        outbase.append(normal((x1+xplus,y1+yplus),(x1+xplus*2,y1+yplus*2)))
        outbase.append((x1+xplus*2,y1+yplus*2))
    outbase.append(inbase[-1])
    return outbase

def draw(inbase):
    print "#PLOT 2"
    print "o"
    print "e"
    for it in range(len(inbase)-1):
        x1,y1= inbase[it]
        x2,y2= inbase[it+1]
        print "l", x1, y1, x2, y2
    print "x"
    
for i in range(5):
    base= transform(base)

draw(base)
1	teddy	1.1	#!/usr/bin/python
2
3			base = [(142,142), (426,142), (285,389), (142,142)]
4
5			def normal(kord1,kord2):
6			x1,y1= kord1
7			x2,y2= kord2
8			x3= 0.5(x2-x1+(y2-y1)pow(3,0.5))+x1
9			y3= 0.5(y2-y1-(x2-x1)pow(3,0.5))+y1
10			return x3, y3
11
12			def transform(inbase):
13			outbase= []
14			for it in range(len(inbase)-1):
15			(x1,y1)= inbase[it]
16			(x2,y2)= inbase[it+1]
17			xplus, yplus= (x2-x1)/3, (y2-y1)/3
18			outbase.append(inbase[it])
19			outbase.append((x1+xplus,y1+yplus))
20			outbase.append(normal((x1+xplus,y1+yplus),(x1+xplus2,y1+yplus2)))
21			outbase.append((x1+xplus2,y1+yplus2))
22			outbase.append(inbase[-1])
23			return outbase
24
25			def draw(inbase):
26			print "#PLOT 2"
27			print "o"
28			print "e"
29			for it in range(len(inbase)-1):
30			x1,y1= inbase[it]
31			x2,y2= inbase[it+1]
32			print "l", x1, y1, x2, y2
33			print "x"
34
35			for i in range(5):
36			base= transform(base)
37
38			draw(base)