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

André, se eu entendi bem o que você quer fazer, criar uma variável (y) para receber os valores de x a cada iteração resolveria. 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* while (abs(old - x) > tol) { old <- x x <- newtonRaphson(x) # Repetir se delta x maior tol *y <- c(y, x)* } print(paste("Raiz: ", x))

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@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.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@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.
participantes (3)
-
Andre Oliveira
-
Diego Miro
-
Tiago Souza Marçal