from pyx import *
from ent import *
from math import *
       
n = 100000
ca = canvas.canvas()

phi = (1 + sqrt(5)) / 2.0


for j in range(n):   
    i = j + 1 
    r = sqrt(i)
    theta = i * 2 * pi  / (phi*phi)
    x = cos(theta)*r
    y = sin(theta)*r      
    if i <= n:    
       
        factors = factor(i)               
        n_fac = 0
        for f in factors:
            n_fac += f[1]
            
        if(n_fac>1):       
            radius = 0.05*pow(2,(n_fac-1)/2.5)
            ca.fill(path.circle(x,y, radius))
            
           
d = document.document(pages = [document.page(ca, paperformat=document.paperformat.A4, fittosize=1)])
d.writePSfile("spiral_vogel_all.ps")