Re: [Rd] uniroot violates bounds?

2023-02-20 Thread Serguei Sokol via R-devel
Le 18/02/2023 à 21:44, J C Nash a écrit : I wrote first cut at unirootR for Martin M and he revised and put in Rmpfr. The following extends Ben's example, but adds the unirootR with trace output. c1 <- 4469.822 c2 <- 572.3413 f <- function(x) { c1/x - c2/(1-x) }; uniroot(f, c(1e-6, 1)) uniroot(

Re: [Rd] uniroot violates bounds?

2023-02-20 Thread Martin Maechler
> Serguei Sokol via R-devel > on Mon, 20 Feb 2023 10:59:22 +0100 writes: > Le 18/02/2023 à 21:44, J C Nash a écrit : >> I wrote first cut at unirootR for Martin M and he revised and put in >> Rmpfr. >> >> The following extends Ben's example, but adds the unirootR