Olá a todos... Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações lineares. Encontrei uma programação, mas está dando erro. Essa programação me fornece um erro do tipo Erro em a[i, j] <- a[i, j] - m * a[k, j] : substituto tem comprimento zero. Como resolvo??? Agradeço a atenção a = matriz dos coeficientes b <- vetor das constantes x <- chute inicial do vetor das variáveis for (k in 1:length(x)-1){ for (i in k+1:length(x)){ m <- a[i,k]/a[k,k] a[i,k] <- 0 for(j in k+1:length(x)){ a[i,j] <- a[i,j]-m*a[k,j] b[i] <- b[i]-m*b[k] } } } x[length(x)] <- b[length(x)]/a[length(x),length(x)] for (k in length(x)-1:1){ s <- 0 for (j in k+1:length(x)){ s <- s+a[k,j]*x[j] } x[k] <- (b[k]-s)/a[k,k] } -- * Wecsley O. Prates* *Doutorando em Estatística - Universidade Federal de Minas Gerais*