Rotation of Circles

xxxxxxxxxx

 

1

r.svg(filename='"Rplots.svg"',width=6,height=6,onefile=True,

2

    family='"sans"',bg='"white"',pointsize=10,

3

    antialias=r.c('"default"','"none"','"gray"','"subpixel"'))

4

r.eval("""

5

rotate_xy<-function (k,x,y){

6

    i<-1; xyi<-array(c(0,0),c(2*k,2))

7

    while (i<=2*k){

8

        xyi[i,1]<-cos(i*pi/k)*x-sin(i*pi/k)*y 

9

        xyi[i,2]<-sin(i*pi/k)*x+cos(i*pi/k)*y; i<-i+1}

10

    return(xyi)}

11

k<-sample(100:900,1)/1000; n<-sample(5:15,1)

12

m<-sample(5:9,1); r<-sample(5:35,m)

13

cat(paste0('r = ',r))

14

cat(paste0('\n','k = ',k,'; n = ',n,'; m = ',m,'\n'))

15

lm<-3.5+m*k+max(scale(r,4)[3:m,1])

16

for (i in 1:m) {

17

    ci<-rotate_xy(n,.5+k*i,.5+k*i); par(mar=c(0,0,0,0))

18

    col1=rgb(runif(1),runif(1),1)

19

    col2=rgb(runif(1),runif(1),1,alpha=.1) 

20

    plot(ci[,1],ci[,2],type='p',cex=r[i],col=col1,

21

         xlab='',ylab='',xaxt='n',yaxt='n',frame=FALSE,

22

         xlim=c(-lm,lm),ylim=c(-lm,lm))

23

    par(new=TRUE)

24

    plot(ci[,1],ci[,2],type='p',pch=16,cex=r[i],

25

         col=col2,xlab='',ylab='',xaxt='n',yaxt='n',

26

         frame=FALSE,xlim=c(-lm,lm),ylim=c(-lm,lm))

27

    par(new=TRUE)}

28

 """)

29

r.dev_off()

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages