Não sei se eu entendi certo, mas acho que o problema é simples é só usar a
integrate duas vezes.
Talvez isso resolva seu problema.
## Primeiro fiz uma função pra integrar o y.
my.function <- function(y,x){
2*exp(-2*x)*(3*exp(-3*y))}
## A função precisa ser vetorizada eu fiz um for vc pode melhorar isso.
my.function2 <- function(x){
temp <- c()
for(i in 1:length(x)){
temp[i] = integrate(my.function, lower = c(0), upper = x[i],
x=x[i])$value}
return(temp)
}
## Depois apenas use a integrate novamente
integrate(my.function2, lower = 0, upper = Inf)
Não é muito eficiente usar uma integral dentro da outra, mas pra situação
simples parece razoável.
--
Wagner Hugo Bonat
----------------------------------------------------------------------------------------------
Department of Mathematics and Computer Science (IMADA)
University of Southern Denmark (SDU) and
Laboratório de Estatística e Geoinformação (LEG)
Universidade Federal do Paraná (UFPR)
_______________________________________________
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.
