[R-br] Método de Jacobi

[Wecsley Prates]

Obrigado a todos pela ajuda.

Abs

Em 29 de junho de 2015 10:01, Rafael Garcia Cunha <rafareds em gmail.com>
escreveu:

> Corrigindo: Por exemplo, se* x <- 1:10*; *1:lenght(x)-1* retorna o vetor *0
> 1 2 3 4 5 6 7 8 9*, e *length(x)-1:1* retorna o valor *9*
>
> Em 29 de junho de 2015 09:58, Rafael Garcia Cunha <rafareds em gmail.com>
> escreveu:
>
>> 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 <
>> tiagosouzamarcal em 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 em 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 em 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
>>
>>
>
>
> --
> Rafael Garcia Cunha
>
> YNWA
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>



-- 
*                       Wecsley O. Prates*


*Doutorando em Estatística - Universidade Federal de Minas Gerais*
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150629/cb521920/attachment.html>