ダブル/トリプル積分in R

Question

私は、下記の（イメージ）の方程式をRの二重積分のために解決しようとしています。私は次のポストの助けを借りて、以下のコードを使用しています： Triple integral in R (how to specifying the domain)

わからない、私はAの値を変更するかのように私は何をしないのです場合、それはまだ私に同じ結果が得られます。次のコードは、私が解決しようとしている方程式を正しく記述していますか？

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) 

おかげで、 Krina

Answer 1

fxyz影あなたは上記で定義されたグローバルwの定義に対してローカルであるw。 Aを変更すると、グローバルwに影響しますが、それはfxyzにローカルなwとはまったく関係ありません。 Aを変更しても積分値に影響を与えない理由は、積分範囲も積分限界もAに依存しないためです。

あなたの質問に答えて： "私は解決しようとしている方程式を正しく記述していますか？"：積分自体のコードは正しいようですが、fxyzのコードは明らかにそうではありません。あなたはその関数とパラメータPw, Cw, Aとの間にある程度の依存関係があるように思えますが、その依存関係が何であるべきかは説明していません。おそらく、あなたはここか別の質問のいずれかで明らかにすることができます。