- 1 名前：１３２人目の素数さん mailto:sage [04/03/06 09:51]
- _....._{{ 〃
, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,）)
! | ! .,'-.{ ! !|; |`､.}ﾞ!.! |. ! ヽ.
', ', |Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l ! | ﾞ',
ヽ､', l:!Kノ}. f:_.)iﾞi: ﾘ ! l ル
| l!iヽﾞ- ' , .!__:ﾉ ﾞ ,ﾘ l ﾘ'´
. ',|!!､ r‐┐ ｀ ノ'. /,ｲ ／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣
'i!ﾞ､ヽ､ ﾞｰ' _, ィ,:',:''´ ＜ Mathematica に関する話題はここに書いてね！
ﾞ:､ｨ､jヾｰ::: 'iﾍ .ノ',ﾘ. ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
,､- '´ ヽ､ﾞ､ { ｀＞"､
/＼＼ ', } //｀ヽ
★前スレ
〓Mathematica〓
science2.2ch.net/test/read.cgi/math/985023298/
- 552 名前：550 [04/11/23 21:50:19]
- Plot[ArcTan[(Sin[x] - 2.767 * Sin[0.2173 * x]) / (Cos[x] - 2.767 * Cos[0.2173 * x])], {x, 0, 10Pi}]
です．
- 553 名前：クラ [04/11/23 21:52:13]
-
押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない押部谷中学校の奴 皆死んでしまえ
私はロバなんかじゃない
- 554 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/23 21:55:43]
- Re:>552
Plot[ArcTan[1/x],{x,-1,1}]
をやってみてくれ。
これでジャンプの原因が分かる。
- 555 名前：１３２人目の素数さん [04/11/23 21:59:24]
- 484 ： ：04/11/23 19:11 HOST:YahooBB219174040245.bbtec.net<8080>
152 ：依頼 ：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 ◆〜〜〜〜〜」における共同荒らしが２ヶ月ほど絶え間なく続いていて数学の議論ができない状態です。
このキーワードでレスを摘出していただければわかります。
science3.2ch.net/test/read.cgi/math/1097495449/
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>
- 556 名前：550 [04/11/23 22:06:18]
- >554
0で-pi/2からpi/2にジャンプしています．関係ないとおもって書きませんでしたが，versionは3です．
- 557 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/23 22:11:47]
- Re:>556 念のために説明すると、lim_{t→∞}Arctan(t)=π/2,lim_{t→∞}Arctan(t)=-π/2.
- 558 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/23 22:12:25]
- ↑lim_{t→-∞}Arctan(t)=-π/2.
- 559 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/23 22:13:30]
- とにかくArctan(+∞)=π/2,Arctan(-∞)=-π/2.
[>558]の上矢印は数学記号ではない。
- 560 名前：550 [04/11/23 22:14:52]
- >554
arctanは多価ですが，(-pi/2, pi/2)に限定せず，とにかくグラフが連続になるようにスマートにやってもらいたいってことです．
- 561 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/23 22:20:33]
- Re:>560 そういう意味か。
とりあえず、先ずはCos[x] - 2.767 * Cos[0.2173 * x]の零点および、前後の符号を調べて、場合分けするといいだろう。
符号が-から+に変わるところでPiを引いて、+から-に変わるところでPiを足せば連続になる。
- 562 名前：”削除”板荒らされる [04/11/23 22:29:24]
- qb5.2ch.net/test/read.cgi/saku/1033142451/l50
qb5.2ch.net/test/read.cgi/saku/1027349232/l50
- 563 名前：550 [04/11/23 22:30:07]
- >561
なるほど．どうもありがとうございます．ところで，それはMathematicaのコマンドでできるんですか?
- 564 名前：563 [04/11/23 22:36:17]
- f(x) = cos x - 2.767 * cos(0.2173 * x)
g(x) = f(x) if f(x) < 0
g(x) = f(x) + pi if f(x) > 0
みたいなのが簡単に書けるのか?ってことです．
- 565 名前：564 [04/11/23 22:37:20]
- 563のようなg(x)をどうMathematicaでやるのかってことです．
- 566 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/23 22:42:40]
- Re:>563-565 私は何年か前にMathematicaの使用をやめたのでよく分からないが、If関数で出来るだろう。
- 567 名前：564 [04/11/23 22:58:16]
- >566
ありがとうございます．やれそうですね．
- 568 名前：１３２人目の素数さん mailto:sage [04/11/24 03:02:19]
- >>552
prev = 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]
- 569 名前：１３２人目の素数さん [04/11/24 19:27:53]
- mathematica5.1記念あげ
- 570 名前：１３２人目の素数さん [04/11/24 19:51:34]
- >>565
module
都会魚まじないを使って関数を定義すれば
よし。
- 571 名前：１３２人目の素数さん [04/11/25 15:33:18]
- ∫ g[x] f[z] dz があらかじめ定義されているとして
∫f[z] dz = Fを定義して、∫ g[x] f[z] dzをF g[x]と
変形するような計算はmathematicaではどうやってやるのでしょうか？
積分変数でない変数しか含まない関数を外に追い出す命令があればあとは
replaceで簡単なのですが・・・
- 572 名前：１３２人目の素数さん mailto:sage [04/11/25 16:28:15]
- ずっと考えてやっと閃きました
こうすればいいようです。
Unprotected[Integrate]
Integrate[f[x],{x,-inf,inf}]:=F
これで、この積分は自動的に評価されます
- 573 名前：１３２人目の素数さん mailto:sage [04/11/25 16:29:58]
- >>571
> 積分変数でない変数しか含まない関数を外に追い出す
簡単に作れるよ。
- 574 名前：１３２人目の素数さん mailto:sage [04/11/25 16:32:09]
- >>572
って自動で追い出されるじゃないか…
- 575 名前：１３２人目の素数さん mailto:sage [04/11/25 17:09:27]
- >>573 >>574
なにか設定が辺なのか追い出されない。
どういうことなんでしょうねぇ・・・
- 576 名前：卒業を助けてくださいm(_ _)m [04/11/25 18:15:18]
- 三次元上に幾つかの球を配置し、それをあらゆる平面で切り、断面積もしくは断面（円）
の半径を求めたいのですがどうしたらいいでしょうか？
これができないと卒業があやうくなります。ちなみに球の座標は
(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)です。
よろしくお願いします。
- 577 名前：１３２人目の素数さん [04/11/26 15:36:22]
- >>569
5.1でたのか、知らんかった。
アップグレードのお知らせとか送ってくれないのかな。
ところで5.0からのアップグレードはいくらくらいぼったくられるのだろう？
- 578 名前：１３２人目の素数さん mailto:sage [04/11/26 15:37:52]
- >>576
留年しろ
- 579 名前：１３２人目の素数さん [04/11/27 14:36:41]
- Mathematica と LaTex ではどちらが数式の記述能力が上ですか？
たとえば、あまり一般的でない記法とかを自由に記述できたりしますか？
(どちらがきれいな表示かは問題にしていません。)
- 580 名前：伊丹公理 [04/11/27 15:42:24]
- 初等的なものなら、Mathematica
編集ならLaTeX
自由自在、縦横無尽に使うなら TeX
- 581 名前：１３２人目の素数さん [04/11/27 15:43:07]
- >>579
普通にそれぞれの用途を考えたらわかると思うが？
- 582 名前：１３２人目の素数さん mailto:sage [04/11/27 15:46:04]
- >>580 >>581 ありがとう。
- 583 名前：１３２人目の素数さん [04/11/30 20:07:54]
-
- 584 名前：１３２人目の素数さん [04/12/02 15:37:05]
- Mathematica5なんですが、EmacsのParen Match Highlightingのような機能
(対応する括弧を強調表示する機能)は無いでしょうか?
- 585 名前：１３２人目の素数さん [04/12/02 15:44:17]
- Mathnoteいい
homepage2.nifty.com/autumn-soft/
- 586 名前：１３２人目の素数さん mailto:sage [04/12/03 00:17:14]
- >>584
Edit -> Check Balance [Shitf Ctrl B]
- 587 名前：１３２人目の素数さん mailto:sage [04/12/03 10:58:08]
- >>586
ありがとう。
自動的に強調表示してくれる機能は無いのね。残念。
- 588 名前：１３２人目の素数さん mailto:sage [04/12/03 12:21:52]
- >>586
面白いね。
- 589 名前：１３２人目の素数さん [04/12/03 21:26:21]
- Mathematica上での置換について質問します。
Mathematicaで置換するやり方はわかるのですが、
テキストファイルなどと違って置換する速度がかなり遅いです。
数個程度の置換ならば良いのですが、置換する数が数万個単位になるような
長い数式だと普通の置換のやり方ではさすがに時間がかかりすぎてしまいます。
今回困っているパターンについて
例えば数式で ax+by に対して"ax+by"のように「"」がついています
これではax+by自体を文字列とみなしてしまい中の計算は行ってくれません
Mathematica上で「"」→「」の置換を行えば問題は解決するのですが、これが数万個単位であるのです
特別なテキスト形式で保存しても文字列と数列では単純に「"」のあるなしではないので
「"」だけを削除すればよいという問題でもありません
よろしければ
Mathematica上で置換をテキストレベルに早く行うことができるやり方
もしくはMathematicaの表記どおりに「"」を削除を行うだけで解決するようなやり方
知っている方がいればお願いします
ツールの有無、シェアの有無にはこだわりません
- 590 名前：１３２人目の素数さん mailto:sage [04/12/04 02:17:27]
- >>589
数万個どうやって作ったわけ？
手入力してしまったとかいうのでなければ、そのデータ入力のところからやりなおすのがいい。
参考>>402,>>404。
試してみるとバグだと言いたくなる遅さだね。フロントエンドを表示せずにはNotebook操作系の関数も使えないのが痛い。
- 591 名前：１３２人目の素数さん [04/12/06 02:08:41]
- \!\(For[i = 1,
i <= 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}]\)
のどこが悪いですか
- 592 名前：１３２人目の素数さん mailto:sage [04/12/06 04:46:54]
- >>591
DM2 = Array[DM, {C\_A + 1, C\_A + 1}]\)
にしたらわかる。
- 593 名前：１３２人目の素数さん mailto:sage [04/12/07 16:19:29]
- solution=
NDSolve[
{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方向無限遠点から見たような図になってしまいました。
- 594 名前：１３２人目の素数さん [04/12/07 19:14:48]
- v=d(a+p/q)^(-c)*y
このような非線形の等式をpについて解いて、
p=の格好にしたいのですが、
こういう式変形をさせるにはどうしたらいいのでしょうか？
教えて下さい。
- 595 名前：550 [04/12/07 20:17:28]
- >568
遅くなりましたが，ありがとうございます．コードの意味も大体分かりました．
- 596 名前：１３２人目の素数さん [04/12/07 20:35:23]
- Plot3Dでグラフを書くと、外側に直方体の枠が出るんですが
あの枠を出さない方法を教えてつかあさい。
- 597 名前：１３２人目の素数さん [04/12/07 20:50:25]
- x^y=y^xのグラフを表示させたいのですが、
Mathematicaで陰関数を表示する方法はありますか？
GCalcやfunctionViewなどでできるんですが、
Mathematicaでもできるのかなと思いまして。
- 598 名前：１３２人目の素数さん mailto:sage [04/12/07 21:06:31]
- >>594
First 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
- 599 名前：伊丹公理 [04/12/07 21:09:06]
- >>597
(log x)/x の重複値を調べる
- 600 名前：１３２人目の素数さん mailto:sage [04/12/07 21:17:42]
- >>593
お願いします。
x[t],y[t],z[t]の式出してそこにy[t]の値を代入、とかできるのかな？
何にせよ今mathematica使えないからなぁ…
- 601 名前：１３２人目の素数さん [04/12/07 22:23:09]
- >598
返信ありがとうございます。
でも、もちろん、私も手では計算できるんですけど、
もっと実はもっと複雑な非線形の等式で似たような事を
したいので、それをmathematicaに計算させたいのです。
それってどういうコマンドを使えば良いんですかねぇ？
- 602 名前：１３２人目の素数さん mailto:sage [04/12/07 22:29:26]
- >>597
<< Graphics`ImplicitPlot`
ImplicitPlot[ {x^y == y^x}, {x, 0.1, 10}, {y, 0.1, 10}]
- 603 名前：593 mailto:sage [04/12/08 00:07:41]
- お答えお待ちしてます… orz dogeza
- 604 名前：１３２人目の素数さん mailto:sage [04/12/08 00:24:22]
- >>603
しつこいな(w 今まで見た奴は誰も知らなかったから答えてないんだろうよ。
俺も簡単な方法は知らない。
思い付くのは、
1) とりあえず3Dのグラフを書いて、グラフデータの線分
(FullFormすればLine[List[...],List[...]]のようなものでグラフが書かれていることがわかる)
と平面y == 0との交点を求める。
2) y[t] == 0 をFindRootで解く。解がたくさんあるので
FindRootの初期値も十分たくさん取って、重複して出てきた解をのぞく。
最後にそのtに対してx[t], z[t]を計算する。
- 605 名前：１３２人目の素数さん mailto:sage [04/12/08 06:54:01]
- >>604
ありがとう。色々やってみます。
1)の交点を求めることさえ難しいわけですが。
- 606 名前：１３２人目の素数さん mailto:sage [04/12/08 07:13:42]
- もう少し詳しく…自分でもやりますが。
しつこいですね、ごめんなさい
- 607 名前：593 mailto:sage [04/12/08 15:46:38]
- y[t]をz軸に設定して、PlotRangeで-0.001,0.000と指定して無理やり断面を求めました。
皆さん色々アリガトウございました。
- 608 名前：１３２人目の素数さん [04/12/08 16:14:09]
- >>593
いくらマセマティカといえどもかれが
答えを出しやすい所まで、お膳立てしてあげるのが
人間の仕事でしょう？
- 609 名前：１３２人目の素数さん mailto:sage [04/12/08 18:47:48]
- >>608
so wht
- 610 名前：１３２人目の素数さん mailto:sage [04/12/08 18:52:10]
- Mathematicaで萌え萌えはにゃ〜んヽ(´ー`)ﾉ
- 611 名前：１３２人目の素数さん [04/12/09 00:25:59]
- グラフ表示で質問です
私が自分で定義した関数をPlot3Dでx軸,y軸の範囲を指定してあげると問題なくグラフが出ます
その点では問題はないのですが、一つだけ行いたいことがあるのでお聞きします
たとえばz軸の値が0〜100で表示されたとします
それを無理やり0〜10の範囲だけに絞って表示することができるのでしょう？
(x,yの範囲の変更を行わないことが条件です)
何がやりたいかといいますと
そのグラフは0〜10が最も多くごく一部が10~90の値を示しているので
そこに0~10の範囲に絞ってメモリ幅をもっと拡大し細かい表示が見たいのです
よろしくお願いします
- 612 名前：611 mailto:sage [04/12/09 00:42:14]
- いろいろ調べたところ何とかできそうです
板汚しすいませんでした
- 613 名前：１３２人目の素数さん mailto:sage [04/12/09 03:07:44]
- >>612
> いろいろ調べたところ何とかできそうです
> 板汚しすいませんでした
できそうなら方法を書け。書かないから板汚しなんだよ。
PlotRange -> {0, 10}
- 614 名前：１３２人目の素数さん [04/12/09 18:36:07]
- Solve関数を使って多項方程式でベクトル*行列＝ベクトル
という形で問題を解いていますが、方程式を形つくるある変数を整数から小数
にすると
"行列\!\({\(\[LeftSkeleton] 1 \[RightSkeleton]\)}\)に含まれる悪条件より，結果\
\!\(RowReduce\)には重大な数値的エラーが含まれている可能性があります．"
とでます。解が求まりません。
どうしたらいいのですか、教えてください。整数ならうまくいくのですが小数になった
とたんこのようになります。つまり、行列の中の要素に、小数が入っていることです。
- 615 名前：614の者 [04/12/09 18:47:07]
- あと、”原始的ではない式が想定されます”
とは、どういう意味ですか
- 616 名前：１３２人目の素数さん mailto:sage [04/12/09 19:29:32]
- 基本的にSolveは多倍長数や記号を用いた厳密な解、NSolveは浮動小数を用いた近似解を
想定しているけど、Solveで浮動小数を用いると、NSolveと同じ結果になって、Solveの本来
意図した結果とは異なることがありますよ、と言う意味だと思う。
浮動小数による数値解を求める時には、最初からNSolveを使った方が早いらしい。
（少なくとも昔は） Mathematicaでは式の中で小数があると、式のデータを全て浮動小数と
して扱うので計算が速くなる。 Mathematica bookのどっかに書いてあった気がする。
- 617 名前：１３２人目の素数さん mailto:sage [04/12/09 21:08:06]
- >>614
行列には"条件数"というものがあって、それが「悪い」と
浮動小数点演算では数値的誤差のために逆行列等が正しく計算できない。
詳しくは数値解析の本を。
>>615
多分そこに想定されていた式は別の形(式の組合せ)だったんだろう。
間違った式がわたってるだろうからプログラムを見直せ。
- 618 名前：617 mailto:sage [04/12/09 21:09:29]
- つまり>>616は間違ってる。
- 619 名前：614の者 [04/12/09 23:56:13]
- 分数で表すことで解決しました。
ご返答ありがとうございます。mathmaticaについて少し理解が深まりました
- 620 名前：１３２人目の素数さん [04/12/10 07:51:43]
-
おめーら全員逝けよ
- 621 名前：１３２人目の素数さん [04/12/11 16:51:47]
- >>619
よかったね。
- 622 名前：ChaosicSoul ◆/yaJbLAHGw [04/12/11 17:42:19]
- Re:>620 いいから私にMathematica5を贈呈しろ。
- 623 名前：１３２人目の素数さん mailto:sage [04/12/11 17:55:38]
- >>622
5.1でもe**leで探せば直ぐに見つかるよ。 鍵メカも前のものが使えるし。
日本語の尻メカはないので、英語しか使えないが。
- 624 名前：１３２人目の素数さん [04/12/12 01:46:38]
- >>622
> Re:>620 いいから私にMathematica5を贈呈しろ。
残念ながら4.2しか持っていない。
- 625 名前：１３２人目の素数さん mailto:sage [04/12/14 07:52:33]
- だれか、Publiconって使ったことのある人います？ いれば、レポお願いします。
- 626 名前：１３２人目の素数さん [04/12/16 00:11:15]
- age
- 627 名前：１３２人目の素数さん [04/12/22 22:27:01]
- 262
- 628 名前：１３２人目の素数さん mailto: [04/12/23 19:02:33]
- Mathematica と LaTex ではどちらが数式の記述能力が上ですか？
たとえば、あまり一般的でない記法とかを自由に記述できたりしますか？
(どちらがきれいな表示かは問題にしていません。)
- 629 名前：１３２人目の素数さん [04/12/23 21:59:47]
- >>628
Latexでしょ
- 630 名前：とらぬ狸 mailto:sage [04/12/23 22:01:42]
- むしろLaTeXでは？
- 631 名前：１３２人目の素数さん mailto:sage [04/12/23 22:12:43]
- いや、LaTeX だろ。
- 632 名前：１３２人目の素数さん mailto:sage [04/12/23 22:55:03]
- Publicon使って見てくれよ。
- 633 名前：１３２人目の素数さん [04/12/23 23:19:55]
- めんどくさいので例は挙げないが、
俺のマセマティは、ループ構文に弱い。だいたい
バグメッセージを吐くか、言われたとおりやらない。
ちなみにマセマティ４とウィンドー98です。
FORもしくはDoを使って、ループさせ、逐一最後の結果を
あるリストに書き足していくという作戦なんですが
あまり成功したことない。リストに書き足さないとか
FORでくるむと今まで行っていた計算をやらなくなるとか
（マシンリミットをこえたとかメッセージ)。
なぜ？
Print関数は、ロバストなので最近はそれで結果を画面に
吐かせてます。そのあと、こぴぺでテキストエディターにはりつける。
悲しいよね。
- 634 名前：１３２人目の素数さん mailto:sage [04/12/24 02:25:40]
- 5.1なんだけれど、
PolarPlotでグラフが表示されないのはなぜでしょう
In[37]:=
<<Graphics`Graphics`
In[38]:=
PolarPlot[theta,{theta,0,5*Pi}]
Out[38]=
PolarPlot[theta,{theta,0,5 π}]
となってしまうのですが…
- 635 名前：１３２人目の素数さん mailto:sage [04/12/24 03:39:07]
- 当方、５．１ですが、ちゃんと貝殻ができたよ。 理由として考えられるのは
<<Graphics`Graphics` でパッケージを読み込む前にPolarPlotをうっかり使ってしまったこと。
そうすると、名前が隠されてしまう。 この場合、
Remove[PolarPlot]
として、現在のコンテクストからPolarPlotを消し去るか、（そうすれば、Graphics'Graphics'PolarPlot
にアクセスできる）
Graphics'Graphics'PolarPlot[theta,{theta,0,5 π}]
として、直接目的の関数にアクセスすれば、プロットできるはず。
- 636 名前：１３２人目の素数さん mailto:sage [04/12/24 04:25:05]
- >>635
できました。
おそらくパッケージを読み込む前に使ってしまったことだと思います。ありがとうございました
- 637 名前：１３２人目の素数さん [04/12/24 05:29:44]
- age
- 638 名前：１３２人目の素数さん mailto:sage [04/12/24 07:43:59]
- >>633
C言語やJavaなどのプログラミングをやってみれ。Javaなら無料だ。
- 639 名前：１３２人目の素数さん mailto:sage [04/12/24 09:44:52]
- >>633
> 悲しいよね。
悲しいね。
> めんどくさいので例は挙げないが、
とか
> あまり成功したことない。リストに書き足さないとか
...
> なぜ？
とか他人に理解できない愚痴書いて
逃避してないで解決しようよ。
>>634
CとかJavaなんかやったらますますForとかDoとかループを使うようになるじゃないか。
Schemeでもやるのがいいだろう。Mathematicaは本来関数型言語なので、Forとか
Doとか、書き換えたりとかいうのは、遅いし読みにくい。
- 640 名前：１３２人目の素数さん [04/12/24 10:32:26]
- >>639
アドバイスと励ましのお言葉ありがとう。
C言語はやったことないけどそういう汎用性のある
言語でコーディングした方がいいんだろうね。
マセマティを呼び出すようなコードね。
- 641 名前：１３２人目の素数さん mailto:sage [04/12/24 11:07:30]
- コンピュータ言語は、人間側が機械に合わせてあげるのがポイントなんだよ。
だから、For文でエラーが出たなら、即自分の書き込んだものをチェックするの!
- 642 名前：１３２人目の素数さん [04/12/24 12:57:43]
- z=x^3 - 18 x y + 8y^3 グラフを等高線プロットで描きたいんです。
きれいに描く方法を教えてくれませんか。
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 は無視してください。
- 643 名前：１３２人目の素数さん mailto:sage [04/12/24 17:59:00]
- ContourPlot[x^3 - 18 x y + 8y^3, {x, -5.5, 5.5}, {y, -5.5, 5.5},
Contours -> 15, PlotRange -> {-30, 50},
ColorFunction -> (Hue[Log[1 + #]] &), PlotPoints -> 200]
ではどう？
- 644 名前：１３２人目の素数さん mailto:sage [04/12/24 21:02:11]
- ColorFunction -> (Hue[# Log[2]] &)
のほうがいいかも。
- 645 名前：642 [04/12/25 09:42:34]
- >643 >644
ありがとうございました。助かりました。
……指紋みたいなグラフだな。
- 646 名前：１３２人目の素数さん [04/12/25 10:30:30]
- z=(x-2y+2)(4x-y-4) の等高線プロットの上に
x^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}]
- 647 名前：１３２人目の素数さん mailto:sage [04/12/25 10:37:31]
- v5.1がリリースされてから約1ヶ月が経ちましたが、
v5.0と比べてどう違いますか？
- 648 名前：１３２人目の素数さん [04/12/25 13:20:09]
- A={{1,2},{3,4}}
B={{5,6},{7,8}}
とやって、ここから対角ブロック行列 X=[[A 0][0 B]] を定義したいのですが、
どうしたらよいでしょうか？
- 649 名前：１３２人目の素数さん mailto:sage [04/12/25 13:26:52]
- A={{1,2},{3,4}}
B={{5,6},{7,8}}
Z={{0,0},{0,0}}
X={{A,Z},{Z,B}}
これでどう？
ゼロ行列は関数化されてたと思うけど忘れた。
- 650 名前：１３２人目の素数さん mailto:sage [04/12/25 15:39:05]
- もう5.1にしているやつ結構いるんだな。
- 651 名前：１３２人目の素数さん mailto:sage [04/12/25 18:17:38]
- mathematica 5.0 とか 5.1 は割れ物として見つかる？
- 652 名前：１３２人目の素数さん mailto:sage [04/12/25 20:51:19]
- 例えば、{{10,3},{11,2},{13,0}}のようなリストから
最初の要素が一番大きい組を取り出す（この場合は｛13，0｝）
方法はあるのでしょうか？
Maxを使えば13となりますが、｛13，0｝という形で取り出したいのです。
もしよろしければ教えてください。
|
 |
