[R-br] Algoritimo newtonRaphson

Andre Oliveira andreolsouza em yahoo.com.br
Sábado Maio 2 16:39:41 BRT 2015


Pessoal já tentei e não consegui.   Como colocar para exibir os valores das iterações passo a passo? E não só a estimativa da raiz. Já tentei colocar i=NULL para receber as iterações ..ma não consegui. 


newtonRaphson <-function (x){
x - ((exp(-x^2) - x^3 + 2) / (-2 * exp(-x^2) - 3 * x^2))
  }
x <- 0.5    # valor inicial 
old <- 0    # valor anterior 
tol <- 0.0001

while (abs(old - x) > tol){    
old <- x
x <- newtonRaphson(x)  # Repetir se delta x maior  tol
}

print(paste("Raiz: ", x))
   obrigado 


André Oliveira Souza. 
Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo.  IFES

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150502/8f7c807f/attachment.html>


Mais detalhes sobre a lista de discussão R-br