chaotic_attractor.R

gistfile1.txt

rm(list =ls())
dev.off()
for(j in 1:4) {

n=1000
x=numeric(n)
y=numeric(n)
x[1] = rnorm(1,0,0.2)
y[1] = rnorm(1,0,0.2)
a=0.51
for(i in 1:(n-1)) {
  x[i+1] = 1+y[i] - a*abs(x[i])
  y[i+1] = -x[i]
}
if (j==1) {
  plot(x,y,col='red', xlim=c(-1,2), ylim=c(-2,1))
} else {
  points(x,y,col='blue')
}
}