Oi David, não entendi o que vc questionou no último e-mail... Mas simplificando apenas escrevi a função considerando as duas variáveis X e Y da seguinte forma: f <- function(x,y) (2*exp(-2*x))*(3*exp(-3*y)) e em seguida, pedi para integrar usando a função do pacote rmutil, ou seja: int2(f, a=c(0,x), b=c(0,Inf)) Vc acha que precisa de mais alguma complementação?? Assumi que quando escrevi a função as variáveis já teriam sido definidas... X e Y... Valeu!