〓 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/
0552550
NGNGです．
0553クラ
NGNG押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない
0554ChaosicSoul ◆/yaJbLAHGw
NGNGPlot[ArcTan[1/x],{x,-1,1}]
をやってみてくれ。
これでジャンプの原因が分かる。
0555１３２人目の素数さん
NGNG152 ：依頼 ：04/10/11 15:38:30 HOST:33.93.215.220.ap.yournet.ne.jp
執拗なまでのコピペ荒しです。「うんち食いたい」や、某コテハンのアドレスを各スレにコピペしながら回っているようです。
これでもまだ1/5ぐらいの量です。よければ削除お願いします。永久アク禁してもらいたいぐらいですが。
153 ：依頼２ ：04/10/20 23:24:48 HOST:14.91.215.220.ap.yournet.ne.jp
名前「********@yahoo.co.jp」（名前がメールアドレスなので一応隠しました）と、
名前「LettersOfLiberty ◆〜〜〜〜〜」（〜〜〜はいろいろと）、
名前「FeaturesOfTheGod ◆〜〜〜〜〜」における共同荒らしが２ヶ月ほど絶え間なく続いていて数学の議論ができない状態です。
このキーワードでレスを摘出していただければわかります。
http://science3.2ch.net/test/read.cgi/math/1097495449/
http://science3.2ch.net/test/read.cgi/math/1095390340/
この２スレを見ていただければ、続いている荒らしについての議論がされています。
挙げた例はほんのわずかな例です。上から順にスレを開けばほとんどのスレが荒されているのがわかります。
いくつかすでにレスが削除されている様子ですが、それは荒らしレスの1/100ほどです。
尋常じゃないです、どうにかしていただきたい。
157 ：∂ ：04/11/20 05:40:32 HOST:65.98.66.20
163 ：π ：04/11/20 21:11 HOST:tetkyo024225.tkyo.te.ftth2.ppp.infoweb.ne.jp<80><8080>
171 ： ：04/11/23 16:10 HOST:glass.ipe.tsukuba.ac.jp<80><8080><3128><8000><1080>
0556550
NGNG0で-pi/2からpi/2にジャンプしています．関係ないとおもって書きませんでしたが，versionは3です．
0557ChaosicSoul ◆/yaJbLAHGw
NGNG0558ChaosicSoul ◆/yaJbLAHGw
NGNG0559ChaosicSoul ◆/yaJbLAHGw
NGNG[>558]の上矢印は数学記号ではない。
0560550
NGNGarctanは多価ですが，(-pi/2, pi/2)に限定せず，とにかくグラフが連続になるようにスマートにやってもらいたいってことです．
0561ChaosicSoul ◆/yaJbLAHGw
NGNGとりあえず、先ずはCos[x] - 2.767 * Cos[0.2173 * x]の零点および、前後の符号を調べて、場合分けするといいだろう。
符号が-から+に変わるところでPiを引いて、+から-に変わるところでPiを足せば連続になる。
0562”削除”板荒らされる
NGNGhttp://qb5.2ch.net/test/read.cgi/saku/1027349232/l50
0563550
NGNGなるほど．どうもありがとうございます．ところで，それはMathematicaのコマンドでできるんですか?
0564563
NGNGg(x) = f(x) if f(x) < 0
g(x) = f(x) + pi if f(x) > 0
みたいなのが簡単に書けるのか?ってことです．
0565564
NGNG0566ChaosicSoul ◆/yaJbLAHGw
NGNG0567564
NGNGありがとうございます．やれそうですね．
0568１３２人目の素数さん
NGNGprev = Null;
YourArcTan[x_] := If[prev === Null, prev = ArcTan[x],
prev = ArcTan[x] + Pi * IntegerPart[(prev - ArcTan[x] + 0.5) / Pi]]
ListPlot[Table[
YourArcTan[(Sin[x] - 2.767*Sin[0.2173*x])/(Cos[x] -
2.767*Cos[0.2173*x])], {x, 0, 10 Pi, 0.1}], PlotJoined -> True]
0569１３２人目の素数さん
NGNG0570１３２人目の素数さん
NGNGmodule
都会魚まじないを使って関数を定義すれば
よし。
0571１３２人目の素数さん
NGNG∫f[z] dz = Fを定義して、∫ g[x] f[z] dzをF g[x]と
変形するような計算はmathematicaではどうやってやるのでしょうか？
積分変数でない変数しか含まない関数を外に追い出す命令があればあとは
replaceで簡単なのですが・・・
0572１３２人目の素数さん
NGNGこうすればいいようです。
Unprotected[Integrate]
Integrate[f[x],{x,-inf,inf}]:=F
これで、この積分は自動的に評価されます
0573１３２人目の素数さん
NGNG> 積分変数でない変数しか含まない関数を外に追い出す
簡単に作れるよ。
0574１３２人目の素数さん
NGNGって自動で追い出されるじゃないか…
0575１３２人目の素数さん
NGNGなにか設定が辺なのか追い出されない。
どういうことなんでしょうねぇ・・・
0576卒業を助けてくださいm(_ _)m
NGNGの半径を求めたいのですがどうしたらいいでしょうか？
これができないと卒業があやうくなります。ちなみに球の座標は
(0.0.0)(2.0.0)(2.2.0)(0.2√2.0)(0.0.2√2)(2√2.0.2√2)(2√2.2√2.2√2)(0.2√2.2√2)
(√2.0.√2)(2√2.√2.√2)(√2.2√2.√2)(0.√2.√2)(√2.√2.0)(√2.√2.2√2)です。
よろしくお願いします。
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｝という形で取り出したいのです。
もしよろしければ教えてください。
■ このスレッドは過去ログ倉庫に格納されています