Как я могу решить нелинейное уравнение в R и функции UNI root - PullRequest
0 голосов
/ 09 апреля 2020

Как я могу найти значение x, используя R для уравнения

a(x^b)+k*log(1+(x^c))+log(1-(u^(1/d)))=0 для значений

u=0.1,c=0.8,k=1.2,d=1.5,a=0.9,b=1 Они получили значение как 0,0539. Но я не получил этого ответа.

1 Ответ

2 голосов
/ 09 апреля 2020

Вот мой код и мои проверки решения

f <- function(x, a, b, c, d, k, u) {
  a * (x^b) + k * log(1 + (x^c))+ log(1 - (u^(1/d)))
}

res <- uniroot(f, interval = c(0, 1), u=0.1,c=0.8,k=1.2,d=1.5,a=0.9,b=1)
res
#> $root
#> [1] 0.08930995
#> 
#> $f.root
#> [1] 1.234389e-06
#> 
#> $iter
#> [1] 5
#> 
#> $init.it
#> [1] NA
#> 
#> $estim.prec
#> [1] 6.103516e-05

# Check solution
f(res$root, u=0.1,c=0.8,k=1.2,d=1.5,a=0.9,b=1)
#> [1] 1.234389e-06

# Plot
x <- seq(0, .1, length.out = 100)
plot(x, f(x, u=0.1,c=0.8,k=1.2,d=1.5,a=0.9,b=1))

Создано в 2020-04-09 пакетом Представить (v0.3.0)

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...