cvs/fract/julia.py

#!/usr/bin/python

def plot(x,y):
    print "p", x, y

print "#PLOT 2"
print "o"

c=-1+0.25j

for x in range(0,570):
    for y in range(0,570):
        z= (x-285)/190.0
        z= z+(y-285)*(0+1j)/190.0
        i=0
        while i<499 and abs(z)<4:
            z= z**2+c
            i=i+1
        if i%2:
            plot(x,y)

print "x"

1	#!/usr/bin/python
2
3	def plot(x,y):
4	print "p", x, y
5
6	print "#PLOT 2"
7	print "o"
8
9	c=-1+0.25j
10
11	for x in range(0,570):
12	for y in range(0,570):
13	z= (x-285)/190.0
14	z= z+(y-285)*(0+1j)/190.0
15	i=0
16	while i<499 and abs(z)<4:
17	z= z**2+c
18	i=i+1
19	if i%2:
20	plot(x,y)
21
22	print "x"
23