Data605_hw1_jr

Letters J and R

library(animation)

x <- c(rep(-1,1000),
       seq(-1.5,-1,length.out=500), 
       seq(0,1,length.out=1000), 
       rep(0, 500),
       rep(.75, 500),
       rep(1, 500),
       seq(0,1,length.out=1000))

y <- c(seq(-1,1, length.out=1000),
       rep(-1,500), 
       rep(0,1000), 
       seq(-1,1,length.out=500),
       seq(-1, 0, length.out=500),
       seq(0, 1, length.out=500),
       rep(1,1000))

z <- rbind(x,y)

plot(y~x, xlim=c(-3,2), ylim=c(-3,3))

dev.control('enable')

myani=ani.record(reset = TRUE, replay.cur = FALSE)


a=diag(2)
x11()
for (i in seq(-3,3,length.out=100)){
  a[1,1]=i
  newmat=apply(z,2,function(x) a%*%x)
  plot(newmat[2,]~newmat[1,], xlim=c(-3,3), ylim=c(-3,3))
  ani.record()
  #Sys.sleep(1)
}

a=diag(2) 
for (i in seq(-3,3,length.out=100)){
  
  a[1,2]=i
  newmat=apply(z,2,function(x) a%*%x)
  plot(newmat[2,]~newmat[1,], xlim=c(-3,3), ylim=c(-3,3))
  ani.record()
  #Sys.sleep(1)
}

a=diag(2) 
for (i in seq(-3,3,length.out=100)){
  
  a[2,1]=i
  newmat=apply(z,2,function(x) a%*%x)
  plot(newmat[2,]~newmat[1,], xlim=c(-3,3), ylim=c(-3,3))
  ani.record()
  #Sys.sleep(1)
}

a=diag(2) 
for (i in seq(-3,3,length.out=100)){
  
  a[2,2]=i
  newmat=apply(z,2,function(x) a%*%x)
  plot(newmat[2,]~newmat[1,], xlim=c(-3,3), ylim=c(-3,3))
  ani.record()
  #Sys.sleep(1)
}

myf=function(x){matrix(c(cos(x), -sin(x), sin(x), cos(x)), byrow=TRUE, nrow=2)}

for (i in seq(0,20,length.out=100)){
  a=myf(i)
  newmat=apply(z,2,function(x) a%*%x)
  plot(newmat[2,]~newmat[1,], xlim=c(-3,3), ylim=c(-3,3))
  ani.record()
  #Sys.sleep(1)
}