0

この関数をMathematica 9と統合しようとしています

NIntegrate[
 c4[u, v, w] Cos[ 2*w] , {u, -\[Infinity],
   0}, {v, -\[Infinity], u}, {w, -\[Infinity], v]

with
c4[x_, y_, z_] := 
  E^(-((x^2 + y^2 + z^2)/
    4 )) (E^((x y)/2 ) + E^((x z)/2 ) + E^((y z)/2));

このエラーがあります：

NIntegrate::ncvb: {u,v,w} = {0.0000318293,2.37147*10^-13,0.0000318293} 付近で v を 27 回再帰的に二分した後、NIntegrate は所定の精度に収束できませんでした。NIntegrate は、積分と誤差の推定値に対して -6.039603923653673*10^7640189323695 と 6.039603923653673`15.954589770191005*^7640189323695 を取得しました。>>

NIntegrate は、積分領域の外で被積分関数を評価しようとします。それはどのように可能ですか？

4

1 に答える 1

1

そうではないと思いますが、 NIntegrateは内部的に正のドメインに変数を代入しているため、エラー メッセージがわかりにくくなっています。

この変更を自分で行うと、まったく同じエラー メッセージが表示されることに注意してください。

 NIntegrate[c4[-u, -v, -w] Cos[-2 w], {u, 0, \[Infinity]}, {v,  u, \[Infinity]}, {w, v, \[Infinity]}]

有限境界で何が起こるかを見ると：

ListPlot[Table[ {s,
  NIntegrate[
   c4[u, v, w] Cos[2 w],
            {u, -s, 0}, 
             {v, -s, u}, 
             {w, -s, v}]} ,   {s, 1, 8, .25}], Joined -> True]

ここに画像の説明を入力

このことが無限に収束する可能性が低いことがわかると思います..

于 2013-06-06T20:35:39.310 に答える