Confira se - for (k in 1:length(x)-1) - for (k in length(x)-1:1)
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 < [email protected]> 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 <[email protected]> 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 > [email protected] > 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
_______________________________________________ R-br mailing list [email protected] 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.
