教えて下さい!
2元4次連立方程式{f(x,y)=0, g(x,y)=0}
(関数fとgはxとyの4次方程式です。)について、
例えば、[x,-10,10]のような限定された範囲で
実数解(x,y)を
Mathematicaを使って、30桁精度で数値的に求め、
それをx,yの2次元プロットしたいと思っております。
Mathematica のどのような関数を組み合わせれば
これができるか、ご教示下さいませんでしょうか?
よろしくお願いいたします。
(自分でC言語でプログラムした場合、解けるには解け
るのですが、有効数字の桁数が十分でなく、部分的に
数値が丸まってしまい、プロットがとぎれてしまう問題
があったので、有効数字を自在に調節できるMathematica
でやってみようと思ったのですが、例えば、安直に
Plot[N[Solve[{f == 0, g == 0}], 30], {x, -10, 10, 0.01}]
としてもダメでした。NRootなども検討しましたがうまく
行きません。)
No.3
- 回答日時:
#2で定数の例が悪かったので補足します.
Mathematicaはπ,eだけでなく,主要な物理定数もパッケージMiscellaneous`PhysicalConstants`に持っています.そちらを使うのもお勧めします.
No.2
- 回答日時:
Mathematicaで小数点を扱う場合はそこで有効桁数が制限されます.
なるべく少数は使わず,1/100のように分数で記述しましょう.そうすれN[]の計算まで,有効桁数が保留されます.
定数も3.1415とか記述せずに,31415 10^(-4)という具合です.
No.1ベストアンサー
- 回答日時:
何か勘違いされていると思います。
二元四次方程式には最大で16個の解があります。
NSolve[{f==0,g==0},{x,y},30]
で解を求めてその中から条件に合うものを選び出して
その値をx,yとすると
ListPlot[{x,y}]
で点を表示すればよいものと思われます。
プログラムをうまく作れば方程式を代入するだけで表示できるようになりますが、そのためにはMathematica
をよく理解しないといけません。
どうもf(x,y)=0とg(x,y)=0のグラフを描こうとされているようにも思えますが、それならば
ImplicitPlot[{f==0,g==0},{x,-10,10}]
で描くことができます。ただ標準のPackageには組み込まれていないので、使用するには下準備が必要です。
Mathematicaのヘルプをよく見て行ってください。
この回答への補足
おかげさまで望みのプロットを得ることができました。適切なアドバイス有り難うございました。
1点途中で自分のミスに気づいたのですが、fとgの関数中に外積の1成分を係数に持たせる部分があったのですが、そこの所をミスって、ベクトルのまま関数に放り込まれておりましたため、計算が先に進んでいなかったということもありました。ベクトルAの成分抽出には、A[[n]]を使いました。
根のうち、純粋に実数を抽出すると言うところまでできてはおりませんが、それは今後の課題としてさらに勉強してゆこうと思います。
大変有り難うございました。
アドバイスありがとうございました。
NSolve[{f==0,g==0},{x,y},30]と
ListPlot[{x,y}]でやってみたいと思います。
欲しい根はx,y共に実数の組み合わせですが、
これにはちょっとしたプログラミングが必要と言うことですね。
これからMathematicaを勉強させていただきます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 高校数学I 2次関数 2つの2次方程式の共通の実数解の問題についての質問です。以下の写真を見てもらえ 4 2022/05/13 11:47
- 数学 x1+3x2+2x3=4 2x1+x2-3x3=2 -5x1+5x2+18x3=a 次の連立1次方程 2 2023/07/02 03:15
- 数学 2次方程式の「(x-3)^2=4」を解くとき、 そのまま解くことも可能ですが A=x-3と置いて、A 3 2023/01/27 18:20
- 数学 高2 数2 3 2022/06/20 21:39
- 数学 『因数に分解するということ』 9 2022/06/27 06:14
- 高校 対数方程式につきまして 4 2022/05/05 07:55
- 数学 大学数学の定期テストの直しを行っているのですがこの線形代数の問題が分かりません。 次の連立一次方程式 1 2022/08/22 13:48
- 数学 数2Bの数列の問題です。 自分は、 まず数列 an=ar^(n-1)と置き こちらの問題の、y= の 1 2022/07/07 16:26
- 数学 3次方程式の解で実部が正のものが存在する条件の調べ方 0 2023/03/23 15:07
- 数学 【 数I 二次方程式の実数解 】 問題 ※写真の(2) 解答 いずれか一方のみが実数解を持つため に 1 2022/06/25 17:36
このQ&Aを見た人はこんなQ&Aも見ています
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
因数分解って何に役立つの?
-
2次関数 y=ax2+bx+cのxを求め...
-
未知数の数と必要な方程式の数...
-
高一数学 整数 〔 チャート 387...
-
何年生で習う範囲ですか?
-
与えられた2数が和と積のとき...
-
3次方程式x^3+x^2-2x-1=0の解
-
解析力学で困ってます
-
エクセルでxを求めたいのですが!
-
線形代数の平面についての問題...
-
数学の3大分野、代数・幾何・解析
-
円の方程式?円の関数じゃないの?
-
xの5乗=1 の答えを教えてく...
-
二次方程式? 2次方程式?
-
数学での文字の消去について
-
数学IIの問題です。 kを定数と...
-
3つの文字の方程式から比を出...
-
円柱と円の方程式
-
数学のアイソクライン法につい...
-
数式は私よりも賢いという意味...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報