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

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