png("cut_a_circle.png", width = 73, height = 73)<br />
rot = function(x,y,theta){<br />
x1 = sin(theta)*x + cos(theta)*y<br />
y1 = -cos(theta)*x + sin(theta)*y<br />
list(x=x1,y=y1)<br />
}<br />
x = runif(1000,-7,7)<br />
y = x + runif(1,3,5)<br />
plot(x,y, xlim = c(-14,14), ylim= c(-14,14), type = "n")<br />
n = 99 # change as you will<br />
# copy color function from Xie Laoda!<br />
color = apply(replicate(2 * n, sample(c(0:9, LETTERS[1:6]),<br />
8, replace = TRUE)), 2, function(x) sprintf("#%s", paste(x,<br />
collapse = "")))<br />
for(i in 1:n){<br />
points(rot(x,y,runif(1,-pi,pi)),col = color,pch = 19,cex = 0.2)<br />
}<br />
dev.off()
</p>