jdigittl / simple.R secret

Created December 07, 2011


    
    View simple.R 
  
  
      simple.R
    
    
        Raw
      
    
    File suppressed. Click to show.
  
  
            1
            2
            3
            4
            5
            6
            7
            8
            9
            10
            11
            12
            13
            14
            15
            16
            17
            18
            19
            20
          
          
            wave <- function (a, R, p, w) {
    x <- (0.5 + 0.5*R*sin(p+w*a))*cos(a)
    y <- (0.5 + 0.5*R*sin(p+w*a))*sin(a)
    return (cbind(x,y))
}
 
blue <- '#3F8492' 
red <- '#E57660'
 
t <- seq(0, 2*pi, length.out=100)
M<-0.8
plot(c(-1, 1)*M, c(-1, 1)*M, pch='', bty='n')
LWD=7
pinc=2*pi/5
amp <- 0.2
lines(wave(t,amp,0,2), lwd=LWD, col=sprintf('%s%2x', blue, 1.0 * 255))
lines(wave(t,amp,pinc*1,2), lwd=LWD, col=sprintf('%s%2x', red, 1.0 * 255))
lines(wave(t,amp,pinc*2,2), lwd=LWD, col=sprintf('%s%2x', blue, 0.8 * 255))
lines(wave(t,amp,pinc*3,2), lwd=LWD, col=sprintf('%s%2x', red, 0.8 * 255))
lines(wave(t,amp,pinc*4,2), lwd=LWD, col=sprintf('%s%2x', blue, 0.6 * 255))