import turtle import random def rule0(xold,yold,zold): x=0 y=.2*yold z=0 return (int(x),int(y),int(z)) def rule1(xold,yold,zold): x=0.85*xold y=.85*yold+.11*zold +70 z=-.1*yold+0.85*zold return (int(x),int(y),int(z)) def rule2(xold,yold,zold): x=0.31*xold-0.41*yold y=0.1*xold+0.21*yold+21 z=0.3*zold return (int(x),int(y),int(z)) def rule3(xold,yold,zold): x=-0.29*xold +0.4*yold y=0.1*xold+0.19*yold+56 z=0.3*zold return (int(x),int(y),int(z)) turtle.speed("fastest") turtle.up() x=y=z=0 rules=[rule0,rule1,rule2,rule3] for i in range(50000): r=random.randint(0,3) x,y,z = rules[r](x,y,z) turtle.goto(x-500,y-500) turtle.dot(1) input("done")