[R] Newton-RaphsonMethod

[Bert Gunter]

This list has a no homework policy.

Bert Gunter

"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )


On Sun, Feb 24, 2019 at 9:22 AM malika yassa via R-help <
r-help using r-project.org> wrote:

> HELLOplease I want to approximate the solution of the equation
> f(x)=x*(x-2)+log(x)=0
> for that i did this program
>
> f <- function(x){x*(x-2)+log(x)}
> x <- c(1 : 2)
> f(x)
> h <- 1e-7
> df.dx <- function(x){(f(x + h) - f(x)) / h}
> df.dx(3/2);df.dx(2)
>  newton <- function(f, tol = 1e-7, x0 = 3/2, N = 100){
>  h = 1e-7
>  i = 1; x1 = x0
>  p = numeric(N)
>  while (i <= N) {
>  df.dx = (f(x + h) - f(x)) / h
>  x1 = (x0 - (f(x0) / df.dx))
>  p[1] = x1
>  i = i + 1
>  if (abs(x1 - x0) < tol) break
>  x0 = x1
>  }
>  return(p[1 : (i-1)])
>  }
>  app <- newton(f, x0 = 3/2)
>
> but i cann't find this approximation
> please can you help me?
>
>
>
>
>
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]