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
_______________________________________________
R-br mailing list
[email protected]
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.
