23 Mai
2014
23 Mai
'14
21:22
Supondo que a função enviada no email anterior esteja correta, o
WolframAlpha dá o mesmo resultado que tu obteve:
http://www.wolframalpha.com/input/?i=integrate%28+0.23%2B%280.64-0.23%29%2F%28%281%2B%281.22*x%29%5E1.34%29%5E%281-1%2F1.34%29%29%2C+x%3D10+to+1500%29
2014-05-23 18:18 GMT-03:00 Alisson Lucrecio <alissonluc@gmail.com>:
> 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
> R-br@listas.c3sl.ufpr.br
> 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.
--
Marcus Nunes
http://marcusnunes.me/