Pourquoi mon double intégral ne fonctionne-t-il pas?

Question

J'ai suivi cet article ici double integral in R et j'ai changé la fonction et les limites pour faire correspondre ci-dessous, mais cela ne fonctionne pas.

InnerFunc = function(x) { x + (y^2) } 
InnerIntegral = function(z) { sapply(y, 
     function(z) { integrate(InnerFunc, 0, 2)$value }) } 
integrate(InnerIntegral, 0, 1) 

Je reçois cette erreur:

Erreur dans l'intégration (InnerFunc, 0, 2): évaluation de la fonction a donné un résultat de mauvais type

Answer 1

vos variables sont hors de question. Cela devrait faire ce que vous voulez

InnerFunc <- function(x, y) { x + (y^2) } 
InnerIntegral <- function(y) { sapply(y, 
    function(z) { integrate(InnerFunc, 0, 2, y=z)$value }) } 
integrate(InnerIntegral, 0, 1) 
# 2.666667 with absolute error < 3e-14 

Remarquez comment InnerFunc est maintenant une fonction à la fois x et y. Notez également comment nous appliquons sur le vecteur passé à InnerIntegral via le paramètre y (plutôt que d'ignorer la valeur z transmise). Nous passons également dans la valeur actuelle de y à .

Bien que vous avez saisi

InnerFunc <- function(x, y) { x + (y^2) } 

les mathématiques que vous avez dessiné les regards devraient se traduire vraiment

donc je ne suis pas sûr de ce que vous étiez vraiment après.