#
# R
# Figura Espectral
#
# Abril 2026 Qui, 2 — Sex, 3
# by
# j
# 

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

det3 <- function(a,b,c,d,e,f,g,h,ii) {
    ddd <- a*(e*ii - f*h) - b*(d*ii-f*g)+c*(d*h-e*g)
    return(ddd)
}

det4 <- function(a,b,c,d,e,f,g,h,ii,j,k,l,m,n,o,p) {
    ddd <- +a*det3(f,g,h,j,k,l,n,o,p)-b*det3(e,g,h,ii,k,l,m,o,p)+c*det3(e,f,h,ii,j,l,m,n,p)-d*det3(e,f,g,ii,j,k,m,n,o)
    return(ddd)
}

graph <- function() {

    lrmin <- -3
    lrmax <- +3
    limin <- -3
    limax <- +2

    m <- 18

    lrdelta <- (lrmax - lrmin)/m
    lidelta <- (limax - limin)/m
    
    i1min <- -31
    i1max <- +31
    i2min <- -31
    i2max <- +31

    n <- 18

    i1delta <- (i1max - i1min)/n
    i2delta <- (i2max - i2min)/n

    eps <- 10

    i1 <- i1min

    while (i1 <= i1max ) {

        i2 <- i2min

        while (i2 <= i2max ) {

            lr <- lrmin

            while (lr < lrmax) {

                li <- limin

                while (li <= limax) {

                    lbd <- lr + li*1i

                    a <- 0.3 - 0.5i - lbd
                    b <- -0 - 0.9i
                    c <- i1 
                    d <- -1i
                    
                    e <- -1i
                    f <- i2 - lbd
                    g <- 0.00
                    h <- 0.1 - 1i

                    ii <- 0.2
                    j <- 0.0
                    k <- 0.3 - lbd
                    l <- 0.2

                    m <- -0 -1.9i
                    n <- 0.2i
                    o <- 0.5 -0.7i
                    p <- -0-0.9i - lbd

                    rrr <- det4(a,b,c,d,e,f,g,h,ii,j,k,l,m,n,o,p)

                    if (Mod(rrr) <= eps) {

                        x <- c(x,lr)
                        y <- c(y,li)

                    }
                    
                    li <- li + lidelta
                    
                }
                
                lr <- lr + lrdelta

            }

            i2 <- i2 + i2delta
            
        }

        i1 <- i1 + i1delta

    }
    
    plot(x, y, main = 'Graph', pch=19, col="blue", cex = 0.1)
    
}

graph()