• ベストアンサー

maximaで積分の計算

maximaで integrate(R*exp(-(5/(10^(-9)))*(R^2)-(1/2)*R), R, 0, inf); を実行したのですが、ずっと処理中のままです。 どこかまちがっていますか? それとも、この積分は解けないのでしょうか? ちなみに画像の上の式を計算しているつもりです。

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

  • ベストアンサー
  • alice_44
  • ベストアンサー率44% (2109/4759)
回答No.2

教えて! の都合で URL にカレット「 ^ 」が入らないので、 http://www.wolframalpha.com/input/?i=integrate(R*exp(-(5%2F(0.000000001))*(R*R)-(1%2F2)*R)%2C+R%2C+0%2C+inf) に変更してみたら、Maxima と同じになった(ようだ)。 どういう内部計算をしているのだろうね?

その他の回答 (2)

  • inara1
  • ベストアンサー率78% (652/834)
回答No.3

Mapleという数式処理ソフトで計算してみると、誤差関数 erf を含む結果が出ました。a = 5e-9、b = 1/2 とすれば数値解が得られます。

  • alice_44
  • ベストアンサー率44% (2109/4759)
回答No.1

Wolfram では、直ぐに答えが出るけど? http://www.wolframalpha.com/input/?i=integrate%28R*exp%28-%285%2F%2810^%28-9%29%29%29*%28R^2%29-%281%2F2%29*R%29%2C+R%2C+0%2C+inf%29

関連するQ&A

  • maximaの積分に関する質問です

    maximaで積分を行っているのですが、 integrate(x,x,sqrt(y),y);と入力してその計算結果が返ってくるはずなのですが、 defint: lower limit of integration must be real; found sqrt(y) -- an error. To debug this try: debugmode(true); とエラー表示がされてこれ以上計算が進みません。どのように対処すれば うまくいくのでしょうか。

  • 積分区間と積分するもの、具体的な計算

    F=∫x dx (0~2まで積分) があったとします。 これを F=∫lnx d(lnx)としたとき積分範囲は(0~2)のままなのでしょうか? また、この2式は同じこと(どちらも値がF)を表しているようなのですが、なぜそうなるのでしょうか? 次に具体的な計算について質問ですが F=[1/{√(2π)*ln2}] *∫(積分区間0~x) exp[-(lnx-ln1)^2/2*(ln2)^2] d(lnx) この計算はどうするのでしょうか? 積分区間0~xについては、x=0.4で計算お願いします。 ちなみに上の式は、元の式で与えられてる値は代入したもので、もとの式の形をわかって頂くために、あえて値を代入しただけで計算してまとめておりません。

  • ガウス積分ぽいんですが…

    ∫∫{x: -Inf to Inf,y: -Inf to Inf} exp(-(x^2+y^2))/(xy) dxdy  --(*)の値を求めようと思っております。以前こちらの質問で見かけたように、x=rcosθ、y=rsinθとおいてr、θについての積分を行ってみました。 (*)=∫{0 to Inf} exp(-r^2)/r dr ∫{0 to 2π} (cosθ sinθ)^-1 dθ となり、θの積分のところを置換などで求めようとしたのですが、結局どうにもできませんでした。アドバイス等いただけたら幸いです。よろしくお願いします。

  • 積分の計算

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

  • 広義積分の計算

    エルミート多項式が完全正規直交系であることを示す途中の式で、      ∫exp(-t^2/2)dt 積分区間は t∈[-∞,∞] の積分計算の方法が分かりません。どなたか教えていただけませんか?よろしくお願いします。

  • maxima複素数の指数表示の加算

    maximaで複素数を指数表示して加算を行いたいのですが、 exp[%i*(1/4)*%pi] + exp[%i*(3/4)*%pi]; と入力しても、そのままの式が出てくるだけで計算してくれません。 どうすればexp一つにできるでしょうか?

  • 積分計算について

    こんにちは、ここでは始めて質問させていただきます。 大学の講義で固体比熱の平均エネルギー<ε>の積分計算をやっているのですが、計算方法がわからず苦戦しております。 ∫∫∫m(u^2+v^2+w^2)exp{-m(u^2+v^2+w^2)/2kT}/2 dudvdw 参考書に載っている計算方法は途中式をだいぶ省略しており、私にはさっぱり理解できません。 どうかよろしくお願いします。

  • 積分の具体的な計算ができるかどうかの判定?

    たとえば、 ∫(-∞,∞)exp(-t^2)dt って具体的に計算してきれいな格好になるじゃないですか。 同じ被積分関数で積分範囲だけ変えて任意のa<bについて ∫(-a,b)exp(-t^2)dt を計算しようとしても、一般には数値計算しかできませんよね? あるいは、原始関数が見つかるような被積分関数とか、上の例のような二重積分を利用するもの、留数定理が使えるもの、特殊関数やうまく極限値が計算できるものなんかに帰着されるなど、何か特別の条件があるものじゃないと、うまく計算できませんよね?ここでいううまく計算できるっていうのは、積分の結果を表す数の表示が具体的に与えられるという意味で、収束判定とかコンピュータで計算できるとかじゃないです。具体的な計算はできないけど級数の形になっているものみたいに部分的に計算できるものは含めても(弱形?)含めなくてもいいです。 うまくいえないですけど、こんな風に積分を数値計算を使わずに計算できるかできないかを判定するための何らかの条件というのは、上に書いたような個別ケース以外にある程度知られているのでしょうか?ガロワ理論とかパンルヴェ方程式に相当するような体系化されたものはあるんでしょうか? その周辺について知られていることが何かあればおしえてください。

  • 積分の計算

     C・R・d(θ0-θ)/dt+(θe-θ)=0 CとRが一定のとき、上式を積分して  (θe-θ)=(θe-θ)exp{-t/(CR)} となるというのですが、導出過程を教えて下さい。

  • 不定積分

    ベッセル関数のPDF資料を見ていて、以下の不定積分がでてきました。 int_ { exp(-t*t) } dt = exp(-t*t) / (-2t) - int_ { exp(-t*t) / (2 * t * t) } dt ここでint_は積分記号を表し、expは指数を求める計算です。 不定積分の式 u'v = uv - int_ { u v' } を使って解こうとしたのですが、うまくいきませんでした。 よろしくお願いいたします。