Double/Triple Integral dans R

Question

J'essaie de résoudre l'équation (image) mentionnée ci-dessous dans R pour une intégrale double sur une surface. J'utilise le code suivant à l'aide du message suivant: Triple integral in R (how to specifying the domain)

Je ne sais pas, si quelque chose me manque comme si je change la valeur de A, cela me donne toujours les mêmes résultats. Le code suivant décrit-il correctement l'équation que j'essaie de résoudre?

library(cubature) 

PW <- 0.0144 
Cw <- 400/(0.91*70) 
A <- 1960     
w <- as.numeric(PW, Cw, A) 

lower <- rep(0,3) 
upper <- rep(1,3) 

fxyz <- function(w) { 
    x <- w[1] 
    y <- w[2] 
    z <- w[3] 
    as.numeric(x <= y)*as.numeric(y <= z) 
} 

Int <- adaptIntegrate(f=fxyz,lowerLimit=lower,upperLimit=upper,doChecking=TRUE, 
     maxEval=2000000,absError=10e-5,tol=1e-5) 

grâce, Krina

Answer 1

Le w qui est locale à la définition de l'ombre fxyzw global défini ci-dessus. Changer A aurait un effet sur le w global, mais cela est tout à fait hors de la w qui est local à fxyz. La raison pour laquelle le changement A n'a aucun effet sur la valeur de l'intégrale est que ni l'intégrande ni les limites d'intégration ne dépendent de A. En réponse à votre question: "Le code suivant décrit-il correctement l'équation que j'essaie de résoudre?": Le code de l'intégrale elle-même semble correct, mais le code pour fxyz ne l'est apparemment pas. Vous semblez avoir une certaine dépendance entre cette fonction et les paramètres Pw, Cw, A, mais vous n'avez pas expliqué quelle devrait être cette dépendance. Peut-être pourriez-vous clarifier, ici ou dans une autre question.