Plotting with Random Coefficients

xxxxxxxxxx

 

1

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

2

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

3

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

4

r.eval("""

5

a<-sample(5:11,1); b<-sample(13:19,1); c<-2*sample(3:6,1)

6

n<-sample(18:36,1); t<-seq(1,2*n+1,.1)

7

fx0<-sin(pi*t/n)-sin(a*pi*t/n)+sin(b*pi*t/n)

8

fy0<-cos(pi*t/n)+cos(a*pi*t/n)+cos(b*pi*t/n)

9

for (k in 1:c){

10

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

11

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

12

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

13

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

14

    plot(fx,fy,type='o',cex=.1,lwd=.5,

15

         col=rgb(1-.4*runif(1),.5*runif(1),0),

16

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

17

    par(new=TRUE)}

18

plot(fx0,fy0,type='o',cex=.2,lwd=10,

19

     col=rgb(1-.4*runif(1),.5*runif(1),0,alpha=.05),

20

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

21

grid(col='slategray')

22

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

23

""")

24

r.dev_off()

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

---

Messages