Wednesday, March 21, 2007

Sierpinsky gasked

# Sierpinsky gasked, Chaos game -algorithm:

X11(width=6,height=6)

plot(0,0,xlim=c(0,1),ylim=c(0,1), col="white")
A=c(0,0)
B=c(1,0)
C=c(0.5,0.8)

D=sample(0:1000,2)/1000


for(i in 1:10000)
{
Corner=sample(1:3,1)
if(Corner==1){CornerS=A}
if(Corner==2){CornerS=B}
if(Corner==3){CornerS=C}
D=(D+CornerS)/2
points(D[1],D[2],pch=".",cex=1)
}

No comments: