#
# R
# Ilhas de Drant
# 
# 2026 março qui,5 — sex,6
# by
# j
# 

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

Tipo <- 4

i <- 0

maxL <- 300
eps <- 0.5

xmin <- -pi/2
xmax <- +pi/2
ymin <- -pi/2
ymax <- +pi/2

dx <- (xmax - xmin)/maxL
dy <- (ymax - ymin)/maxL

xk <- xmin

i <- 0

while (xk <= xmax) {

   xk <- xmin + i*dx

   yk <- ymin

   j = 0

   while ( yk <= ymax ) {

      yk <- ymin + j*dy
    
      z1 <- xk^2*sin(xk^4) + yk^3*cos(yk^6) - 1
      z2 <- xk^2*sin(xk^3) + yk^3*cos(yk^2) - 1
      z3 <- xk^2*sin(xk^4) + yk^3*cos(yk^3) - 1
      z4 <- xk^2*sin(xk^6) + yk^3*cos(yk^11) - 1
       
      z <- switch(Tipo,
                  z1,
                  z2,
                  z3,
                  z4
                 )
       
      if ( abs(z) < eps ) {
    
         x <- c(x,xk)
         y <- c(y,yk)

      }

      j <- j + 1

    }
   
    i <- i + 1
}

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