| Off[General::spell]; |
| Off[General::spell1]; |
| iters = 6 |
| (* number of iterations - no higher than 10, unless you are very patient *) |
| SeedRandom[123332]; |
| (* Rerunning the program with the same random number seed produces the same picture. Change the seed for a different picture. |
| {a,b}=N[{4/9,2/3}]; {c,d}=N[{5/9,1/3}]; |
| (* These are the turning points *) |
| {dx[1],dy[1]}={a,b}; {dx[2],dy[2]}={c-a,d-b}; {dx[3],dy[3]}={1-c,1-d}; |
| ptlst = {{0,0},{a,b},{c,d},{1,1}}; |
| permlst = Permutations[{1,2,3}]; |
| Do[{newptlst = {},Do[{{x1,y1}=ptlst[[j]], {x2,y2}=ptlst[[j+1]], |
| delx = x2-x1, dely = y2-y1, |
| tk = Random[Integer,{1,6}],ind1=permlst[[tk]][[1]],ind2=permlst[[tk]][[2]], |
| {nx1,ny1} = {x1,y1} + {delx*dx[ind1],dely*dy[ind1]}, |
| {nx2,ny2} = {nx1,ny1} + {delx*dx[ind2],dely*dy[ind2]}, |
| newptlst = Flatten[{newptlst,{x1,y1},{nx1,ny1},{nx2,ny2}}]}, |
| {j,1,Length[ptlst]-1}],ptlst = Partition[Flatten[{newptlst,{1,1}}],2]}, |
| {i,1,iters}]; |
| Show[Graphics[Line[ptlst]],PlotRange->All,AspectRatio->1/2.5]; |