André tente este adaptação:

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
y <- x
i <- 0

while (abs(old - x) > tol) {
  i <- i + 1
  old <- x
  x <- newtonRaphson(x)  # Repetir se delta x maior  tol        
       cat("\nIteração:",i,"\n")
       cat("\nRaiz:",x,"\n")       
}

Att.

Tiago.


#################################################################
 
Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES
 
Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas
 
################################################################# 



Date: Sat, 2 May 2015 19:39:41 +0000
From: andreolsouza@yahoo.com.br
To: r-br@listas.c3sl.ufpr.br
Subject: [R-br] Algoritimo newtonRaphson

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


_______________________________________________ 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.