- 締切済み
maximaで微分方程式の数値的解法
- みんなの回答 （2）
- 専門家の回答
関連するQ&A
- maximaで複雑な微分ができない
皆様よろしくお願いいたします。 フリーの数式処理ソフトmaximaで複雑な関数を微分してみましたが、 添付画像のように多項式近似なのでしょうか。 よくわからない式が出力されてしまいます。 これは無料ソフトの限界なのでしょうか。 それとも使い方が間違っているからなのでしょうか。 ご存知の方、いらっしゃいましたらご教示ください。 画像のT(x,t)は偏微分方程式 ∂T/∂t=a*∂^2T/∂x^2 の解で、maximaで解の正しさを確認できるか試しているところ、 基本的な微分ができずに困ってます。 因みに、微分の答えは以下です。 -(q0/(kA))*erfc(x/√(4at))
- ベストアンサー
- 数学・算数
- Maximaでの方程式の解き方
お世話になります。 フリーの数式処理ソフトMaximaで、次のような方程式を解きたいのですがどのようにすればよいのか詳しい方ご教示ください。 x^5+ln(x)+1=0
- 締切済み
- その他(ソフトウェア)
- 5元連立方程式を解きたい
複雑な方程式を解きたいと考えています。 方程式は、未知数は5つ、関係式も5つあるので、 解析的に解くことが可能な5元連立方程式だと思います。 複雑、といっても、高度な関数などが登場するのではなく (せいぜい三角関数程度)、項の数がかなり多いだけです。 とりあえず手計算を試みたのですが、 項を書き下すことすら難しいような状況で、 そのまま計算を進めても、結局、解まで辿り着けたことはありません (数人で挑戦したのですが、みんな駄目でした)。 フリーの数式処理ソフト Maxima も導入して計算させてみたのですが、 エラーを吐いてしまい、やはり解は求まりませんでした (誤植等が無いか何度も確認しましたが、駄目でした)。 他の、有償のソフトならば解ける可能性はあるでしょうか? また、方程式が解析的に解けるものなのかどうかを 判断する方法などがありましたら、お教え下さい。
- ベストアンサー
- 数学・算数
- 偏微分方程式の数値解法
偏微分方程式の込み入った質問です。 2次元(x,y)の空間で2つの関数f(x,y),g(x,y)を考えます。 そこで、それぞれにラプラス方程式を立てました。 fxx+fyy = 0 (1) gxx+gyy = 0 (2) です。これは境界値問題で、差分式からSOR法を使って収束計算によって数値解を求めることができます。f, gはそれぞれ独立という形にはなります。 そこにもう1つ式が出てきました。 fxfy + gxgy = 0 (3) というものです。f,gをx,yで1回微分してできる式です。 都合3つの式が出てきました。 この数値解を求めるにはどのような方法があるでしょうか。 数値解ですから近似解です。 3つ目の拘束条件の下でのラプラス方程式とみると、ペナルティ関数とかラグランジュの未定係数法とかいろいろあるかもなと思いますが。 3つ目の式は完全に満たすというより、できるだけ満足するようにしたいというものです。 よろしくお願いします。
- 締切済み
- 数学・算数
- 高階連立常微分方程式の数値計算
4次のルンゲクッタ法を用いた数値計算を勉強しています． 1階連立常微分方程式と高階常微分方程式は理解でき，プログラムも作成することができました． 次に高階の連立常微分方程式を解こうと思ったら，頭が混乱してしまいました． 4次のルンゲクッタ法を用いて高解連立常微分方程式を解く考え方を教えて頂ければ嬉しいです． また何か良い参考書があれば教えて頂きたいと思います． よろしくお願いします．
- ベストアンサー
- 数学・算数
- maximaで方程式解の文字式のままの計算
数式処理ソフトmaximaで、連立1次方程式を解いたのち、その解を別の式の中に文字式のまま使いたい場合どのようにすればよいのでしょうか。 たとえば、 solve([x+y=2、x-y=3],[x,y]) を解いてx=2.5,y=-0.5と出ている場合に、続けてz=x+2yと記述して計算させ、結果を1.5と表示させたいということです。 普通に記述するだけでは欲しい結果が得られません。 詳しい方ご教示ください。
- 締切済み
- 数学・算数
- 連立微分方程式の解法について
Ｃ．Ｒ．ワイリー著の工業数学という本で勉強しているものです。 ミルンの方法で数値解析はできるのですが、真値が知りたくて質問しました。連立微分方程式は、以下に示すものです。 dy/dx = y^2 + xz dz/dx = x^2 + yz ただしx=0のときy=0,z=1 誰かこの方程式が解ける人は、回答よろしくお願いします。
- ベストアンサー
- 数学・算数
- 偏微分方程式の数値計算について
偏微分方程式の解の数値計算の仕方についての質問です。 例えば、偏微分方程式の解が以下のような形だったとします。 Ｔ(x,y)=Σ(A1n*x^n+A2n*x^(-n))*cos(ny) nは変数分離定数です。A1n,A2nが積分定数です。またn=1～∞とします。 今、あるxにおけるTとyの数値計算をするとします。 つまりxを定数として扱います。 そして、上式を条件に代入して解を求めるのですが、 n=3の場合の結果が欲しい場合は、 n=1,2,3のとき連立方程式をそれぞれで解いて、 数値結果は、 （n=1の結果）＋（n=2の結果）＋（n=3の結果） という考え方であってますか？ アドバイスをお願いします。 プログラムで数値計算をするのですがなかなかうまくいかずこまっています。
- ベストアンサー
- 数学・算数
- 数式処理ソフト[maple]で微分方程式を解く
数式処理ソフト[maple]を使い始めました． 連立微分方程式を解こうとしているのですが、条件をどう入れればいいか困っています．よろしくおねがいします。 dx/dt=y-2ax; dy/dt=-x; この微分方程式は手で解く事も可能ですが、練習としてmapleで数式処理したいと思っています． 手で解く際には、特性方程式を解くと-a±sqrt(a^2-1)が出てきます。 これをmapleのdsolveで解くと，(a^2-1)が正の値として計算されてしまいます． ここで質問なのですが，(a^2-1)を負の値とした場合の連立微分方程式の解法を教えていただきたいです。 よろしくおねがいします．
- ベストアンサー
- 数学・算数
