Caro Colegas,

Boa noite.

Eu estou tentando solução a integral dessa equação dvg com a função integrate, mas o resultado é muito superior ao esperado, algum colega poderia me ajudar apontando possíveis erros?

Obrigado.

> dvg <- function(x, ts, tr, a, n){
+   tr+(ts-tr)/((1+(a*x)^n)^(1-1/n))
+ }
> dvg <- function(x){
+   0.23+(0.64-0.23)/((1+(1.22*x)^1.34)^(1-1/1.34))
+ }
> dvg(10) - dvg(1500)
[1] 0.1417496
> integrate(dvg, lower = 10, upper = 1500)
412.4841 with absolute error < 0.0011

--

Alisson Lucrecio da Costa