Polygon Plotting with Random Coefficients 2

xxxxxxxxxx

 

1

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

2

    family='"courier"',bg='"oldlace"',pointsize=10,

3

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

4

r.eval("""

5

a<-sample(5:11,1); b<-sample(13:19,1)

6

c<-2*sample(3:10,1); n<-sample(18:144,1)

7

t<-seq(1,2*n+1,2)

8

col=rgb(runif(1),runif(1),1,alpha=.3)

9

plot(0,0,type='n',xlim=c(-3,3),ylim=c(-3,3),

10

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

11

for (k in 1:c){

12

    fx<-sin(pi*t/n+k*2*pi/c)-

13

        sin(a*pi*t/n+k*2*pi/c)+sin(b*pi*t/n+k*2*pi/c)

14

    fy<-cos(pi*t/n+k*2*pi/c)+

15

        cos(a*pi*t/n+k*2*pi/c)+cos(b*pi*t/n+k*2*pi/c)

16

    polygon(fx,fy,lwd=1,xaxt='n',yaxt='n',frame=FALSE,

17

            border=rgb(1,runif(1),runif(1),alpha=.8),

18

            xlab='',ylab='',xlim=c(-3,3),ylim=c(-3,3))

19

    par(new=TRUE)}

20

paste0('a = ',a,'; b = ',b,'; c = ',c,'; n = ',n)

21

""")

22

r.dev_off()

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

---

Messages