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 coeficientesb <- vetor das constantesx <- 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] <- 0for(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 <- 0for (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.