〓 Mathematica 3 〓
■ このスレッドは過去ログ倉庫に格納されています
0001132人目の素数さん
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/
0307132人目の素数さん
2006/08/01(火) 00:52:56多分そうだと思う。
初等幾何の図を書くなら、ドローソフトを使った方が良いと思います。
0308298
2006/08/01(火) 04:16:47キタ*・゜゚・*:.。..。.:*・゚(≧∀≦)゚・*:.。. .。.:*・゜゚・*!!!!
キタ*・゜゚・*:.。..。.:*・゚(≧∀≦)゚・*:.。. .。.:*・゜゚・*!!!!
キタ*・゜゚・*:.。..。.:*・゚(≧∀≦)゚・*:.。. .。.:*・゜゚・*!!!!
キタ*・゜゚・*:.。..。.:*・゚(≧∀≦)゚・*:.。. .。.:*・゜゚・*!!!!
キタ*・゜゚・*:.。..。.:*・゚(≧∀≦)゚・*:.。. .。.:*・゜゚・*!!!!
どうみても0-1区間です。
本当にありがとうございました。
0309132人目の素数さん
2006/08/01(火) 23:03:59そういう目的なら、シンデレラが良いかも。フリーじゃないけど、
書籍扱いで安い。
http://www.amazon.co.jp/gp/product/4431711201/
http://www.amazon.co.jp/gp/product/4431709665/
0311132人目の素数さん
2006/08/01(火) 23:18:460312132人目の素数さん
2006/08/02(水) 01:01:09wintpicは図がめちゃくちゃ汚い。直線が角度の取り方によってギザギザに
なってしまい、見た目が直線ではなくなってしまう。
よって、使えない。
0313132人目の素数さん
2006/08/03(木) 19:27:370314132人目の素数さん
2006/08/05(土) 09:20:30どのようなアルゴリズムを使っているか何となく分かる方いらっしゃいませんか?
0315314
2006/08/05(土) 10:08:56初期値のとり方などが難しいので、
自分でアルゴリズムを考えて作ろうとすると
答えが全部出なかったりするのです。
0316132人目の素数さん
2006/08/05(土) 12:00:10> 多項式の根はジェンキンス・トラウブ(Jenkins‐Traub)アルゴリズムに基づく.
0317132人目の素数さん
2006/08/05(土) 14:35:30ぬおっ、どうもです!
Jenkins‐Traubアルゴリズムを知りたいのですが、
載っている資料をご存知でしたら教えていただけませんか?
よろしくおねがいします。
0318132人目の素数さん
2006/08/05(土) 16:36:340319314
2006/08/06(日) 00:01:27そんくらい調べてますって。
それでも具体的なやり方や掲載されている教科書名も無いのです。
論文集の
Algorithm 419,Comm. ACM,vol.15,(1972), p. 97
にはあるみたいですけど、一般的な日本の教科書には
まだ紹介されていないのでしょうか。
さっき本屋で調べたけど見当たりませんでした。
(ベアストウ法、DKA法で大体終わってます)
0320132人目の素数さん
2006/08/06(日) 00:12:34本当に調べたか?
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
0322314
2006/08/06(日) 00:18:47値段高すぎwwww
0323132人目の素数さん
2006/08/06(日) 14:11:150324132人目の素数さん
2006/08/07(月) 13:56:19定義などもろもろを書いた file が << でゲット出来て楽なのだが、
Windows や MacOSX だとそのままではではそれが出来ない。
どうしたら良いかな?
クリックで起動したらカレントは無理としても、
どこか適当な場所を常に参照できる様にしたい。
0325132人目の素数さん
2006/08/07(月) 14:18:410326132人目の素数さん
2006/08/07(月) 16:42:21だからカレントではなく、ホームの Math とかを常にパスに含めたい、って事。
0327132人目の素数さん
2006/08/07(月) 21:40:00この辺読んで、初期化ファイルに SetDirectory["...."]みたいに入れておけばできるんじゃね?
やったことないけど。
0328132人目の素数さん
2006/08/20(日) 15:39:580329132人目の素数さん
2006/08/22(火) 20:09:270330132人目の素数さん
2006/08/22(火) 22:58:08エクセルで1回作業してから、その度数見ながらPlotで点を書いているうちに
無駄な作業をしている気がしてきたので質問。
オフィシャルブックにも出ていないんで。
0331132人目の素数さん
2006/08/23(水) 00:11:21そんなあなたに割れまてぃか
つかまっちゃってくだたい
0332132人目の素数さん
2006/08/23(水) 22:54:20<<Graphics`Graphics`
Histogram[data]
0333132人目の素数さん
2006/08/24(木) 11:41:28�dクス
実行してきます。
0334132人目の素数さん
2006/08/30(水) 17:25:510335132人目の素数さん
2006/08/30(水) 17:59:24教えてください。
0336132人目の素数さん
2006/08/31(木) 18:39:57一発でやるコマンドは知らない。線形代数パッケージを探せばあるかも。
0337132人目の素数さん
2006/09/01(金) 11:08:50どうもありがとうございます。
0338132人目の素数さん
2006/09/02(土) 21:34:36このスレッド初めて見ました。
半年以上前の投稿なのでもう解決してる事でしょうし、わざわざレスするのも失礼かと思ったのですが
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}]を{ }で囲んだだけです。
0339132人目の素数さん
2006/09/03(日) 03:05:530340132人目の素数さん
2006/09/09(土) 02:50:26Show\[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:54Show[\[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]]\)\)
0342132人目の素数さん
2006/09/09(土) 09:17:47YY= \!\({{1\/w\^2 - 1\/w\((\[ImaginaryI]\ \((\(-\((\[ImaginaryI]\ w\^9\ \ ......
の絶対値を求めようとAbs[Y] とすると
\!\({{\(Abs\)\([\)\(1\/w\^2 - 1\/w\((\[ImaginaryI]\ \((\(-\((\(\[ImaginaryI]\)\(\ \)\(w\^9\)\
と、そのまま出てきちゃうんですが、
どうやったら、ちゃんと絶対値を表示してくれるでしょうか?
0343132人目の素数さん
2006/09/11(月) 03:16:54http://www.kurokichi.net/inside/up/src/up1849.png
↑こうなるようにしたいのですが、
↓こうなってしまいます。
http://www.kurokichi.net/inside/up/src/up1848.png
初歩的な質問ですが、解決方法を教えてください。
よろしくお願いします。
0344132人目の素数さん
2006/09/11(月) 03:31:17A = a^2+2*a*b+b^2
FullSimplify[A]
とした場合,A=(a+b)^2となりますが,
A=a(a+2b+b^2)のように特定の変数でくくることはできますか?
どなたか宜しくお願いします
0345132人目の素数さん
2006/09/11(月) 12:30:31無理じゃね?
下付きのスモールエックス(x)をラージエックス(X)
にするぐらいしか思いつかん
0346132人目の素数さん
2006/09/11(月) 15:30:18そうゆう問題にはMathematicaは不適格だと思う。
LiveMath Makerを使うのが良いと思う。数式を自分の好みに合わせて変形する
のはこれが一番だと思う。
0347343
2006/09/11(月) 18:09:30無理ですか・・・。
ラージエックスを使う方法は良さそうですね。
自分でいろいろ試行錯誤した結果
http://www.kurokichi.net/inside/up/src/up1860.png
最初に↑これをやっとけば一応希望通りの動作にすることができました。
でもこの割り当て、Clear["Global`*"]をやっても消えないんだよね。
どこに残ってるんだろう。
0348132人目の素数さん
2006/09/11(月) 23:19:43そもそもその式が違うじゃねーか(w
Factor[a^2+2a b+a b^2]
=>a(a + 2b + b^2)
なんかはできるけどね。細かい指定をするなら>>346のようなソフトが
良いかと。
0350132人目の素数さん
2006/09/15(金) 07:31:240351132人目の素数さん
2006/09/15(金) 07:33:58ttp://www.varsitywave.co.jp/campaign.html
Mathematicaでググるとなかなかヒットしないもんで。
0352132人目の素数さん
2006/09/16(土) 00:46:00A=Table[Random[]/2,100] (重複がないように細工)
Position[A,Min[A]]
とした時に後で最小値の位置を使いたいのですがPositionの時点で
結果がレベル2のリストの中に埋もれてしまいますよね、{{16}}のように。
取り出すのに %[[1,1]] とか思いつくのですがどうにもキモイです。
他のすっきりした方法はないですか?
0353132人目の素数さん
2006/09/16(土) 01:13:58別にリストから取り出す必要などなくそのままでいいんですね。
0354132人目の素数さん
2006/09/30(土) 20:46:270355132人目の素数さん
2006/10/03(火) 02:43:150356132人目の素数さん
2006/10/04(水) 03:08:39"重大な数値的エラーが含まれている可能性があります"
という警告が出ても、最後まで計算して結果を返してきますよね。
この場合の計算結果というのは、やはりほとんど信用出来ないのですか?
0357132人目の素数さん
2006/10/13(金) 02:39:51極座標でパラメータを与えてベクトル場のグラフ表示をする
方法を教えて下さい
0358132人目の素数さん
2006/10/17(火) 09:57:50こんな感じで関数を定義してみたらどう?
<< 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)
が指定できます。
0359132人目の素数さん
2006/10/17(火) 22:28:30「マケマシタ6」
0360132人目の素数さん
2006/10/17(火) 23:06:07t3=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}]
という連立方程式を解きたいのですが、エラーとなってしまいます。どこがおかしいのでしょうか?
0361132人目の素数さん
2006/10/17(火) 23:19:40ヒント:=は代入,==は比較
0362132人目の素数さん
2006/10/17(火) 23:27:47t3:=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動かせないです
0363132人目の素数さん
2006/10/18(水) 00:58:22何がやりたいのかさっぱり見えないので,
解きたい連立方程式を,普通の数学の記号で書いてもらえないか?
0364132人目の素数さん
2006/10/18(水) 01:04:32t4 == 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に解かしてみても解なしとなった。
0365362
2006/10/18(水) 02:12:11t2 == t1 - u2
t1 == t0 - u1
t3 == 0
t0 == 1000000
v == 0.05*t0 + u1
v == 0.05*t1 + u2
v == 0.05*t2 + u3
すいません、式はこうでした><
0366132人目の素数さん
2006/10/18(水) 02:57:14それなら単に
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以外の文字が全て文字定数と見なされるため,
未知変数が1個,式が8本ということになって,解無しになってしまう。
また,解を近似値じゃなくて厳密値で求めたい場合は0.05を1/20に変えること。
0367132人目の素数さん
2006/10/18(水) 17:13:31グラフに値をプロットすることはできないのでしょうか?
たとえば
y=x+5のグラフにx=1のときの値をグラフ上にプロットできないのでしょうか?
0368362
2006/10/18(水) 18:04:58ほんとうにありがとうございました!
0369132人目の素数さん
2006/10/18(水) 20:19:53Clear["Global`*"]
y[x_]=x+5;
x1=1;
gridstyle={Dashing[{0.01,0.01}]};
Plot[y[x],{x,0,5},
GridLines->{{{x1,gridstyle}},{{y[x1],gridstyle}}},
Epilog->{
Text["x="<>ToString[x1],{x1,y[x1]+1},{0,0},Background->GrayLevel[1]],
Text["y="<>ToString[y[x1]],{x1+1,y[x1]},{0,0},Background->GrayLevel[1]]
}
];
0370132人目の素数さん
2006/10/18(水) 20:45:17おまけ
Clear["Global`*"]
y[x_] = x + 5;
x1 = 1;
Plot[y[x], {x, 0, 5},
Epilog -> {
{PointSize[0.02], RGBColor[1, 0, 0], Point[{x1, y[x1]}]},
Text["{ x=" <> ToString[ x1] <> ",y=" <> ToString[y[x1]] <> "}", {x1, y[x1] + 0.5}, {0, -1}]
}
];
0371132人目の素数さん
2006/10/21(土) 10:12:58ありがとうございます.
返事が遅れて申し訳ありません.
早速月曜日にでも試したいと思います.
ありがとうございました.
0372132人目の素数さん
2006/10/25(水) 00:59:44例えば
Integrate[Exp[-(x - a)^2/2c^2], {x, -Infinity, Infinity}]
とすると
If[Re[c^2] > 0, (答1), (答2)]
と、ご丁寧に条件つきで答を教えてくれるわけですが、
あらかじめこういった条件を前提として与えるにはどうすれば良いのでしょう?
Integrate[...] where Re[c^2]>0
みたいな感じのことをしたいのです。
0373132人目の素数さん
2006/10/25(水) 02:01:54極簡単な例として,y=2xを書きたいグラフの例とします。
このとき鉛直下向きを正のx,水平左向きをyの正の方向
としてグラフを書きたいのですが,どうすればいいのでしょうか?
よろしくお願いします。m(__)m
0374132人目の素数さん
2006/10/25(水) 12:15:25Integrate[Exp[-(x-a)^2/2c^2], {x, -Infinity, Infinity}, Assumptions->{c>0}]
0376132人目の素数さん
2006/10/25(水) 19:26:52なかなか難しい問題だった。
MirroringPlot[functions_, {x_, min_, max_}, (opts___)?OptionQ] :=
Module[{old, newTicks, newLines, new},
old = Plot[functions, {x, min, max}, opts, DisplayFunction -> Identity];
newTicks = AbsoluteOptions[old, Ticks][[1, 2]] /.
{position_, label_, length_, style_} -> {-position, label, length, style};
newLines = old[[1]] /.
{pointX_Real, pointY_Real} -> {-pointX, -pointY};
new = ReplacePart[old, newLines, 1] /.
{RuleDelayed[DisplayFunction, Identity] -> RuleDelayed[DisplayFunction, $DisplayFunction],
(Ticks -> _) -> (Ticks -> newTicks)};
Show[new]]
(* 使用例 *)
MirroringPlot[2x, {x, -2, 2}]
MirroringPlot[{x^2, x^4}, {x, -2, 2}, AspectRatio->Automatic]
上記のように,一度グラフィックスプリミティブを取得して,
その目盛りの出力位置の座標と,曲線のサンプリング点の座標の全てを正負反転させた
グラフィックスプリミティブを作り直し,それを出力するという方法で
関数 MirroringPlot を作成しました。
参考:http://support.wolfram.com/mathematica/graphics/2d/twoaxisgraph.html
0377376
2006/10/25(水) 19:35:41>>376 の関数 MirroringPlot は,鉛直下向きをyの正方向,
水平左向きをxの正方向とするプロットだった。
鉛直下向きをxの正方向,水平左向きをyの正方向としてプロットするには,
>>376 の定義を改変してもいいけど,プロットしたい関数の逆関数を,
必要ならば InverseFunction を使って求めて,それを >>376 の
MirroringPlot でプロットするのが簡単かな。
0378132人目の素数さん
2006/10/25(水) 20:41:160379132人目の素数さん
2006/10/25(水) 22:04:33ありがとうございます。
おもいっきり,活用&勉強(プログラムを)させていただきます。
>>378
すみませんm(__)m
0380132人目の素数さん
2006/10/31(火) 17:20:080381132人目の素数さん
2006/11/03(金) 23:53:46このサイトのように3Dプロットで複数表示させるにはどのようにすればよいのでしょうか?
0382132人目の素数さん
2006/11/04(土) 01:48:22g1 = Plot3D[x^2 + y^2, {x, -2, 6}, {y, -2, 6},
ColorFunction->(RGBColor[0, 1, #] &), DisplayFunction->Identity];
g2 = Plot3D[4(x-2) + 4(y-2) + 8, {x, -2, 6}, {y, -2, 6},
ColorFunction->(RGBColor[1, #, 0] &), DisplayFunction->Identity];
Show[g1, g2, DisplayFunction->$DisplayFunction];
0383132人目の素数さん
2006/11/06(月) 10:03:000384383
2006/11/06(月) 11:04:410386132人目の素数さん
2006/11/09(木) 10:15:28どっちがいいと思う?
0387132人目の素数さん
2006/11/12(日) 04:12:28アップグレード代を考えると,待った方がいいんじゃないかな?
0388132人目の素数さん
2006/11/13(月) 13:07:46このプレミアサービス1年付きっての買っとけ
買ってから1年以内に新verが出た場合ただでもらえる
ttp://www.hulinks.co.jp/shop/camp_mathe52v4.html
0389132人目の素数さん
2006/11/13(月) 19:11:460390132人目の素数さん
2006/11/15(水) 12:09:040391132人目の素数さん
2006/11/15(水) 12:42:550392132人目の素数さん
2006/11/15(水) 18:17:11A,Bは三角関数を含み、数個の文字を使っています。
0393132人目の素数さん
2006/11/15(水) 20:03:59Expand
0394132人目の素数さん
2006/11/15(水) 21:01:07ありがとうございました。 もし等式の形のまま右辺をExpandできる方法があれば教えてください。
0395132人目の素数さん
2006/11/15(水) 21:07:51シツモンノイミガリカイデキナイ。
0396392
2006/11/15(水) 21:57:22初心者がおさわがせしました。
0397132人目の素数さん
2006/11/15(水) 23:41:14eq=(F==(A+B)/C)
eq2=(eq[[1]]==Expand[eq[[2]]])
これでeqの右辺だけExpandした式eq2が得られるけど、こんなことする必要性がわからんw
0398132人目の素数さん
2006/11/15(水) 23:44:23まあ、まあ、まあー
そう興奮しないで
0399132人目の素数さん
2006/11/17(金) 13:31:13Mathematicaのnotebookで論文書くとかしてるんじゃね?
0400132人目の素数さん
2006/11/17(金) 13:35:5512/12(火)東京ビッグサイト
ttp://www.wolfram.co.jp/news/events/JMC2006/
04010519
2006/11/18(土) 11:52:10Mathematica5.2をjavaで遠隔操作しようとしているのですが、行き詰ってしました・・・。
そもそも、あらかじめ作っておいたnbファイルをjavaから実行させて結果を出させることはできるんでしょうか?
もし、ご存知の方がいましたら教えてください。
よろしくお願いします。
0402132人目の素数さん
2006/11/18(土) 23:53:0504030519
2006/11/19(日) 00:21:47スクリプトを調べてはみたのですが、分かりやすい物がみつかりませんでした・・・。
すいません。
スクリプトとはどんなものなのでしょうか?
0404132人目の素数さん
2006/11/19(日) 17:17:31「javaで遠隔操作」ってのは,JSP/サーブレット?それともアプレットとか?
0405132人目の素数さん
2006/11/21(火) 11:14:08今のところ、Webサーバ上からを考えています。
0406132人目の素数さん
2006/11/21(火) 12:54:17そのWebサーバがイントラネット上ならともかく,インターネット上のサーバの場合,
「nbファイルを実行して結果を得る」行為は,Mathematica のカーネルを不特定多数に使わせることになるので,それは Mathematica のライセンス的にまずい。
そのサイトに自分しか使えないような認証をかけておくなら大丈夫だろうが。
Wolfram的には,そういう用途には webMathematica を使ってくれってことなんだろうね。
■ このスレッドは過去ログ倉庫に格納されています