〓 Mathematica 2 〓
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
NGNG, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .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/
0577１３２人目の素数さん
NGNG5.1でたのか、知らんかった。
アップグレードのお知らせとか送ってくれないのかな。
ところで5.0からのアップグレードはいくらくらいぼったくられるのだろう？
0578１３２人目の素数さん
NGNG留年しろ
0579１３２人目の素数さん
NGNGたとえば、あまり一般的でない記法とかを自由に記述できたりしますか？
(どちらがきれいな表示かは問題にしていません。)
0580伊丹公理
NGNG編集ならLaTeX
自由自在、縦横無尽に使うなら TeX
0581１３２人目の素数さん
NGNG普通にそれぞれの用途を考えたらわかると思うが？
0582１３２人目の素数さん
NGNG0583１３２人目の素数さん
NGNG0584１３２人目の素数さん
NGNG(対応する括弧を強調表示する機能)は無いでしょうか?
0585１３２人目の素数さん
NGNGhttp://homepage2.nifty.com/autumn-soft/
0586１３２人目の素数さん
NGNGEdit -> Check Balance [Shitf Ctrl B]
0587１３２人目の素数さん
NGNGありがとう。
自動的に強調表示してくれる機能は無いのね。残念。
0588１３２人目の素数さん
NGNG面白いね。
0589１３２人目の素数さん
NGNGMathematicaで置換するやり方はわかるのですが、
テキストファイルなどと違って置換する速度がかなり遅いです。
数個程度の置換ならば良いのですが、置換する数が数万個単位になるような
長い数式だと普通の置換のやり方ではさすがに時間がかかりすぎてしまいます。
今回困っているパターンについて
例えば数式で ax+by に対して"ax+by"のように「"」がついています
これではax+by自体を文字列とみなしてしまい中の計算は行ってくれません
Mathematica上で「"」→「」の置換を行えば問題は解決するのですが、これが数万個単位であるのです
特別なテキスト形式で保存しても文字列と数列では単純に「"」のあるなしではないので
「"」だけを削除すればよいという問題でもありません
よろしければ
Mathematica上で置換をテキストレベルに早く行うことができるやり方
もしくはMathematicaの表記どおりに「"」を削除を行うだけで解決するようなやり方
知っている方がいればお願いします
ツールの有無、シェアの有無にはこだわりません
0590１３２人目の素数さん
NGNG数万個どうやって作ったわけ？
手入力してしまったとかいうのでなければ、そのデータ入力のところからやりなおすのがいい。
参考>>402,>>404。
試してみるとバグだと言いたくなる遅さだね。フロントエンドを表示せずにはNotebook操作系の関数も使えないのが痛い。
0591１３２人目の素数さん
NGNGi <= C\_A + 1, \(i++\), \[IndentingNewLine]Which[\[IndentingNewLine]i ==
1, DM[1, 1] = \(-\((C\_B*\[Mu]\_B)\)\);
For[l = 2, l <= C\_A + 1, \(l++\),
DM[1, l] = 0], \[IndentingNewLine]i == \((C\_A + 1)\),
For[l = 1, l <= C\_A - 1, \(l++\), DM[C\_A + 1, l] = 0];
DM[C\_A + 1, C\_A] = C\_A*\[Mu]\_A;
DM[C\_A + 1, C\_A + 1] = \(-\[Mu]\), \[IndentingNewLine]i\ != 1\ && \
i\ != \((C\_A + 1)\), \
For[l = 1, l <= i - 2, \(l++\), DM[i, l] = 0];
DM[i, i - 1] = \((i - 1)\) \[Mu]\_A;
DM[i, i] = \(-\((C\_B*\[Mu]\_B + \((i - 1)\)*\[Mu]\_A)\)\);
For[l = i + 2, l <= C\_A + 1, \(l++\),
DM[i, l] = 0]]]\[IndentingNewLine]
DM = Array[DM, {C\_A + 1, C\_A + 1}]\)
のどこが悪いですか
0592１３２人目の素数さん
NGNGDM2 = Array[DM, {C\_A + 1, C\_A + 1}]\)
にしたらわかる。
0593１３２人目の素数さん
NGNGNDSolve[
{Derivative[1][x][t]==-3(x[t]-y[t]),
Derivative[1][y][t]==-x[t]z[t]+27x[t]-y[t],
Derivative[1][z][t]==x[t]y[t]-z[t],
x[0]==0,z[0]==0,y[0]==2.1},
{x,y,z},
{t,0,100},
MaxSteps->5000,
Method->ExplicitEungeKutta];
ParametricPlot3D[
Evaluate[{x[t],y[t],z[t]}/.solution],
{t,0,100},
PlotPoints->5000];
でローレンツ方程式の解を3Dで描いたのですが、
この立体図形のy[t]=0での断面を描くにはどうしたら良いでしょうか。
Evaluate[{x[t],y[t]=0,z[t]}/.solution],
とやったら、断面にはならずにy方向無限遠点から見たような図になってしまいました。
0594１３２人目の素数さん
NGNGこのような非線形の等式をpについて解いて、
p=の格好にしたいのですが、
こういう式変形をさせるにはどうしたらいいのでしょうか？
教えて下さい。
0595550
NGNG遅くなりましたが，ありがとうございます．コードの意味も大体分かりました．
0596１３２人目の素数さん
NGNGあの枠を出さない方法を教えてつかあさい。
0597１３２人目の素数さん
NGNGMathematicaで陰関数を表示する方法はありますか？
GCalcやfunctionViewなどでできるんですが、
Mathematicaでもできるのかなと思いまして。
0598１３２人目の素数さん
NGNGFirst move d and y to the left hand side.
second, take the power of -1/c and then move a to the left hand side
0599伊丹公理
NGNG(log x)/x の重複値を調べる
0600１３２人目の素数さん
NGNGお願いします。
x[t],y[t],z[t]の式出してそこにy[t]の値を代入、とかできるのかな？
何にせよ今mathematica使えないからなぁ…
0601１３２人目の素数さん
NGNG返信ありがとうございます。
でも、もちろん、私も手では計算できるんですけど、
もっと実はもっと複雑な非線形の等式で似たような事を
したいので、それをmathematicaに計算させたいのです。
それってどういうコマンドを使えば良いんですかねぇ？
0602１３２人目の素数さん
NGNG<< Graphics`ImplicitPlot`
ImplicitPlot[ {x^y == y^x}, {x, 0.1, 10}, {y, 0.1, 10}]
0603593
NGNG0604１３２人目の素数さん
NGNGしつこいな(w 今まで見た奴は誰も知らなかったから答えてないんだろうよ。
俺も簡単な方法は知らない。
思い付くのは、
1) とりあえず3Dのグラフを書いて、グラフデータの線分
(FullFormすればLine[List[...],List[...]]のようなものでグラフが書かれていることがわかる)
と平面y == 0との交点を求める。
2) y[t] == 0 をFindRootで解く。解がたくさんあるので
FindRootの初期値も十分たくさん取って、重複して出てきた解をのぞく。
最後にそのtに対してx[t], z[t]を計算する。
0605１３２人目の素数さん
NGNGありがとう。色々やってみます。
1)の交点を求めることさえ難しいわけですが。
0606１３２人目の素数さん
NGNGしつこいですね、ごめんなさい
0607593
NGNG皆さん色々アリガトウございました。
0608１３２人目の素数さん
NGNGいくらマセマティカといえどもかれが
答えを出しやすい所まで、お膳立てしてあげるのが
人間の仕事でしょう？
0609１３２人目の素数さん
NGNGso wht
0610１３２人目の素数さん
NGNG0611１３２人目の素数さん
NGNG私が自分で定義した関数をPlot3Dでx軸,y軸の範囲を指定してあげると問題なくグラフが出ます
その点では問題はないのですが、一つだけ行いたいことがあるのでお聞きします
たとえばz軸の値が0〜100で表示されたとします
それを無理やり0〜10の範囲だけに絞って表示することができるのでしょう？
(x,yの範囲の変更を行わないことが条件です)
何がやりたいかといいますと
そのグラフは0〜10が最も多くごく一部が10~90の値を示しているので
そこに0~10の範囲に絞ってメモリ幅をもっと拡大し細かい表示が見たいのです
よろしくお願いします
0612611
NGNG板汚しすいませんでした
0613１３２人目の素数さん
NGNG> いろいろ調べたところ何とかできそうです
> 板汚しすいませんでした
できそうなら方法を書け。書かないから板汚しなんだよ。
PlotRange -> {0, 10}
0614１３２人目の素数さん
NGNGという形で問題を解いていますが、方程式を形つくるある変数を整数から小数
にすると
"行列\!\({\(\[LeftSkeleton] 1 \[RightSkeleton]\)}\)に含まれる悪条件より，結果\
\!\(RowReduce\)には重大な数値的エラーが含まれている可能性があります．"
とでます。解が求まりません。
どうしたらいいのですか、教えてください。整数ならうまくいくのですが小数になった
とたんこのようになります。つまり、行列の中の要素に、小数が入っていることです。
0615614の者
NGNGとは、どういう意味ですか
0616１３２人目の素数さん
NGNG想定しているけど、Solveで浮動小数を用いると、NSolveと同じ結果になって、Solveの本来
意図した結果とは異なることがありますよ、と言う意味だと思う。
浮動小数による数値解を求める時には、最初からNSolveを使った方が早いらしい。
（少なくとも昔は） Mathematicaでは式の中で小数があると、式のデータを全て浮動小数と
して扱うので計算が速くなる。 Mathematica bookのどっかに書いてあった気がする。
0617１３２人目の素数さん
NGNG行列には"条件数"というものがあって、それが「悪い」と
浮動小数点演算では数値的誤差のために逆行列等が正しく計算できない。
詳しくは数値解析の本を。
>>615
多分そこに想定されていた式は別の形(式の組合せ)だったんだろう。
間違った式がわたってるだろうからプログラムを見直せ。
0619614の者
NGNGご返答ありがとうございます。mathmaticaについて少し理解が深まりました
0620１３２人目の素数さん
NGNGおめーら全員逝けよ
0621１３２人目の素数さん
NGNGよかったね。
0622ChaosicSoul ◆/yaJbLAHGw
NGNG0623１３２人目の素数さん
NGNG5.1でもe**leで探せば直ぐに見つかるよ。 鍵メカも前のものが使えるし。
日本語の尻メカはないので、英語しか使えないが。
0624１３２人目の素数さん
NGNG> Re:>620 いいから私にMathematica5を贈呈しろ。
残念ながら4.2しか持っていない。
0625１３２人目の素数さん
NGNG0626１３２人目の素数さん
NGNG0627１３２人目の素数さん
NGNG0628１３２人目の素数さん
NGNGたとえば、あまり一般的でない記法とかを自由に記述できたりしますか？
(どちらがきれいな表示かは問題にしていません。)
0629１３２人目の素数さん
NGNGLatexでしょ
0630とらぬ狸
NGNG0631１３２人目の素数さん
NGNG0632１３２人目の素数さん
NGNG0633１３２人目の素数さん
NGNG俺のマセマティは、ループ構文に弱い。だいたい
バグメッセージを吐くか、言われたとおりやらない。
ちなみにマセマティ４とウィンドー98です。
FORもしくはDoを使って、ループさせ、逐一最後の結果を
あるリストに書き足していくという作戦なんですが
あまり成功したことない。リストに書き足さないとか
FORでくるむと今まで行っていた計算をやらなくなるとか
（マシンリミットをこえたとかメッセージ)。
なぜ？
Print関数は、ロバストなので最近はそれで結果を画面に
吐かせてます。そのあと、こぴぺでテキストエディターにはりつける。
悲しいよね。
0634１３２人目の素数さん
NGNGPolarPlotでグラフが表示されないのはなぜでしょう
In[37]:=
<<Graphics`Graphics`
In[38]:=
PolarPlot[theta,{theta,0,5*Pi}]
Out[38]=
PolarPlot[theta,{theta,0,5 π}]
となってしまうのですが…
0635１３２人目の素数さん
NGNG<<Graphics`Graphics` でパッケージを読み込む前にPolarPlotをうっかり使ってしまったこと。
そうすると、名前が隠されてしまう。 この場合、
Remove[PolarPlot]
として、現在のコンテクストからPolarPlotを消し去るか、（そうすれば、Graphics'Graphics'PolarPlot
にアクセスできる）
Graphics'Graphics'PolarPlot[theta,{theta,0,5 π}]
として、直接目的の関数にアクセスすれば、プロットできるはず。
0636１３２人目の素数さん
NGNGできました。
おそらくパッケージを読み込む前に使ってしまったことだと思います。ありがとうございました
0637１３２人目の素数さん
NGNG0638１３２人目の素数さん
NGNGC言語やJavaなどのプログラミングをやってみれ。Javaなら無料だ。
0639１３２人目の素数さん
NGNG> 悲しいよね。
悲しいね。
> めんどくさいので例は挙げないが、
とか
> あまり成功したことない。リストに書き足さないとか
...
> なぜ？
とか他人に理解できない愚痴書いて
逃避してないで解決しようよ。
>>634
CとかJavaなんかやったらますますForとかDoとかループを使うようになるじゃないか。
Schemeでもやるのがいいだろう。Mathematicaは本来関数型言語なので、Forとか
Doとか、書き換えたりとかいうのは、遅いし読みにくい。
0640１３２人目の素数さん
NGNGアドバイスと励ましのお言葉ありがとう。
C言語はやったことないけどそういう汎用性のある
言語でコーディングした方がいいんだろうね。
マセマティを呼び出すようなコードね。
0641１３２人目の素数さん
NGNGだから、For文でエラーが出たなら、即自分の書き込んだものをチェックするの!
0642１３２人目の素数さん
NGNGきれいに描く方法を教えてくれませんか。
ContourPlot[x^3 - 18 x y + 8y^3, {x, -5.5, 5.5}, {y, -5.5, 5.5},
Contours -> 60,PlotRange -> {-30, 300},
ColorFunction -> Hue]
とやってみたんですが、極地がわからんやら端っこが汚いやらで……。
もう少し何とかできないでしょうか。
いま書いてあるPlotRange や ColorFunction は無視してください。
0643１３２人目の素数さん
NGNGContours -> 15, PlotRange -> {-30, 50},
ColorFunction -> (Hue[Log[1 + #]] &), PlotPoints -> 200]
ではどう？
0644１３２人目の素数さん
NGNGのほうがいいかも。
0645642
NGNGありがとうございました。助かりました。
……指紋みたいなグラフだな。
0646１３２人目の素数さん
NGNGx^2+2y^2-4=0 のグラフを重ねて描きたいのですがどうしたらよいでしょうか。
ContourPlot[(x - 2 y + 2) (4 x - y - 4), {x, -10, 10}, {y, -10, 10},
Contours -> 15, PlotRange -> {-300, 100},
ColorFunction -> (Hue[# Log[1 + #]] &), PlotPoints -> 200]
ParametricPlot[{Sqrt[2] Sin[t], 2 Cos[t]}, {t, 0, 2Pi}]
0647１３２人目の素数さん
NGNGv5.0と比べてどう違いますか？
0648１３２人目の素数さん
NGNGB={{5,6},{7,8}}
とやって、ここから対角ブロック行列 X=[[A 0][0 B]] を定義したいのですが、
どうしたらよいでしょうか？
0649１３２人目の素数さん
NGNGB={{5,6},{7,8}}
Z={{0,0},{0,0}}
X={{A,Z},{Z,B}}
これでどう？
ゼロ行列は関数化されてたと思うけど忘れた。
0650１３２人目の素数さん
NGNG0651１３２人目の素数さん
NGNG0652１３２人目の素数さん
NGNG最初の要素が一番大きい組を取り出す（この場合は｛13，0｝）
方法はあるのでしょうか？
Maxを使えば13となりますが、｛13，0｝という形で取り出したいのです。
もしよろしければ教えてください。
0653１３２人目の素数さん
NGNG変化が大きいんじゃない？
この機能が楽しそうだなと思ったり。
http://www.wolfram.com/products/mathematica/newin51/equationtrekker.html
0654１３２人目の素数さん
NGNG最大数が複数個ある場合を考慮に入れるなら
f[mlist_] := Select[#, Function[x, Part[x, 1] >= Part[#, 1, 1]]] &[Sort[mlist, Part[#2, 1] < Part[#1, 1] &]]
とし、一つだけでならば
f[mlist_]:=Part[Sort[mlist, Part[#2, 1] < Part[#1, 1] &], 1]
仕組みは簡単で、各セルの左側の値で大きい順にソートして最も大きいものを選ぶだけ。
純関数を使っているから見にくいが、やっていること自体は単純。
0655１３２人目の素数さん
NGNG5.1はWinでは x86版、 Linuxでは ｘ86版、x86-64bit版(athlon64)、Itanium版が探せばある。
0656１３２人目の素数さん
NGNGありがとうございます。まだ使い始めて日が浅いので。
こうやればいいんですね。よく考えてみると単純ですね。
0657１３２人目の素数さん
NGNGどなたかご教授願います。
0658１３２人目の素数さん
NGNG0659１３２人目の素数さん
NGNGtest1 = ContourPlot[(x - 2 y + 2) (4 x -
y - 4), {x, -10, 10}, {y, -10, 10}, Contours -> 15, PlotRange -> {-300,
100}, ColorFunction -> (Hue[# Log[1 + #]] &), PlotPoints -> 200,
AspectRatio -> Automatic]
test2 = ParametricPlot[{Sqrt[2] Sin[t], 2 Cos[
t]}, {t, 0, 2Pi}, PlotRange -> {{-10, 10}, {-10, 10}},
AspectRatio -> Automatic]
Show[test1, test2]
0660１３２人目の素数さん
NGNG0661１３２人目の素数さん
NGNGxご教授 o御教示
mathよりjapaneseやれよｗ
0662１３２人目の素数さん
NGNG0663１３２人目の素数さん
NGNGありがとうございます。
Show[test1,test2]で表された図を印刷したいんですが、図だけをコピーすることはできないでしょうか。
枠線が邪魔なんです。
>>661
＿|￣|○
0664１３２人目の素数さん
NGNGを中古で安く譲ってくれる方はいますか？
もし合法でないなら無視して下さい。
合法的に安く入手したいだけですので。
0665１３２人目の素数さん
NGNG仕事で使うのならまずいけど。
0666１３２人目の素数さん
NGNGへんなシト。 英語版ならクレカを持っていれば、ダウソ版が1880ドルで買える。 これが
一番安いと思う。 Mathematicaは業務用ソフトではないし、値段はこんなものか。
日本語版は確か３５万くらいする。
0667１３２人目の素数さん
NGNGgnuです。mathematicaとの互換性は？
0668１３２人目の素数さん
NGNG0669１３２人目の素数さん
NGNGMathematica恐怖症になってしまった。
やる気がおきない。
0670１３２人目の素数さん
NGNG0671１３２人目の素数さん
NGNG0672１３２人目の素数さん
NGNGのグラフの輪っかになってる部分の面積を求めたいんですが、どうしたらよいでしょうか
0673１３２人目の素数さん
NGNGグラフはこれで見れるね。
<< Graphics`ImplicitPlot`
ImplicitPlot[18x y - x^3 - y^3 == 0, {x, -10, 10}]
で、面積は、
18 x y - x^3 - y^3 < 0
の範囲であることは分かりますか？
0674673
NGNGSin[x]やCos[x]を使うんだっけ…
0675１３２人目の素数さん
NGNG>0 の間違いでは？
>>674
version5以降なら、↓で計算できる。でも、式の簡略化に失敗してるなあ。
<< Calculus`Integration`
Integrate[Boole[18 x y - x^3 - y^3 > 0], {x, 0, 10}, {y, 0, 10}]
0676673
NGNG訂正サンクス、やり直さなきゃいけないのは自分の方だった。
N[Integrate[Boole[18 x y - x^3 - y^3 > 0], {x, 0, 10}, {y, 0, 10}] ]
で値がでましたね。
54.
■ このスレッドは過去ログ倉庫に格納されています