#
# R
# Ikeda Map
# 
# Domingo 22 fevereiro 2026
# by
# j
# 

x <- c()
y <- c()

tipo <- 1

a <- switch(tipo,
             4.0,
            2.16,
            2.98,
            6.41
            )

b <- switch(tipo,
             0.9,
            0.65,
            0.99,
            0.91
            )

w <- switch(tipo,
             3.5,
            4.69,
            3.26,
            3.02
            )

d <- switch(tipo,
            36.0,
           37.65,
           47.74,
           47.74
            )

xk <- 0
yk <- 0

i <- 0

maxL <- 30000

while (i<=maxL) {

   tk <- w - (d/(1 + xk*xk + yk*yk))

   xkp1 <- a + b*(xk*cos(tk) + yk*sin(tk) )
   ykp1 <- b*(xk*sin(tk) - yk*cos(tk))
    
   xk <- xkp1
   yk <- ykp1

   x <- c(x,xk)
   y <- c(y,yk)

   # cat(xk,yk,"\n")
   
   i <- i + 1
   
}

plot(x, y, pch=19, col="darkblue", cex = 0.15)