Confira se 
estão fazendo o que é necessário. Por exemplo, se x <- 1:10; o primeiro "for" retorna o vetor 0 1 2 3 4 5 6 7 8 9, e o segundo retorna o valor 9. Talvez o que você esteja interessado seja
  • for (k in 1:(length(x)-1))
  • for (k in (length(x)-1):1)

Em 29 de junho de 2015 09:47, Tiago Souza Marçal <tiagosouzamarcal@hotmail.com> escreveu:
Da uma olhada neste site e verifique se os códigos te atendem.

http://morotalab.org/Mrode2005/mme/mme.html#section0002

Att.

Tiago.

Wecsley Prates <woprates@gmail.com> escreveu:

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
 
 

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.



--
Rafael Garcia Cunha

YNWA