トップページ⇒math
1001コメント292KB

〓 Mathematica 3 〓

■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん2005/07/15(金) 16:37:20
　　　　　　　　_....._{{ 〃
　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-<
　 　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）)
　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ.
　　　',　',　|Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　 ﾞ',
　　　 ヽ､', l:!Kノ}. 　 　 f:_.)iﾞi: ﾘ !　l　ル　
　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉ ﾞ ,ﾘ l ﾘ'´ 　 　　
.　　　　 ',|!!､ 　 　r‐┐ 　 ｀ ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣ 　
　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　 ＜　 Mathematica に関する話題はここに書いてね！
　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　 　 　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､ 　
　　/＼＼　　　　',　　 }　　 //｀ヽ

　★前スレ

〓Mathematica〓
http://science2.2ch.net/test/read.cgi/math/985023298/

〓 Mathematica 2 〓
http://science3.2ch.net/test/read.cgi/math/1078534285/
0268１３２人目の素数さん2006/06/27(火) 02:18:25
>>267
x={{x[1]},...,{x[20]}};
y={{y[1]},...,{y[20]}};
ListPlot[Flatten[Transpose[x_List]],Flatten[Transpose[y_List]]]
でいいんじゃない？
02692672006/06/27(火) 12:48:02
>>268さん
レスありがとうございます。
やってみたのですが、うまくいきませんでした・・・。
前の書き込みでは説明不足だったのですが、
x={{x[1]},...,{x[20]}};
y={{y[1]},...,{y[20]}};
のところは、私がやろうとしていたのは、
Excelのデータから抜き出したデータの羅列、例えば
x={{1},{2},{3}};
y={{3},{4},{5}};
のような感じで、番号を振ることができないのでうまくいかないのかなと思います。
ヘルプから、これらのようなxとyを、
n={{1,3}{2,4}{3,5}};
といった行列に直せば良いことは分かっていたのですが、この変換がクセモノで。。

勉強不足と言われればそれまでですが、
妙案がございましたら、レスしていただけると幸いです。
02702682006/06/27(火) 13:00:27
>>269
そういうnを得たいんだったら、
n=Transpose[{Flatten[x],Flatten[y]}]
でどう？
02712672006/06/27(火) 13:28:27
>>268さん
素早いレスありがとうございます！
希望の通りの行列ができました！！
本当に、ありがとうございます！！

もっと精進していこうと思います。
0272１３２人目の素数さん2006/06/28(水) 03:08:19
アルキメデスのらせんを書く方法を教えて下さい。
もうさっぱりです。
0273１３２人目の素数さん2006/06/28(水) 08:23:47
極座標で
r=nθ
を表現すればいい。
0274１３２人目の素数さん2006/06/28(水) 21:09:38
すいません、極座標の表現の仕方がわかりません。

x:=Cos[θ]
y:=Sin[θ]

とおいて、置換で書いた見たんですが直線になってしまいます。
02752682006/06/29(木) 03:02:00
>>272
<<Graphics`Graphics`
PolarPlot[t,{t,0,2 Pi}];
0276１３２人目の素数さん2006/06/29(木) 19:31:46
>>275
書き方がわかりました。
ありがとうございました。
0277１３２人目の素数さん2006/06/30(金) 12:55:28
等高線に対しての流線の描き方がわかりません。
よろしくお願いします。
0278１３２人目の素数さん2006/07/04(火) 12:10:11
正の整数 N に対して，√N を根号の中ができるだけ小さな整数に
なるように変形して a√b の形に直すとき，この整数 a と b の
求め方を教えて下さい．よろしくお願いします．
0279１３２人目の素数さん2006/07/04(火) 23:50:17
>>277,278
かんがれ!
0280１３２人目の素数さん2006/07/06(木) 14:23:48
>>278
root[27]で実行すれば、3root[3]とかなるが、こういうことじゃないのか？
0281１３２人目の素数さん2006/07/12(水) 22:08:19
あるセルを評価したときに結果がとても膨大になってしまうことってありますよね。
そういうときヘルプの中のセルのように、見出しだけ表示しておいて、
横のところをダブルクリックすると、中身が出てくるようにしたいのですが、
どなたかやり方をご存知の方いらっしゃいますでしょうか。
どうかよろしくお願いいたしますm(_ _)m
0282１３２人目の素数さん2006/07/13(木) 23:07:49
>>280
その 3 と 3 とを取り出すにはどうしたらいいですか？
0283１３２人目の素数さん2006/07/13(木) 23:49:23
「Print」で出力した数値の羅列を、Excelに出力する方法が分からないよ～え～ん。。。
0284１３２人目の素数さん2006/07/14(金) 01:14:15
>>282
私も初心者なので分かりませんが。
f[n_]:=If[Head[root[n]===Integer,root[n],If[Head[root[n]===Power,1,root[n][[1]] ]]
でよいはず。

root[27]//FullForm は Times[3,Power[... なので root[27][[1]]=3.
0285１３２人目の素数さん2006/07/15(土) 10:59:08
>>283
tsvとかで別ファイルに書き出したらいいのでは？
0286１３２人目の素数さん2006/07/22(土) 01:53:35
Piのリストを作りたいのですがどういう関数使えばいいですか？
3.141592653を[3,1,4,1,5,9,2,6,5,3]みたいにしたいです
0287１３２人目の素数さん2006/07/22(土) 04:37:25
List
0288１３２人目の素数さん2006/07/22(土) 14:46:38
RealDigits
0289１３２人目の素数さん2006/07/22(土) 15:58:10
対数関数のグラフの表示方法を教えてください。
0290１３２人目の素数さん2006/07/22(土) 16:12:38
Plot
0291１３２人目の素数さん2006/07/22(土) 16:17:34
LogPlot[Exp[-x] + 4 Exp[-2x], {x, 0, 6}]
で表示できないのですが…
02922862006/07/22(土) 16:24:56
>>288
どもです！
0293１３２人目の素数さん2006/07/22(土) 16:45:05
対数関数のグラフの表示方法を教えてください。
LogPlot[Exp[-x] + 4 Exp[-2x], {x, 0, 6}]
で表示できないのですが…
0294１３２人目の素数さん2006/07/22(土) 16:58:21
>>293

<< Graphics`Graphics`
LogPlot[Exp[-x] + 4 Exp[-2x], {x, 0, 6}]
0295１３２人目の素数さん2006/07/22(土) 17:22:50
>>294
ありがとう御座います！

<< Graphics`Graphics`はパッケージを読み込むってことですよね？

LogLinearPlot[60-Log[x], {x, 10^0, 10^30}]
y=60-Log[x]
縦軸y（線形）と横軸x（対数）
で表示したいのですが、できませんでした。
どこが問題なのでしょうか？
0296１３２人目の素数さん2006/07/27(木) 23:52:08
複素数の計算とかはできないのでしょうか？
たとえば|z|^2をzで微分するとconj(z)になったりとかそういうやつです。
0297１３２人目の素数さん2006/07/28(金) 18:41:40
471
0298１３２人目の素数さん2006/07/31(月) 01:56:04
「ListPlot3Dはデフォルトで各データ点のx座標とy座標の値が、1から始まる連続する整数であると解釈する」
とあるのですが（ヘルプ）、これを変更する方法をご存知の方はいらっしゃいますか？
教えて頂けるとすごくうれしいです。

具体的に問題点を指摘しておきます。
Plot3D[x^2+y^2,{x,0,1},{y,0,1}]
と同じグラフィックスをListPlot3Dで描くということが課題です。
partition=24;ListPlot3D[Table[x^2+y^2,{x,0,1,1/partition},{y,0,1,1/partition}]]
を描くと、x軸・y軸ともに1から25までの整数値になります。これを0から1までの実数値にすることが課題です。
ちなみにTableのx^2+y^2を{x,y,x^2+y^2}にすればうまくいくかなと思ったのですがダメでした。
２次元（ListPlot）では{x,x^2}とすると思ったような描画ができるのに・・・拡張性が無いですね（改善してほしいです）。
0299１３２人目の素数さん2006/07/31(月) 10:03:02
i=1->25⇔i-1=0->24⇔x=(i-1)/24=0->1
0300１３２人目の素数さん2006/07/31(月) 10:38:17
>>299
>>298が指摘してる問題点はそういうことではないだろ
具体的にListPlot3Dで課題を実現してみようとすればわかる
03012992006/07/31(月) 11:54:24
>>300
機械は苦手なもので申し訳ありません。
ところで、定義域を{x,0,24}にするのと{x,1,25}にするのでは
目盛りは一緒なのですが、値は違うようです。
たとえば{x,-24,24}なんかにすると分かりやすいですが。

目盛りを変える方法まではわからないなぁ。
0302１３２人目の素数さん2006/07/31(月) 15:50:02
>>298
いろいろ試してみたけどむずかしいね
でも「デフォルト」って書いてあるくらいだから変更できるはず

>>301
日本語でおｋ
0303１３２人目の素数さん2006/07/31(月) 22:31:08
うーん、あとはTicks->ほにゃららで目盛りをいじる
ことぐらいしか思いつかん。
0304１３２人目の素数さん2006/07/31(月) 22:53:17
あー、わかった。

Help で MeshRange をみれ
0305１３２人目の素数さん2006/08/01(火) 00:38:00
楕円曲線における解（楕円曲線暗号）の求め方ってどうやるんでしょうか？
教えてください。
0306１３２人目の素数さん2006/08/01(火) 00:40:40
教員をやってます。バージョン４です。
初等幾何で表示する図形をマセマティカで表示使用と思っているのですが、
やはり、2次元の座標上に図形の方程式を使ってえがくしかないんでしょうか？
よろしくお願いいたします。
0307１３２人目の素数さん2006/08/01(火) 00:52:56
>>306
多分そうだと思う。
初等幾何の図を書くなら、ドローソフトを使った方が良いと思います。
03082982006/08/01(火) 04:16:47
>>304
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!

どうみても0-1区間です。
本当にありがとうございました。
0309１３２人目の素数さん2006/08/01(火) 23:03:59
>>306
そういう目的なら、シンデレラが良いかも。フリーじゃないけど、
書籍扱いで安い。
http://www.amazon.co.jp/gp/product/4431711201/
http://www.amazon.co.jp/gp/product/4431709665/
03103092006/08/01(火) 23:09:59
>>309
すまん、一番目のURLのやつはソフトは付いていないようだ。
二番目のやつは付いている。
0311１３２人目の素数さん2006/08/01(火) 23:18:46
てふが良い。winTpic
0312132人目の素数さん2006/08/02(水) 01:01:09
>>311
wintpicは図がめちゃくちゃ汚い。直線が角度の取り方によってギザギザに
なってしまい、見た目が直線ではなくなってしまう。
よって、使えない。
0313１３２人目の素数さん2006/08/03(木) 19:27:37
Mathematica 5.2の解説書はありませんか？
0314１３２人目の素数さん2006/08/05(土) 09:20:30
NSolveで、虚数を含む高次式の複素解なども早く簡単に導けますが、
どのようなアルゴリズムを使っているか何となく分かる方いらっしゃいませんか？
03153142006/08/05(土) 10:08:56
多分、ニュートン法を使っていると思うのですが、
初期値のとり方などが難しいので、
自分でアルゴリズムを考えて作ろうとすると
答えが全部出なかったりするのです。
0316１３２人目の素数さん2006/08/05(土) 12:00:10
>>314
＞ 多項式の根はジェンキンス・トラウブ(Jenkins‐Traub)アルゴリズムに基づく．
0317１３２人目の素数さん2006/08/05(土) 14:35:30
>>316
ぬおっ、どうもです！
Jenkins‐Traubアルゴリズムを知りたいのですが、
載っている資料をご存知でしたら教えていただけませんか？
よろしくおねがいします。
0318１３２人目の素数さん2006/08/05(土) 16:36:34
ttp://www.google.com/search?ie=Shift_JIS&oe=Shift_JIS&q=Jenkins%81%5DTraub%83A%83%8B%83S%83%8A%83Y%83%80
03193142006/08/06(日) 00:01:27
>>318
そんくらい調べてますって。
それでも具体的なやり方や掲載されている教科書名も無いのです。
論文集の
Algorithm 419,Comm. ACM,vol.15,(1972), p. 97
にはあるみたいですけど、一般的な日本の教科書には
まだ紹介されていないのでしょうか。
さっき本屋で調べたけど見当たりませんでした。
（ベアストウ法、DKA法で大体終わってます）
0320１３２人目の素数さん2006/08/06(日) 00:12:34
>>319
本当に調べたか？
http://66.102.7.104/search?q=cache:mNFcla9WEwUJ:oshiete1.goo.ne.jp/kotaeru.php3%3Fq%3D32341+%E3%82%B8%E3%82%A7%E3%83%B3%E3%82%AD%E3%83%B3%E3%82%B9%E3%83%BB%E3%83%88%E3%83%A9%E3%82%A6%E3%83%96&hl=ja&ct=clnk&cd=3&lr=lang_ja
03213142006/08/06(日) 00:17:35
>>320
おおっ！こんなところに！！！
明日本屋に行ってみます。thx
03223142006/08/06(日) 00:18:47
ちょｗｗｗｗｗ
値段高すぎｗｗｗｗ
0323１３２人目の素数さん2006/08/06(日) 14:11:15
tanE*cosE
0324１３２人目の素数さん2006/08/07(月) 13:56:19
Unix だとカレントディレクトリがパスに入っていて、
定義などもろもろを書いた file が << でゲット出来て楽なのだが、
Windows や MacOSX だとそのままではではそれが出来ない。
どうしたら良いかな？
クリックで起動したらカレントは無理としても、
どこか適当な場所を常に参照できる様にしたい。
0325１３２人目の素数さん2006/08/07(月) 14:18:41
そもそもクリックで起動したらカレントなんてないだろ
0326１３２人目の素数さん2006/08/07(月) 16:42:21

だからカレントではなく、ホームの Math とかを常にパスに含めたい、って事。
0327１３２人目の素数さん2006/08/07(月) 21:40:00
http://documents.wolfram.com/v5/GettingStarted/System-SpecificInformation/FileLayout.ja.html
この辺読んで、初期化ファイルに SetDirectory["...."]みたいに入れておけばできるんじゃね？
やったことないけど。
0328１３２人目の素数さん2006/08/20(日) 15:39:58
保守
0329１３２人目の素数さん2006/08/22(火) 20:09:27
Ver6の情報ｸﾚ
0330１３２人目の素数さん2006/08/22(火) 22:58:08
ヒストグラムを書くコマンドってあるんですか？

エクセルで１回作業してから、その度数見ながらPlotで点を書いているうちに
無駄な作業をしている気がしてきたので質問。
オフィシャルブックにも出ていないんで。
0331１３２人目の素数さん2006/08/23(水) 00:11:21
>>322
そんなあなたに割れまてぃか

つかまっちゃってくだたい
0332１３２人目の素数さん2006/08/23(水) 22:54:20
>>330
<<Graphics`Graphics`
Histogram[data]
0333１３２人目の素数さん2006/08/24(木) 11:41:28
>>332
㌧クス

実行してきます。
0334１３２人目の素数さん2006/08/30(水) 17:25:51
816
0335１３２人目の素数さん2006/08/30(水) 17:59:24
あるベクトルの直交空間の基底を求める方法ってあるでしょうか？
教えてください。
0336１３２人目の素数さん2006/08/31(木) 18:39:57
そりゃ普通は直交条件を線型方程式として解いてそこから基底を抜き出す。
一発でやるコマンドは知らない。線形代数パッケージを探せばあるかも。
0337１３２人目の素数さん2006/09/01(金) 11:08:50
Solveで普通にやればよいということですね。
どうもありがとうございます。
0338１３２人目の素数さん2006/09/02(土) 21:34:36
>>162
このスレッド初めて見ました。
半年以上前の投稿なのでもう解決してる事でしょうし、わざわざレスするのも失礼かと思ったのですが

Clear["Global`*"]
sty1={Dashing[{0.01,0.01}]};
x1=Pi/2;
x2=3Pi/2;
Plot[Sin[x],{x,0,10},GridLines->{{{x1,sty1},{x2,sty1}},{Automatic,sty1}}];

で行けました。違いは、2行目のDashing[{0.01,0.01}]を｛　｝で囲んだだけです。
0339１３２人目の素数さん2006/09/03(日) 03:05:53
〓 Mekosujica 69 〓
0340１３２人目の素数さん2006/09/09(土) 02:50:26
\!\(TraditionalForm\`\(<< Graphics`ImplicitPlot`\[IndentingNewLine]
Show\[IndentingNewLine]
\([\[IndentingNewLine]Plot[{x, \(-x\), 0}, {x, 3, 4}], \
\[IndentingNewLine]Plot[{x, \(-x\), 0}, {x, \(-4\), \(-3\)}], \
\[IndentingNewLine]ImplicitPlot[{x\^2 + y\^2 == 1, x\^2 + y\^2 ==
4, x == 0}, {x, \(-2\), 2}, {y, \(-4\),
4}], \[IndentingNewLine]\[IndentingNewLine]Axes -> False, AspectRatio ->
1, ImageSize -> 400\[IndentingNewLine]]\)\)\)
0341こうだった2006/09/09(土) 02:51:54
\!\(TraditionalForm\`\(<< Graphics`ImplicitPlot`\[IndentingNewLine]
Show[\[IndentingNewLine]Plot[{x, \(-x\), 0}, {x, 3, 4}], \
\[IndentingNewLine]Plot[{x, \(-x\), 0}, {
x, \(-4\), \(-3\)}], \[IndentingNewLine]ImplicitPlot[{x\^2 + y\^2 == 1,
x\^2 + y\^2 == 4, x == 0}, {x, \(-2\), 2}, {
y, \(-4\), 4}], \[IndentingNewLine]\[IndentingNewLine]Axes -> False,
AspectRatio -> 1, ImageSize -> 400\[IndentingNewLine]]\)\)
0342１３２人目の素数さん2006/09/09(土) 09:17:47
非常に複雑な複素数　

YY= \!\({{1\/w\^2 - 1\/w\((\[ImaginaryI]\ \((\(-\((\[ImaginaryI]\ w\^9\ \ ......
の絶対値を求めようとAbs[Y]　とすると

\!\({{\(Abs\)\([\)\(1\/w\^2 - 1\/w\((\[ImaginaryI]\ \((\(-\((\(\[ImaginaryI]\)\(\ \)\(w\^9\)\

と、そのまま出てきちゃうんですが、

どうやったら、ちゃんと絶対値を表示してくれるでしょうか？
0343１３２人目の素数さん2006/09/11(月) 03:16:54
下付き文字を含んだシンボルを使いたいのですが、どうもうまくいきません。

http://www.kurokichi.net/inside/up/src/up1849.png
↑こうなるようにしたいのですが、
↓こうなってしまいます。
http://www.kurokichi.net/inside/up/src/up1848.png

初歩的な質問ですが、解決方法を教えてください。
よろしくお願いします。
0344１３２人目の素数さん2006/09/11(月) 03:31:17
質問です.

A = a^2+2*a*b+b^2
FullSimplify[A]

とした場合,A=(a+b)^2となりますが,
A=a(a+2b+b^2)のように特定の変数でくくることはできますか？

どなたか宜しくお願いします
0345１３２人目の素数さん2006/09/11(月) 12:30:31
>>343
無理じゃね？
下付きのスモールエックス（ｘ）をラージエックス（X）
にするぐらいしか思いつかん
0346１３２人目の素数さん2006/09/11(月) 15:30:18
>>344
そうゆう問題にはMathematicaは不適格だと思う。
LiveMath Makerを使うのが良いと思う。数式を自分の好みに合わせて変形する
のはこれが一番だと思う。
03473432006/09/11(月) 18:09:30
>>345
無理ですか・・・。
ラージエックスを使う方法は良さそうですね。

自分でいろいろ試行錯誤した結果
http://www.kurokichi.net/inside/up/src/up1860.png
最初に↑これをやっとけば一応希望通りの動作にすることができました。
でもこの割り当て、Clear["Global`*"]をやっても消えないんだよね。
どこに残ってるんだろう。
0348１３２人目の素数さん2006/09/11(月) 23:19:43
>>344
そもそもその式が違うじゃねーか（ｗ

Factor[a^2+2a b+a b^2]
=>a(a + 2b + b^2)

なんかはできるけどね。細かい指定をするなら>>346のようなソフトが
良いかと。
0349１３２人目の素数さん2006/09/12(火) 04:11:06
>>346 >>348様

レスありがとうございます.
式間違ってましたねww

LiveMath Maker,探してみます
0350１３２人目の素数さん2006/09/15(金) 07:31:24
Student版て通販では売ってないのですか。

0351１３２人目の素数さん2006/09/15(金) 07:33:58
すいませんここでいいのですね。
ttp://www.varsitywave.co.jp/campaign.html

Mathematicaでググるとなかなかヒットしないもんで。
0352１３２人目の素数さん2006/09/16(土) 00:46:00
1-11見ていきなりCCさくらの流れなのに噴いてしまった…

A＝Table[Random[]/2，100]　（重複がないように細工）
Position[A，Min[A]]

とした時に後で最小値の位置を使いたいのですがPositionの時点で
結果がレベル2のリストの中に埋もれてしまいますよね、{{16}}のように。
取り出すのに　％［[1,1]］　とか思いつくのですがどうにもｷﾓｲです。
他のすっきりした方法はないですか？
0353１３２人目の素数さん2006/09/16(土) 01:13:58
続き書いてたら自己解決しました、申し訳ない。
別にリストから取り出す必要などなくそのままでいいんですね。
0354１３２人目の素数さん2006/09/30(土) 20:46:27
アニメーション書こうとするとエラーがでるんだけど
0355１３２人目の素数さん2006/10/03(火) 02:43:15
0356１３２人目の素数さん2006/10/04(水) 03:08:39
計算途中で
"重大な数値的エラーが含まれている可能性があります"
という警告が出ても、最後まで計算して結果を返してきますよね。

この場合の計算結果というのは、やはりほとんど信用出来ないのですか？
0357１３２人目の素数さん2006/10/13(金) 02:39:51
質問です
極座標でパラメータを与えてベクトル場のグラフ表示をする
方法を教えて下さい
0358１３２人目の素数さん2006/10/17(火) 09:57:50
>>357
こんな感じで関数を定義してみたらどう？

<< Graphics`PlotField`
PlotPolarVectorField[{f_, g_}, {r_,rmin_,rmax_}, {theta_, thetamin_, thetamax_}, options__] :=
  Module[{h}, h[x_, y_] :=
   If[rmin <= Norm[{x, y}] <= rmax && thetamin <= Arg[x + y I] <= thetamax, 1, 0];
   PlotVectorField[h[x, y](f{Cos[g], Sin[g]} /.
      {r -> Norm[{x, y}],  theta -> Arg[x + y I] }),
        {x, -rmax, rmax}, {y, -rmax, rmax}, options]]

PlotPolarVectorField[{r, r+theta}, {r, 2, 5}, {theta, -Pi, Pi}, PlotPoints->20]

↑ 極座標で表された各点(r,θ)において，長さr，偏角r+θのベクトルが
生えているようなベクトル場を，円環状領域 2≦r≦5, -π≦θ≦π において
表示するという例。

オプションには PlotVecotrField と同じオプション（ScaleFactor,ScaleFunction, MaxArrowLength, ColorFunction, PlotPoints）
が指定できます。
0359１３２人目の素数さん2006/10/17(火) 22:28:30

「マケマシタ６」
0360１３２人目の素数さん2006/10/17(火) 23:06:07
Solve[{t4=0.95*t3-u3,
t3=0.95*t2-u2,
t2=0.95*t1-u1,
t4=0,t1=1000000,
u1+u2+u3=1000000
v=0.05*t1+u1=0.05*t2+u2=0.05*t3+u3},{v}]

という連立方程式を解きたいのですが、エラーとなってしまいます。どこがおかしいのでしょうか？
0361１３２人目の素数さん2006/10/17(火) 23:19:40
>>360
ヒント：=は代入，==は比較
0362１３２人目の素数さん2006/10/17(火) 23:27:47
Solve[{t4:=0.95*t3-u3,
t3:=0.95*t2-u2,
t2:=0.95*t1-u1,
t4:=0,t1:=1000000,
u1+u2+u3:=1000000
v==0.05*t1+u1=0.05*t2+u2=0.05*t3+u3},{v}]

こうですか？今家でちょっとmathematica動かせないです
0363１３２人目の素数さん2006/10/18(水) 00:58:22
>>362
何がやりたいのかさっぱり見えないので，
解きたい連立方程式を，普通の数学の記号で書いてもらえないか？
0364１３２人目の素数さん2006/10/18(水) 01:04:32
ちなみに，もし連立方程式

t4 == 0.95*t3 - u3
t3 == 0.95*t2 - u2
t2 == 0.95*t1 - u1
t4 == 0
t1 == 1000000
u1 + u2 + u3 == 1000000
v == 0.05*t1 + u1
v == 0.05*t2 + u2
v == 0.05*t3 + u3

を v について解きたいということなら，文字が8個でで式が9本だから，
普通は解なしになる。
実際，Mathematicaに解かしてみても解なしとなった。
03653622006/10/18(水) 02:12:11
t3 == t2 - u3
t2 == t1 - u2
t1 == t0 - u1
t3 == 0
t0 == 1000000
v == 0.05*t0 + u1
v == 0.05*t1 + u2
v == 0.05*t2 + u3

すいません、式はこうでした＞＜
0366１３２人目の素数さん2006/10/18(水) 02:57:14
>>365
それなら単に
Solve[{t3 == t2 - u3,
t2 == t1 - u2,
t1 == t0 - u1,
t3 == 0,
t0 == 1000000,
v == 0.05*t0 + u1,
v == 0.05*t1 + u2,
v == 0.05*t2 + u3}]
でOK.
解のリスト
{{t0 -> 1.0*10^6, t3 -> 0,
t1 -> 682791.4353687549,
u1 -> 317208.56463124516,
t2 -> 349722.44250594766,
u2 -> 333068.9928628073,
u3 -> 349722.44250594766,
v -> 367208.56463124504}}
が得られる。

ちなみに，vの値だけが欲しいからといって {v} をつけてはだめ。
v以外の文字が全て文字定数と見なされるため，
未知変数が１個，式が８本ということになって，解無しになってしまう。

また，解を近似値じゃなくて厳密値で求めたい場合は0.05を1/20に変えること。
0367１３２人目の素数さん2006/10/18(水) 17:13:31
初歩的な質問なのですが、
グラフに値をプロットすることはできないのでしょうか？
たとえば
y=x+5のグラフにx=1のときの値をグラフ上にプロットできないのでしょうか？
■ このスレッドは過去ログ倉庫に格納されています