theta=3.14/6 mm <- matrix(c(cos(theta),-sin(theta), sin(theta),cos(theta)),2,2) mm=t(mm) n=200 x=runif(n) y=runif(n) xt=matrix(0,ncol=n,nrow=2) for(i in 1:n){ for(p in 1:2){ for(q in 1:2){ xt[,i]=mm%*%c(x[i],y[i]) } } } plot(x,y,xlim=c(-1.…