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: [email protected]
To: [email protected]
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
[email protected]
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. _______________________________________________
R-br mailing list
[email protected]
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.
