• ベストアンサー

ガウス積分?

次の積分を数値計算無しで計算することは可能でしょうか。   ∫x^2 × exp(-ax^2)dx   積分範囲 [k,∞) ただし、k>0です。 実験のレポートで使用したいのですが、数学には疎いので困っています。よろしくお願いします。

質問者が選んだベストアンサー

  • ベストアンサー
  • Ae610
  • ベストアンサー率25% (385/1500)
回答No.1

aの条件が書かれていないので何とも言えないが・・・ a≦0ならば積分は発散してしまうから取り敢えずa>0で考える。 ∫[k~∞](x^2*exp(-ax^2))dx =1/2a{k*exp(-ak^2)+(1/√a)*erfc(k)} (但し、erfc(k)=∫[k~∞]exp(-t^2)dt) ・・・と計算する事は出来ると思う。

_takuan_
質問者

お礼

ありがとうございました!

その他の回答 (1)

  • Ae610
  • ベストアンサー率25% (385/1500)
回答No.2

ANo.1です。 スミマセン! 下記のように訂正させて頂きます。 ∫[k~∞](x^2*exp(-ax^2))dx =(1/2a){k*exp(-ak^2)+(1/√a)*erfc(k√a)} (但し、erfc(k)=∫[k~∞]exp(-t^2)dt)

_takuan_
質問者

お礼

ありがとうございました!

_takuan_
質問者

補足

回答有難うございます。部分積分を使うのですね! しかし、最後のガウス積分は残ったままです。 やはり計算できないのでしょうか。 因みに、a>0です。

関連するQ&A

  • ガウス積分2

    昨日も質問したのですが、進展があったので改めて質問させていただきます。 範囲が以下     [k,∞)  (k>0) になるガウス積分を求めようと思い、次のように計算してみました。     I = ∫exp(-ax^2)dx     I^2 = ∬rexp(-ar^2)drdθ ここまでは定石どおりです。 積分範囲は、x>k y>k の領域になるので、     r>k    arccos(k/r)>θ>arcsin(k/r) よって、     I^2 = ∫[k→∞]rexp(-ar^2)dr ∫[arccos(k/r)→arcsin(k/r)]dθ        = ∫r{ arccos(k/r)-arcsin(k/r) }exp(-ar^2)dr 部分積分してまとめる。         = (-1/4a){ πexp(-ak^2) + 4∫[k→∞] exp(-ar^2)dr/√(1-r^2) } もう一息で計算できそうなのですが、最後の積分方法が思いつきません。分かる方居りましたら、宜しくお願いします。 又、工科系で数学には疎いので、計算ミスなのどのお叱りも是非お願いします。

  • ガウス積分みたいです。

    ある期待値(または平均値)を計算する中にでてくるんですが、∫[-∞,∞]x・exp(-ax^2)dxの積分ってどうやればいいんですか?部分積分でやると、こんがらがってしまいます。 ガウス積分なんですか? ∫[-∞,∞]x^2・exp(-ax^2)dxの積分は1/2a*(π/a)^(1/2)っていうのは、いろんなサイトや教科書にもでていますが、前者にあげたxの1乗の場合がどうしたらいいかわかりません。ガウス積分に一般式でもあるのでしょうか? 急なお願いになってしまうのですが、お願いします。

  • 指数関数と三角関数の積の積分

    ∫sin(ax)exp{-b(x-c)^2}dx, 積分範囲[-∞, ∞] ∫cos(ax)exp{-b(x-c)^2}dx, 積分範囲[-∞, ∞] これらの定積分はどうやって計算すればいいのでしょうか? 数値計算ではなくて、不定積分を導いて計算する方法を知りたいです。

  • 積分の計算

    積分I=∫[-∞→∞]exp(-ax^2)dx の計算を極座標を用いて計算するらしいのですが、 I^2=∬exp{-a(x^2+y^2)}dxdy =exp(-ar^2)rdrdθ とするまでは分かったのですが、積分範囲がわかりません。 どのようにして考えるのでしょうか。よろしくお願いします。

  • ガウス積分

    ∫[-∞→∞](x^2)*exp(-(x^2)/α^2)dxの答えが、(α^3*√π)/2 となるのですが、1/2が出てきません。 部分積分をしてガウスの法則を用いて、以下のように解きました。 ∫[-∞→∞](x^2)*exp(-(x^2)/α^2)dx =[-α^2/2*exp(-(x^2)/α^2)*x]+∫α^2*exp(-(x^2)/α^2)dx ここで、第一項はゼロになり、第二項にガウスの積分を使いました。 =α^2*√(α^2*π) =α^3*√π どこが間違っているかわかる方いたら教えてください。

  • 有限までのガウス積分

    積分範囲が(0~p)のガウス積分 ∫exp(-a x^2)dx を行いたいのですが、どうしてよいかわかりません。 教えてください。

  • erfc(x)を含む積分

    次の計算がうまく出来なくて困ってます。 ∫ {exp(x^2) * erfc(x)}^2 dx 積分範囲は、∞~a (a>0)です。 初等関数では表せないと思うので、何とか近似して数式で表現したいです。 (数値計算ではなく) よろしくお願いします。

  • ガウス分布の積分

    実際解けるのかも疑問なのですが、 ガウス分布の積分について、範囲が0から無限ならばわかります。しかし、下式のような範囲だとどうなるのかわからないのですが、わかる方教えていただけないでしょうか? ∫[0,a]exp(-x^2)dx よろしくお願いします!

  • ガウス積分公式の証明についてです!

    ガウス積分公式の証明で、 ∫dx exp(-ax^2)=(π/a)^1/2 までは分かるのですが、 多次元の場合の ∫・・・∫dx1・・・dxn exp(-Σxi Aij xj) =(π^n/2)/(detA)^1/2 (積分範囲は-∞~∞、和の範囲はij~nです。) の証明が、どのようにすれば良いのかどうしても分かりません!教科書などもかなり調べたのですが、基本的すぎるのか探し方が悪いのか、どうやっても分かりませんでした。明日までに解かなければならず、最後の頼みの綱としてここに書かせて頂きます!分かるという方いらっしゃいましたら、どうか教えてやって下さい!お願いします!

  • ガウス積分に似た式なのですが...

    ガウス積分に似た式なのですが... ∫[0→∞](x^2)*exp(-(x^4)-(x^2))dx ガウス積分を含んでいるのですが,exp()の()中が多項式になっており,部分積分が出来そうもありません. どなたか解ける方がいらっしゃったらお教えください.