from visual import *
from random import *


scene = display(title='Sierpinski triangle',
                width=800,height=800,center=(0.5,0.433,0),range=0.6)
        
x=0.0
y=0.0

n=0
while n<10000:
    
    RandomNumber=random()

    if RandomNumber<=0.333:
        x=0.5*x
        y=0.5*y
        
    if (RandomNumber>0.333 and RandomNumber<=0.666):
        x=0.5*x+0.5
        y=0.5*y

    if RandomNumber>0.666:
        x=0.5*x+0.25
        y=0.5*y+0.433

    dot=sphere(pos=(x,y), radius=0.002)