2017-12-14 18 views
2

私の質問は「double integral in R」と似ていますが、私の場合、それぞれの限界が外側積分に依存する三重積分を行いたいと思います。それぞれの内部積分が他の積分に依存するときの積分

私は

私の問題は、それぞれの層は、すべての外側の層に依存していることであるトリプル統合のトリプルインテグラ 例を計算したいです。 、ポストからの提案に基づかことをあなたが鍵を添付していると思い

enter image description here

+0

これまでに何を試みましたか? – loki

+1

フラグを立てるチームへの質問:この質問は実際にはカテゴリが「あまりにも広すぎますか」ですか?私は、OPが彼または彼女がやろうとしたことを提示していたはずだが、問題は明らかであると思われる:関数f:R-> Rが上記の形式の三重積分を計算する。私は "広すぎる"カテゴリについて何かを見逃していますか? – storaged

+0

@storagedと合意しました。 RGG:この積分を評価するために、積分の領域を単純化し、SimplicialCubatureパッケージと統合することができます。しかし、私は開発のための答えの空間が必要です... –

答えて

0

Vectorize機能を追加することであり、我々はこのような何かを思い付く:

# Example function 
phi <- function(t){t} 

integrate(Vectorize(function(x) { 
    integrate(Vectorize(function(y) { 
    integrate(function(z) { 
     phi(x)*phi(y)*phi(z) ## function you want to calculate 
    }, -10, 6-x-y)$value  ## 2nd level interval 
    }), -5, 3-x)$value  ## 1st level interval 
}), -5,4)     ## top interval 

# 16480.2 with absolute error < 1.8e-10 

上記のアプローチ非常に一般的なので、他の関数をテストすることができます。 (答えはWolframAlphaと同じです)

+0

私はそれを試み、それは完全に動作するようです。ありがとうございます – RGG

+0

あなたはそのような可能性がある場合は、答えを受け入れてください – storaged