分子が無理関数の分数微分方程式
dw/dx =(-3c/2)(x-b)^(3a/2)/(x^2+bx)
を途中を含めて教えてください。
a, b は定数で常に x>b, x>0, b>0, a>0です。
【補足】自力では積分できなかったので以下のような近似解を求めてみました。
積分できて厳密解が得られるほうが良いですが、下記よりもっと近似の良い級数解でもOKです。
御存じの方、よろしくお願いいたします。
dw/dx =-(3c/2 )(x-b)^(3a/2)/(x^2+bx)=-(3c/2 ){x(1-b/x)}^(3a/2)/{2x^2(1+b/x)}
マクローリン展開よりx>bより1>b/x>0
(1+x)^r=1+rx+r(r-1)x^2/2+・・・+r!x^n/{r-n)!n!)}より
(1-b/x)^(3a/2)≒1+(3a/2)(-b/x)=1-3ab/2x
(1+b/x)^-1≒1+(-1)(b/x)=1-b/x
よって
dw/dx≒-(3c/2 )x^(3a/2)(1-3ab/2x)(1-b/x)/x^2
= -x^(3a/2-2)(3c/2){1-b/x-3ab/(2x)-3ab^2/(2x^2)}
w=∫(-3c/2x)x^(3a/2-2){1-b/x-3ab/(2x)-3ab^2/(2x^2)}dx
w=(3c/2)∫{-x^(3a/2-2)+(b+3ab/2)x^(3a/2-3)-3ab^2x^(3a/2-4)/2}dx
w=(3c/2){-x^(3a/2-1)/(3a-2)+(b+3ab/2)x^(3a/2-2)/(3a-4)-(ab^2/2)x^(3a/2-3)(a-2)}+C
No.1ベストアンサー
- 回答日時:
微分方程式とはいえ、右辺はwによらないので単なる積分で、
簡単に実行できますよ?
http://integrals.wolfram.com/index.jsp?expr=%28x …^%283a%2F2%29%2F%28x^2%2Bbx%29&random=false
リンクは途中で黒くなってしまいますが、
気にせず全部貼りつけてください。
もっとも結果は超幾何関数ですから、
あんまりうれしくないといえばそのとおりなんですけど。
ご回答をありがとうございます。超幾何関数がよくわかりません。2F1[-3a/2,-3a/2;1-3a/2;b/x]とは、どのような関数でしょうか。
実際にエクセルなどで計算させたいのですが、式を教えていただけると助かります。
また、検算のため、得られた解を微分して元の関数になるか確認するには御教示いただいたサイトでどのようにすればよろしいでしょうか。
No.2
- 回答日時:
役に立たない回答ですが。
>超幾何関数がよくわかりません。2F1[-3a/2,-3a/2;1-3a/2;b/x]とは、どのような関数でしょうか。
私にもわかりませんw
ある2階の微分方程式の解で、いろいろ興味深い性質を持つんでしょうけど
おそらく回答者さんの知りたい方向ではないでしょうね・・・
>実際にエクセルなどで計算させたいのですが、式を教えていただけると助かります。
エクセルで計算させるのはまず無理ではないかと。
VBAが使えるならなんとかなるかもしれません。
C言語が使えるなら、この本にアルゴリズムが出ていた気がします。
http://www.amazon.co.jp/C%E8%A8%80%E8%AA%9E%E3%8 …
>また、検算のため、得られた解を微分して元の関数になるか確認するには御教示いただいたサイトでどのようにすればよろしいでしょうか。
それも無理。
本気で検算したいなら、Mathematicaみたいな数式処理ソフトを使いましょう。
今はいい時代で、フリーの数式処理ソフトは結構たくさんあります。
グラフで出力することが目的なら数式処理ソフトにやらせるか、
数値的にルンゲクッタあたりで解くか・・・
でも級数解を出したいなら、超幾何関数が結局一番便利かもしれません。
http://ja.wikipedia.org/wiki/%E8%B6%85%E5%B9%BE% …
もともと級数で定義されてますしね
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
ベクトルの問題
-
x^3-6x-6=0 この三次方程式...
-
次の二次関数をy=(x-p)²+qの...
-
教えて下さい。
-
異なる2つの無理数の積について
-
この不定積分が解けません
-
lim{(a^x+b^x)/2}^1/x x→0 (a...
-
三重根号を簡単にする問題です...
-
数学の問題でわからないものが...
-
【数学】2√3の整数部分をa, 少...
-
6の倍数であることを証明
-
大学受験における相加平均・相...
-
ax + by (a,bは自然数で互いに...
-
アーベル群
-
ab>a+bは常に成り立つでしょうか?
-
数学の問題
-
位数6の群を分類したいです。
-
割り算における商と余り
-
他の式を利用した因数分解 x^3+...
-
2離れた奇数が互いに素なこと...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
異なる2つの無理数の積について
-
(a+b)(a2-ab+b2)=a3+b3の途中式
-
lim{(a^x+b^x)/2}^1/x x→0 (a...
-
三重根号を簡単にする問題です...
-
6の倍数であることを証明
-
「因数定理」は、いつ習います...
-
a^n+b^nの因数分解の仕方
-
高校生です! 数学の質問です。...
-
【数学】2√3の整数部分をa, 少...
-
ベクトルの問題
-
他の式を利用した因数分解 x^3+...
-
群の乗積表の作り方は?
-
異なる自然数a,bがあり, a+bを7...
-
立方完成,N乗完成は存在するの?
-
ab>a+bは常に成り立つでしょうか?
-
不等式の問題
-
楕円の変数変換
-
ax + by (a,bは自然数で互いに...
-
この問題の解き方教えてくださ...
-
2離れた奇数が互いに素なこと...
おすすめ情報