# by http://fractalswithr.blogspot.fr/2007/04/rossler-attractor.html

####################
#Rossler Attractor #
####################

#dx / dt = - y - z
#dy / dt = x + a y
#dz / dt = b + z ( x - c )
#
#where a = 0.2, b = 0.2, c = 5.7

rossler <- function (n = 2000,
                     a = .2, b = .2, c = 5.7,
                     x1 = 0.0001, y1 = 0.0001, z1 = 0.0001) {
  x <- vector()
  y <- vector()
  z <- vector()
  x[1] <- x1
  y[1] <- y1
  z[1] <- z1

  h <- 0.015
  for (i in 2:n){
    x[i]<-x[i-1]-h*(y[i-1]+z[i-1])
    y[i]<-y[i-1]+h*(x[i-1]+a*y[i-1])
    z[i]<-z[i-1]+h*(b+z[i-1]*(x[i-1]-c))
  }

  require(rgl)
  plot3d(x,y,z, col = heat.colors(n), size = 1)
}

rossler(20000)