f <- function(x) { 6 * x[1] * x[2]^2 }
cubature::adaptIntegrate(f, c(0, 0), c(1, 1))
g <- function(x) {
if (x[1] > x[2]) {return(0)} # set value to 0 if X > Y
return(f(x)) # else return joint pdf
}
cubature::adaptIntegrate(
g, c(0, 0), c(1, 1), tol = 0.01) # tol controls precision
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.