[R-br] Algoritimo newtonRaphson

[Tiago Souza Marçal]

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.0001y <- xi <- 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 em yahoo.com.br
To: r-br em 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.0001while (abs(old - x) > tol){    old <- xx <- 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 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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150504/e454082d/attachment.html>