〓 Mathematica 2 〓

[表示 : 全て最新50 1-99 101- 201- 301- 401- 501- 601- 701- 801- 901- 1001- 2chのread.cgiへ]
Update time : 12/19 02:15 / Filesize : 256 KB / Number-of Response : 1002
[このスレッドの書き込みを削除する]
[＋板最近立ったスレ＆熱いスレ一覧 : ＋板最近立ったスレ／記者別一覧] [類似スレッド一覧]


↑キャッシュ検索、類似スレ動作を修正しました、ご迷惑をお掛けしました

〓 Mathematica 2 〓

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/
511 名前：１３２人目の素数さん mailto:sage [04/11/13 09:00:23]: >>510
意味は分からないが…

単純に言うと
a = 1 + 3 (のような)
f[a] = ...
と定義したらaと入力して1 + 3が返り、
f[a]と入力してf[1] + f[3]みたいな表示が返ればいいのか？
a = HoldForm[1 + 3]
f[a:HoldForm[x + y]] = HoldForm[f[x] + f[y]]

いづれにせよ
計算されてなければa, bからc, d,...やPlus, Minusをとり出すことはできる。
512 名前：１３２人目の素数さん mailto:sage [04/11/13 09:10:47]: MinusじゃなくてTimes[-1,..]という形か。
513 名前：510 [04/11/13 20:10:41]: >>511
ありがとうございます
一つ目の質問に関しては「HoldForm」が私が思っていた理想とするものでした
大変参考になりました

これを元に2つ目の
f[a,b] = f[{c,d},{g,h}] - f[{c,d},{i,j}] + f[{e,f},{g,h}] - f[{e,f},{i,j}]
についても考えたいと思います。

これは要は
Outer[ F,{a,b},{c,d} ] = {F(a,b),F(a,c),F(b,c),F(b,d)}
というものがあると思いますが、これと似たような機能で
　h[ F,{ list1,list2 },{ list3,list4 } ]
　={F(list1,list3)},F(list1,list4),F(list2,list3),F(list2,list4)]
のようにリストの中にあるリストを分配してくれるとうなものを探していました。
引き続き調べてみようと思います。
また何かありましたらよろしくお願いいたします。
ありがとうございました
514 名前：１３２人目の素数さん mailto:sage [04/11/13 21:10:06]: >>513
>Outer[ F,{a,b},{c,d} ] = {F(a,b),F(a,c),F(b,c),F(b,d)}
ではなくて、
Outer[F, {a, b}, {c, d}]
{{F[a, c], F[a, d]}, {F[b, c], F[b, d]}}

同様に、リストもOuterでよい。必要ならFlattenする。
Outer[F, {{1, 2}, {3, 4}}, {{5, 6}, {7, 8}}, 1]
{{F[{1, 2}, {5, 6}], F[{1, 2}, {7, 8}]}, {F[{3, 4}, {5, 6}],
F[{3, 4}, {7, 8}]}}
515 名前：510 [04/11/13 21:22:22]: >>514
御指摘とアドバイスをありがとうございます
　Outer[F,list1,list2, ... ,listn, m ]
のmの部分によってどのリストに適用できるかどうかを指定できる
という性質を見逃していました。
おかげさまで一つの悩みが解決し、先に進めそうです。
ありがとうございました。
516 名前：１３２人目の素数さん [04/11/13 21:40:46]: HoldForm
フレッドホルム
517 名前：１３２人目の素数さん [04/11/15 00:18:14]: ｋの変数としての関数ｇ　は、

ｇ（ｋ）＝係数1（ｋ）×ｇ（ｋ-オフセット）＋係数２（ｋ）×ｇ（ｋ+オフセット）

の形に従う。これを逐次近似でとけ
っていわれても。といて
ｇ　vs　ｋ　曲線を導け！といわれても、難しいな。
なんか初期試行関数のあたりをつける（Mathで）いい方法ない？
518 名前：１３２人目の素数さん [04/11/15 11:28:14]: Integrate[E^0.5/(1+Exp[E-eta]),{E,0,infty}]
で定義される特殊関数の
名前は？（変数はetaです）

教えてください。F1／2とかいうシンボルが付加されてます
519 名前：１３２人目の素数さん mailto:sage [04/11/15 14:16:01]: >>518
正しい式を書き込んでね。
520 名前：１３２人目の素数さん mailto:sage [04/11/15 14:35:18]: F_1/2[eta_]:=
(2/Sqrt[Pi])*Integrate[E^0.5/(1+Exp[E-eta]),{E,0,infty}]
521 名前：１３２人目の素数さん mailto:sage [04/11/15 14:40:46]: 分母のところが　１＋指数関数

なので、γ関数ではなさそう。
522 名前：１３２人目の素数さん mailto:sage [04/11/15 14:50:29]: www.itl.nist.gov/div898/software/dataplot/refman2/auxillar/fermdira.htm

これだな。しかし、マセマティカには組み込まれてないようだ。
LerchPhiとかいうのとガンマ関数を
両方使って定義しても良いけど　そんなんするくらいなら
自分で定義した方がすっきり
523 名前：１３２人目の素数さん [04/11/15 16:08:41]: 二項演算ｍを言語とし、次を公理とする理論のモデルで、ドメインが
｛１，２，３｝であり、e＝１となるものをすべてあげよ。

●∀x[m(e,x)=x]
●∀x[m(x,e)=x]
●∀x∃y[x=m(y,y)]

誰か教えてくらさい。
524 名前：１３２人目の素数さん [04/11/15 17:18:19]: 一つお尋ねします

F[ {a,b}+{c,d} , {e,f}-{g,h} ]
を
F[{a,b},{e,f}] - F[{a,b},{g,h}] + F[{c,d},{e,f}] - F[{c,d},{g,h}]
に変えるようなものを作りたいと考えています

Distribute[ F[a + b, c + d] ]
↓
F[a, c] + F[a, d] + F[b, c] + F[b, d]
というのは見つけたのですが目的のものができません
どうも関数の中の値がリストになるとわからなくなってしまいます。
どうかご指導よろしくお願いいたします
525 名前：１３２人目の素数さん mailto:sage [04/11/15 18:18:25]: >>524
元ネタがわからんからこうではないんかもしれんが、

SetAttributes[MyDistribute, HoldAll]
MyDistribute[f_[x___]] := Module[{g}, SetAttributes[g, HoldAll];
Distribute[g[x]] /. g -> f]
F[a_, {- b_, -c_}] := - F[a, {b, c}]

MyDistribute[F[{a, b} + {c, d}, {e, f} - {g, h}]]
526 名前：１３２人目の素数さん mailto:sage [04/11/15 18:24:13]: このほうがいいかもしれん。
SetAttributes[MyDistribute, HoldAll]
MyDistribute[f_[x___]] := Module[{g}, SetAttributes[g, HoldAll];
g[a_, -b_] = - g[a, b];
Distribute[g[x]] /. g -> f]
MyDistribute[F[{a, b} + {c, d}, {e, f} - {g, h}]]
527 名前：１３２人目の素数さん [04/11/15 19:11:31]: 二変数の関数Fujをかんがえる。ここで変数は　iとｊ。

要素をFuj_{i,j}とする正方行列をつくりdeterminentを得たい
とき。どうすればいい？

Det[Flatten[{Table[N[Fuj[i, j]], {i, 1, 4, 1}, {j, 2, 5, 1}]}, 1]]？？
528 名前：１３２人目の素数さん mailto:sage [04/11/15 19:37:03]: >>527
Det[Table[Fuj[i, j], {i, 1, 4}, {j, 1, 4}]]
529 名前：１３２人目の素数さん [04/11/15 20:13:18]: >>528

これだとサブグループを作ってしまって、
正方行列にならないのでは？

list[i, j]=Table[Table[Fuj[i, j], {i, m1, m100}], {j, n1, n100}]

で行列を作っておいて、行列式へ。
そのほうがいくない？
530 名前：１３２人目の素数さん mailto:sage [04/11/15 20:58:56]: >>529
はぁ？それ間違ってるし。
普通Fuj[i,j]ならiが行インデックスだろ。Array[Fuj,{4,4}]と比べてみろ。
っていうか最初っからArray使えばよかったな。
531 名前：１３２人目の素数さん [04/11/16 01:28:53]: NIntegrate[Exp[-t^2], {t, 0, ∞}]
を数値積分すると結果は0.886227と出ます。

NIntegrate[Exp[-(t+a)^2], {t, 0, ∞}])　ここでaは0以外の整数。
この結果を数値(実数×aの形)で求めたいのですがmathematicaでは無理なのでしょうか？
もし可能でしたらその方法、不可能でしたら可能と思われるｿﾌﾄを教えていただきたいです。
よろしくお願いします。
532 名前：１３２人目の素数さん mailto:sage [04/11/16 01:40:49]: >>531
どんなソフトでも不可能。実数×aなんて形にはならないから。
aの関数にしたいならIntegrateがいいだろう。
NIntegrateしてフィッティングもできなくはないが。
533 名前：531 [04/11/16 01:44:53]: 解答ありがとうございます。
Integrateで計算するとインテグラルが残ってしまうんです・・。
フィッティングとは組み込み関数のfitのことでいいんでしょうか？
ヘルプを熟読してみます。
534 名前：532 mailto:sage [04/11/16 02:02:52]: >>533
Integrate[Exp[-t^2], {t, a, \[Infinity]}] といっしょでしょ。
535 名前：531 mailto:sage [04/11/16 02:12:40]: >>532
なんどもありがとうございます。
534の式だと結果に積分記号が残ったままで求めたい結果(数値だけ)が得られないんですよね・・。
他の方法を模索してみます・・。
536 名前：532 mailto:sage [04/11/16 03:12:12]: >>535
バージョンいくつ？4.2だったら
(-1/2)Sqrt[Pi](-1 + Erf[a])
と計算されるんだが。
537 名前：531 [04/11/16 20:16:49]: >>536
返答遅れてすみません・・。
バージョンは５です。Erfは誤差関数でインテグラルが残ってるんですよね・・。
538 名前：１３２人目の素数さん mailto:sage [04/11/17 00:58:51]: 意味が分からん。Erfの内部実装が積分だと知ってるといってるんだろうか。
x^2 = Integral[t, {t, 0, x}]と書けるんだが
x^2にはインテグラルが残っているんだろうか。
539 名前：１３２人目の素数さん mailto:sage [04/11/17 01:00:01]: Integrate[t, {t,0,x}] / 2
540 名前：１３２人目の素数さん mailto:sage [04/11/17 01:00:50]: Integrate[t, {t,0,x}] * 2
541 名前：１３２人目の素数さん mailto:sage [04/11/17 01:22:18]: Integrate[t, {t,0,x}] * Integrate[t, {t,0,2}]
542 名前：１３２人目の素数さん mailto:sage [04/11/17 01:30:45]: よーし、積分してから微分しちゃうぞー。
543 名前：１３２人目の素数さん mailto:sage [04/11/17 16:20:29]: 例えば Mathematica上で
In[2]:= Solve[ a^2 - 1 == 0 , a ]
と入力すると
Out[2]= {{a -> -1}, {a -> 1}}
が返ってきますよね？ですがこの時に a の 1 つ目の解を参照するにはどのような
記述をすれば良いのでしょうか？

上の例で言えば、

一つ目の解[ Solve[ a^2 - 1 == 0 , a ] ] * 10

のようにして、
Out[2]= -10
が返るようにしたいのです。
544 名前：１３２人目の素数さん mailto:sage [04/11/17 18:50:58]: >>543
(a /. Solve[a^2 - 1 == 0, a])[[1]] * 10
545 名前：１３２人目の素数さん [04/11/22 12:56:21]: ウォルフラムマークって何のためにあるの？
546 名前：１３２人目の素数さん mailto:sage [04/11/22 23:53:16]: Mathematica で f[x,y]=K （K=1～10までの等高線）の書き方が分かりません。

ContourPlot[f[x, y], {x, -5, 15}, {y, -5, 15}, Contours -> 10, PlotRange -> {0, 10}]

これを実行すると、等高線が10本出ますが、f[x,y]=1～10 なのかどうか、、、

どのようにしたらよいか教えてください。
547 名前：１３２人目の素数さん mailto:sage [04/11/23 02:14:12]: >>546
ImplicitPlot
548 名前：１３２人目の素数さん [04/11/23 13:44:55]: Mathematica の中・上級者向けで，
Mathematica のプログラミングを本格的に学べるいい本ある？
549 名前：546 mailto:sage [04/11/23 15:22:47]: >>547
ありがとう。ImplicitPlot でできました。
550 名前：１３２人目の素数さん [04/11/23 21:41:44]: arctanのグラフで，pi/2から-pi/2に突然飛ぶんですが，連続になるように描けませんか?
551 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/23 21:43:49]: Re:>550 そんなはずはない。入力した関数(式といったほうが分かりやすい？)を書いてくれ。
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｝という形で取り出したいのです。
もしよろしければ教えてください。
653 名前：１３２人目の素数さん mailto:sage [04/12/26 00:29:25]: 5.1 は、結構いろんな機能が入ってるな。v4->v5の時よりも
変化が大きいんじゃない？

この機能が楽しそうだなと思ったり。
www.wolfram.com/products/mathematica/newin51/equationtrekker.html
654 名前：１３２人目の素数さん mailto:sage [04/12/26 00:36:27]: >>652
最大数が複数個ある場合を考慮に入れるなら

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]

仕組みは簡単で、各セルの左側の値で大きい順にソートして最も大きいものを選ぶだけ。
純関数を使っているから見にくいが、やっていること自体は単純。
655 名前：１３２人目の素数さん mailto:sage [04/12/26 02:10:06]: >>651
5.1はWinでは x86版、　Linuxではｘ86版、x86-64bit版(athlon64)、Itanium版が探せばある。
656 名前：１３２人目の素数さん mailto:sage [04/12/26 20:45:15]: >>654
ありがとうございます。まだ使い始めて日が浅いので。
こうやればいいんですね。よく考えてみると単純ですね。
657 名前：１３２人目の素数さん [04/12/27 16:27:39]: >>646
どなたかご教授願います。
658 名前：１３２人目の素数さん [04/12/30 10:31:23]: age
659 名前：１３２人目の素数さん [04/12/30 14:15:04]: >>657
test1 = 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]
660 名前：１３２人目の素数さん [04/12/30 14:59:44]: 395
661 名前：１３２人目の素数さん mailto:sage [04/12/30 21:04:39]: >>657
xご教授 o御教示

mathよりjapaneseやれよｗ
662 名前：１３２人目の素数さん mailto:sage [05/01/01 00:00:25]: 本年初カキコ。
663 名前：１３２人目の素数さん [05/01/01 15:14:48]: >>659
ありがとうございます。
Show[test1,test2]で表された図を印刷したいんですが、図だけをコピーすることはできないでしょうか。
枠線が邪魔なんです。

>>661
＿|￣|○
664 名前：１３２人目の素数さん [05/01/02 20:26:57]: mathematica for linux
を中古で安く譲ってくれる方はいますか？

もし合法でないなら無視して下さい。
合法的に安く入手したいだけですので。
665 名前：１３２人目の素数さん mailto:sage [05/01/03 00:17:33]: 放送大学入学して、アカデミック（？）が一番安いのかな・・・。
仕事で使うのならまずいけど。
666 名前：１３２人目の素数さん mailto:sage [05/01/03 00:42:07]: >>664
へんなシト。英語版ならクレカを持っていれば、ダウソ版が1880ドルで買える。　これが
一番安いと思う。　Mathematicaは業務用ソフトではないし、値段はこんなものか。
日本語版は確か３５万くらいする。
667 名前：１３２人目の素数さん [05/01/03 14:43:19]: mathomaticaというものを見つけました。
gnuです。mathematicaとの互換性は？
668 名前：１３２人目の素数さん mailto:sage [05/01/03 21:50:48]: 別のソフト
669 名前：１３２人目の素数さん [05/01/03 22:04:31]: ざんねんん。ですね。
Mathematica恐怖症になってしまった。
やる気がおきない。
670 名前：１３２人目の素数さん [05/01/07 04:09:53]: age
671 名前：１３２人目の素数さん mailto:sage [05/01/07 10:55:52]: オークションサイトを探すとか
672 名前：１３２人目の素数さん [05/01/07 11:35:09]: 18xy-x^3-y^3=0
のグラフの輪っかになってる部分の面積を求めたいんですが、どうしたらよいでしょうか
673 名前：１３２人目の素数さん mailto:sage [05/01/07 12:23:17]: 高校の数学の教科書があるなら、是非見て下さい。数IIICのどっちかで理解出来ると思う。

グラフはこれで見れるね。
<< Graphics`ImplicitPlot`
ImplicitPlot[18x y - x^3 - y^3 == 0, {x, -10, 10}]

で、面積は、
18 x y - x^3 - y^3 < 0
の範囲であることは分かりますか？
674 名前：673 mailto:sage [05/01/07 12:37:19]: 陰関数の積分の方法を忘れてしまった＿|￣|○

Sin[x]やCos[x]を使うんだっけ…
675 名前：１３２人目の素数さん mailto:sage [05/01/07 20:32:06]: >>673
>0 の間違いでは？

>>674
version5以降なら、↓で計算できる。でも、式の簡略化に失敗してるなあ。

<< Calculus`Integration`
Integrate[Boole[18 x y - x^3 - y^3 > 0], {x, 0, 10}, {y, 0, 10}]
676 名前：673 mailto:sage [05/01/07 20:42:01]: >>675
訂正サンクス、やり直さなきゃいけないのは自分の方だった。

N[Integrate[Boole[18 x y - x^3 - y^3 > 0], {x, 0, 10}, {y, 0, 10}] ]
で値がでましたね。

54.
677 名前：１３２人目の素数さん [05/01/09 10:57:34]: ParametricPlot

Is there such function in the Mathematica?
678 名前：１３２人目の素数さん [05/01/09 11:24:36]: なんとか２万円くらいにまからんか？
679 名前：１３２人目の素数さん mailto:sage [05/01/09 12:14:23]: >>677
Yes, there is.

such as "ParametricPlot"

ParametricPlot[ { Sin[t],Cos[t] } , { t, 0, 2Pi } ];

It makes the unit circle (sin^2 + cos ^2 = 1)

>>678
学校の研究室で買って貰え。
680 名前：１３２人目の素数さん mailto:sage [05/01/09 15:23:19]: >>679

thnx
681 名前：１３２人目の素数さん mailto:sage [05/01/09 18:06:55]: >>678
学生だったらスチューデント版がたしかそんなもんじゃね？
682 名前：１３２人目の素数さん mailto:sage [05/01/13 19:12:14]: マセマティカにやらせて久しぶりにまともな計算結果を得た

やっとやっと

なんでもやらせようとするから駄目なんだな。
683 名前：１３２人目の素数さん [05/01/16 09:01:14]: 「ぱすじぇね」で出るパスは英語モード用だけど、
日本語モード用のぱすじぇねは無いのかな・・・

今Ver5.1最新版を英語モードで使っているのだが、
そのパスでは当然日本語モードへ出来ない。
少しでも日本語環境へ近づけるため、
プログラム本体の到る所にあるJapaneseフォルダーを
Englishフォルダーへコピーさせる荒業を試みたけど、
入力パッド画面・ヘルプ（中身のみ）は日本語化出来た
けど他はダメでした。
肝心の本体をどうすれば日本語に出来るのか・・・
684 名前：１３２人目の素数さん [05/01/16 10:19:03]: マセマティカ４で計算や表示できる関数が
マセマティカ５【日本語版】になると
きかいせいどがどうのこうのとかIntegrandがどうのおうの
とかうだうだいってきて、
あげくに計算しない。あほか。今度暇なときに
ボッフラムに文句をいうてやるわ。30万円近くしたのに。・・・・
685 名前：１３２人目の素数さん [05/01/17 16:03:54 ]: mathlinkを使って、計算部分をmathematicaからCの実行ファイルに渡して計算させて、
結果をmathematicaに戻すという事をやりたいんですが、

Cソースファイルのコンパイルの部分でつまづいています・・・。

私の環境はWindowsXP、Mathematica4.1、MicrosoftVisualC++5.0です。

①ウルフラムリサーチ
support.wolfram.co.jp/mathematica/mathlink/updates/mldk.html
からWindows用のMathLinkデベロッパキットを持ってくる。

②Files\CompilerAdditions\mldev32以下のファイルを
C:\Program Files\DevStudio\VC
以下にコピペ.。

③コマンドプロンプトから
～～/mprep f.tm -o f.c
と実行する。

というところまでやって、出来たf.cをVisualC++でコンパイルは出来るんですが、
ビルドをしようとすると
f.obj : error LNK2001: 外部ｼﾝﾎﾞﾙ "_MLPutInteger@8" は未解決です
f.obj : error LNK2001: 外部ｼﾝﾎﾞﾙ "_MLPutFunction@12" は未解決です
f.obj : error LNK2001: 外部ｼﾝﾎﾞﾙ "_MLNewPacket@4" は未解決です
・
・
・

というエラーが出て実行ファイルが作れません。
CのソースファイルはMathematicaブックの２．１２．３にあるものをそのまま使ってます。
やり方がわかるひと教えてください。
686 名前：685 [05/01/17 23:08:33 ]: 1 マイクロソフト開発者スタジオを始動します。

2 Fileメニューから、Newを選びます。

3 New箱の中では、Projectsタブをクリックします。

4 位置には、MathLinkExamplesのあるフォルダをタイプします。

5 Project Nameテキストボックスに、名前ADDTWOをタイプします。

6 Projectから、Win32 Applicationを選択します。

7 Win32が選択されると確認します。

8 OKボタンをクリックします。

9 ProjectメニューのProjectへ追加サブメニューから、Filesを選びます。

10 「MathLinkExamples」というディレクトリの中の"addtwo"ディレクトリを選択
687 名前：685 [05/01/17 23:09:52 ]: 11 File Name textboxに、ファイル名addtwo.cをタイプします。

12 OKボタンをクリックします。

13 二度ステップ9～12を繰り返してください: 一回目のファイル名addtwo.tm、２回目のファイル名addtwotm.cとする。

14 クリック、はい、"addtwotm.c"の参照を加える確認するボタン。

15 ProjectメニューのProjectへ追加サブメニューから、Filesをもう一度選びます。

16 今回はファイル名C:\Program Files\DevStudio\VC\lib\ml32i2m.lib。

17 OKボタンをクリックします。

18 Projectメニューから、設定を選びます。

19 Custom Buildタブをクリックします。

20 左の設定リストの中の、ADDTWOフォルダーアイコンを広げます。
688 名前：685 [05/01/17 23:10:31 ]: 21 "addtwo.tm"というファイルを選択します。

22 説明に、「mprepingなど」の注意をタイプしてください…

23 Build Commandsリストの中では、最初の空白のエントリーをクリックしてください、そ
して、タイプする:　mprep addtwo.tm-o addtwotm.c

24 Output Filesリストの中では、最初の空白のエントリーをクリックしてください、そし
て、タイプする:　. \addtwotm.c

25 OKボタンをクリックします。

26 Buildメニューから、Build ADDTWO.EXEを選びます。
689 名前：685 [05/01/17 23:11:35 ]: 685の③は無視してください。

その後色々調べて
①と②の後に
↑
の手順に従ってビルドを試みてみました。

その結果
--------------------構成: ADDTWO - Win32 Debug--------------------
ｶｽﾀﾑﾋﾞﾙﾄﾞを実行中
ｺﾝﾊﾟｲﾙ中...
addtwo.c
addtwotm.c
ﾘﾝｸ中...
..\..\..\..\..\..\..\..\..\DevStudio\VC\lib\ml32i2m.lib : fatal error LNK1106: ﾌｧｲﾙが無効であるか､またはﾃﾞｨｽｸがいっぱいです: 0x3ee7a0e7 にｼｰｸできません
link.exe の実行ｴﾗｰ

ADDTWO.exe - ｴﾗｰ 1、警告 0

というエラーが出てしまいます。
原因が分かる方いませんでしょうか・・・。
690 名前：685 [05/01/18 12:27:55 ]: いろいろやってはみたけれど、結局原因が分からなかったので
Microsoft　Visual　C++を5.0から6.0にしてやってみました。
先に結果を書くとうまくいきました。

上に書いた手順は5.0用なので6.0では少し手順が変わります。
MathLinkデベロッパキットの
Files\Documentation\English
に5.0、6.0共に詳しいやり方が乗っているので見てみるといいかも。

後でMathlinkを使う人のために私の環境と出来た場合と出来なかった場合について書いておきます。
駄目な環境：WinXPSP2　MicrosoftVisualC++5.0　Mathematica4.1
駄目な手順（原因は不明）
①ウルフラムリサーチ
support.wolfram.co.jp/mathematica/mathlink/updates/mldk.html
からWindows用のMathLinkデベロッパキットを持ってくる。

②Files\CompilerAdditions\mldev32以下のファイルを
C:\Program Files\DevStudio\VC
以下にコピペ.。

③>>686-688の手順でビルト失敗。
上にも書いたけど>>686-688の手順はMathLinkデベロッパキットの
Files\Documentation\English
に載ってる。
691 名前：685 [05/01/18 12:28:57 ]: 良かった環境：WinXPSP2　MicrosoftVisualC++6.0　Mathematica4.1
良かった手順
①上と同じ

②Files\CompilerAdditions\mldev32以下のファイルを
C:\Program Files\Microsoft Visual Studio\VC98
以下にコピペ。

③デベロッパキットに載っている6.0の手順に従ってビルト成功。
692 名前：１３２人目の素数さん mailto:sage [05/01/18 13:29:46 ]: おめでとう + グッジョブ！
693 名前：１３２人目の素数さん [05/01/19 02:06:52 ]: Mathematicaでのsubroutine命令Module、Blockの相違点や長短を表と例で示せと
課題であるのですが、まったくわかりません。（ちなみに私は、学校の実験程度で
やっているだけなので）
694 名前：１３２人目の素数さん mailto:sage [05/01/19 04:05:29 ]: BlockとModuleではスコープにおける名前解決の仕方が異なる。

Moduleの中で宣言された変数は、Moduleのスコープ内のみで有効でその外に
名前解決がされることはない。

一方、Blockの場合は、スコープ内でのみ有効な値を取ることも出来、名前解決が
スコープ内で出来ない場合はスコープの外に（例えば、グローバル・スコープに）
名前解決を求めることが出来る。

スコープとはMathematicaで言えば、Module[], Block[] の[]の中や
グローバル・セッションを指す。
変数とはポインタのことで、ポインタは何かオブジェクトを指している。
オブジェクトとは、この場合、整数、浮動小数点数、式、etc。
Mathematicaは実行中、Global Scope（global session）を含む各スコープで
ポインタとオブジェクトの表を持っている。　式の評価中にポインタの
lookupが必要になったとき、Moduleではその評価式のあるスコープの表のみ
を探す。　一方、Blockではそのスコープの表で何らかのオブジェクトを
指していなければ、そのBlock式を評価しているスコープの表でlookupを
行う。

Blockの長所は便利であること。短所は、注意しないと思わぬ結果を招くことが
挙げられる。　インタラクティブなモードで式を評価するときには良い。
Moduleはライブラリなどを書くときに、外のスコープの値によって結果が
異なることを避けることが出来る。
695 名前：１３２人目の素数さん mailto:sage [05/01/20 16:22:55 ]: 5.0から5.1にしてみた。
違いがまだほとんどわからないけど、計算がわずかに早くなったような気がする。
696 名前：１３２人目の素数さん [05/01/21 14:10:40 ]: 『Inverse::luc』
とはどういう意味の警告ですかね？
697 名前：１３２人目の素数さん mailto:sage [05/01/21 16:15:48 ]: 逆行列を定義出来ない行列を、逆行列にしようとしたの？
698 名前：１３２人目の素数さん mailto:sage [05/01/21 19:23:53 ]: General::luc
行列expr に含まれる悪条件より，結果
には重大な数値的エラーが含まれている可能性があります．

・線形方程式の数値解の計算に用いるアルゴリズムの失敗によって表示される．
・もし思わぬ場面でこのメッセージが表示されたら，テクニカルサポートまでご連絡ください．
699 名前：１３２人目の素数さん [05/01/23 17:35:30 ]: ３ＤＰｌｏｔで双曲面を画こうとする時、±√を＋√と－√に分けて作図すると、＋と－のつなぎ目の部分が途切れてしまってうまく画けないのですが、どのようにすればいいのでしょうか？
４．０を使っていて、±ってはじめに出ているパレットには無いのですが、何処からか呼び出せるのでしょうか？
初歩過ぎる質問だとは思いますが、どうか教えていただきたいです。よろしくお願いします。
700 名前：１３２人目の素数さん mailto:sage [05/01/23 21:11:16 ]: こういうこと？
ParametricPlot3D[{x, y, x^2/1.3^2 - y^2/1.5^2}, {x, -1, 1}, {y, -1, 1}];

具体的な式書いてよ。
701 名前：１３２人目の素数さん [05/01/24 12:20:57 ]: 次の論理式について積和標準形を求めよ。トートロジーかどうかを判定せよ。

（１）or[imply[p,q],imply[q,p]]
（２）or[and[p,q],and[q,r],and[r,p]]～imply～and[or[p,q],or[q,r],or[r,p]]
（３）and[or[p,q],or[q,r]or[r,p]]～imply～or[and[p,q],and[q,r],and[r,p]]
（４）imply[or[p,q],and[q,r]]～imply～imply[p,r]

という問題です。mathematicaを使ってみたら、

（１）or[not[p],q,not[q],p]
（２）or[and[not[p],not[q],not[r]],and[not[p],not[r]],and[not[p],not[q]],and[not[q],not[r]],and[not[q],not[p]],and[not[q],not[r],not[p]],and[p,q,r],and[q,r],and[p,r],and[p,q],and[q,p],and[q,r,p]
（３）or[and[not[p],not[q]],and[not[q],not[r]],and[not[r],not[p]],and[[p,q],and[q,r],and[r,p]]
（４）or[and[p,not[q]],and[q,not[q]],and[p,not[r]],and[q,not[r],not[p],r]]

となりました。ここからどうやって積和標準形を求め、トートロジーの判定をすればいいのかわかりません。教えてください。お願いします。
702 名前：１３２人目の素数さん mailto:sage [05/01/24 21:01:27 ]: >>701
積和標準形はなんだかしらないけれど、トートロジーは簡単だ。
どういう真偽値の組合せでも常に真になることを言えばいい。

Or[Implies[p, q], Implies[q, p]] /.
Join @@ Outer[{p -> #1, q -> #2} &, {True, False}, {True, False}] //
Apply[And, #] &
703 名前：１３２人目の素数さん [05/01/24 22:04:55 ]: >>700
レスありがとうございます！
所謂、二葉双曲面で色んなパターンが描きたいのですが…
例えば
Z^2=3X^2-Y^2-3/4とかZ^2=-X^2+3X-2sqrt(3)XY+2Y^2-sqrt(3)y+1/4
この様な式のグラフィックとはどう描けば良いのでしょうか？
ネットで調べてみると難しいパラメータを利用してるみたいでさっぱりわかりませんでした…
704 名前：１３２人目の素数さん [05/01/25 17:43:15 ]: x+2x+3x
のような多項式を

3
∑nx
n=1

と変換させるにはどうしたらよいのでしょうか？
Outに∑の形で出したいのですが可能でしょうか。
705 名前：１３２人目の素数さん mailto:sage [05/01/25 17:53:01 ]: >>704
> Outに??の形で出したいのですが可能でしょうか。
それ自体は可能だ(おれの環境では見えないが、\sumだと仮定して)。
HoldFormを使う。

> x+2x+3x
> のような多項式を

こっちは意味がわからない。
1) 6xと足して欲しくないということか？
2) x + 2 x^2 + 3 x^3 の間違いか？
3) 係数がどのようなnの関数であるか既知なのか？
706 名前：１３２人目の素数さん [05/01/25 17:55:40 ]: 704です。

6xと表示させたいのでは無く、
∑とn (n=1,2,3・・・・)の形で表示させたいんですが。
nがどのような形でかかってくる関数かは既知です。
707 名前：１３２人目の素数さん mailto:sage [05/01/25 18:45:50 ]: >>706
じゃ HoldForm[Sum[n x, {n, 3}]]
708 名前：１３２人目の素数さん [05/01/25 19:09:22 ]: すいません
nx にあたる部分が未知です。

ｎで∑にまとめられそうな多項式を入力として

∑なんたら
n=1

と出力するにはどうしたらいいでしょうか？
∑の形にMathematicaに自動で処理させたいのです。
709 名前：１３２人目の素数さん [05/01/25 19:33:05 ]: 3
In=∑x^n
n=1

Out=x+x^2+x~3

の逆で

Int=x+x^2+x~3

　 3
Out=∑x^n
　n=1

がしたいってことなんですが・・・
710 名前：１３２人目の素数さん mailto:sage [05/01/25 21:41:42 ]: >>709
それは無理だと思うんだが…

あと、あなたの書いてる式は多項式ではないよ。xという1つの項しか使ってないから。
711 名前：１３２人目の素数さん [05/01/26 11:30:02 ]: 調べてもわからなかったのでご意見ください．

Sin[A Cosθ]

を入力としたとき，
A<<１
ならばという条件で近似させ

A Cosθ

を出力として取り出したいのです．
このように，ある条件を与えて近似する方法，ご存知の方ご享受ください．
712 名前：１３２人目の素数さん mailto:sage [05/01/26 14:32:08 ]: >>709
要領を得ない書き込みだな。ミスも多いし。

例えば
x + 9 x + 5 x^ + 2 x^3
ならどうnの式で表せると思う？

係数が未知なら逆に戻すのは不可能。
If[n == 1, 9, If[ ...]]みたいな係数でよければできるがな。
713 名前：１３２人目の素数さん mailto:sage [05/01/26 14:39:00 ]: >>711
SeriesCoefficient[Series[Sin[A Cos[t]], {A, 0, 1}], 1]

>>712
自分でミスってりゃ世話ないな。
1 + 9 x + 5 x^2 + x^3
714 名前：１３２人目の素数さん [05/01/26 18:11:21 ]: ベッセル関数の行列式の絶対値の表示の仕方を教えてください。
ちなみに今は12×12の行列式をmatと置き、それをDet[mat]としています。
715 名前：１３２人目の素数さん [05/01/28 00:18:08 ]: 二葉双曲面で色んなパターン
例えば
Z^2=3X^2-Y^2-3/4とかZ^2=-X^2+3X-2sqrt(3)XY+2Y^2-sqrt(3)y+1/4
この様な式のグラフィックとはどう描けば良いのでしょうか？
ParametricPlot3D[{x, y, Sqrt[3X^2 - Y^2 - 3/4], -Sqrt[3X^2 - Y^2 - 3/4]}, {x,
-3, 3}, {y, -3, 3}];
では描けないのですがやはりパラメータを使わなければ描けない物なんでしょうが？
716 名前：１３２人目の素数さん mailto:sage [05/01/28 06:59:12 ]: << Graphics`
InequalityPlot3D[3x^2 - y^2 - z^2 - 3/4 >= 0, {x, -5, 5}, {y, -5, 5}, {z, -5, 5}]
717 名前：716 mailto:sage [05/01/28 07:00:35 ]: やりなおし
<<Graphics`
InequalityPlot3D[3x^2-y^2-z^2-3/4>=0,{x,-5,5},{y,-5,5},{z,-5,5}]
718 名前：１３２人目の素数さん mailto:sage [05/01/28 13:51:25 ]: >>717
なんか変な図だね。円錐に見える。

こっちのほうがきれいだと思う。
ContourPlot3D[z^2 - ( 3 x^2 - y^2 - 3/4), {x, -5, 5}, {y, -5, 5}, {z, -5, 5},
PlotPoints -> {5, 5}]
719 名前：１３２人目の素数さん [05/01/29 00:41:39 ]: >>717
>>718
本当にありがとうございます！大変参考になりました。
なぜか717の式をコピペしてもグラフィックが出力されないのですが…
4.0では無理って事はないですよね？
720 名前：１３２人目の素数さん mailto:sage [05/01/29 02:07:52 ]: >>719
バージョン4.2か5以降だったような気がする。
721 名前：１３２人目の素数さん [05/01/29 10:50:06 ]: Ver4.1使ってます。
Plot[hoge[x],{x,-100,100}]
とした場合実際には最大で何ポイントで計算されてるんでしょうか？

PlotPoints->25 、 MaxBend->15 、 PlotDivision->30　が標準らしいのですが、

25*30 = 750 ポイントなんでしょうか？
722 名前：１３２人目の素数さん [05/01/30 00:03:01 ]: 誰か自然現象・社会現象のなかからテーマを一つ選び，問題を定義し，ルールを書き，シミュレーションするプログラム(notebookの形)を作成してください．手法として，ルールと記述し，状態を変化させそうなんですが。。。。
723 名前：１３２人目の素数さん mailto:sage [05/01/30 00:19:31 ]: ここまで投げやりなのも珍しいな？　低凶？
724 名前：７２２ [05/01/30 02:09:10 ]: いや一橋のものです。。
725 名前：１３２人目の素数さん mailto:sage [05/01/30 02:36:58 ]: 3次元曲線の線を太くするってどうやればできるんですか？
726 名前：１３２人目の素数さん [05/01/30 03:29:00 ]: Plot[{Sqrt[48x^2 - 0.979], -Sqrt[48 x^2 - 0.979]}, {x, 0, 8},
AspectRatio -> Automatic, PlotRange -> {{-7, 7}, {-7, 7}}]
とプロットすると

　　　Plot::plnr : x = 3.3333333333333333333333`*^-7
において\!$\@\(\(\(48\\ x\^2$\) \
- 0.979`\)\)は機械精度の実数ではありません．"\)

Plot::"plnr" : x = (0.07539936143743177`\)において\!$\@\(\(\(48\\ x\^2$\) - \
0.979`\)\)は機械精度の実数ではありません．"\)

Plot::"plnr" x = \!$0.11815455710265632`$において\!$\@\(\(\(48\\ x\^2$\) - \
0.979`\)\)は機械精度の実数ではありません．"\)

General::"stop": "計算中，Plot :: plnrのこれ以上の出力は表示されません．"

こんな風に表示されて、双曲線の頂点付近がプロットされないのですが、どのように解決すればよいのでしょうか？
727 名前：１３２人目の素数さん mailto:sage [05/01/30 03:34:55 ]: >>725
Plot3D[BesselJ[2,x y],{x,-10,16},{y,-1,1},MeshStyle->{Thickness[.006]}]
728 名前：１３２人目の素数さん mailto:sage [05/01/30 08:43:16 ]: >>726
＞機械精度の「実数」ではありません．"
729 名前：１３２人目の素数さん mailto:sage [05/01/30 09:30:57 ]: >>722
現代の特殊出生率と、過去のものを比較してみれば？

>>724
嘘をつくな嘘をｗ
730 名前：７２２ [05/01/30 12:35:45 ]: >>729
やり方が全くわからないので、ここに載せてもらえないでしょうか？？お願いします。。。
731 名前：１３２人目の素数さん mailto:sage [05/01/30 13:43:38 ]: >>727
MeshStyleか！　どうも㌧ｸｽです
PlotStyleが使えないのどうして？ってずっと思ってたです
732 名前：１３２人目の素数さん mailto:sage [05/01/30 16:42:00 ]: >>731
Options[Plot3D]とかやれば、使えるオプションの一覧が出るよ
733 名前：１３２人目の素数さん [05/01/30 21:05:41 ]: ??Plot3D
でいいんじゃね？

で、質問なのですが、
f[x_,y_,z_,t_]=で関数定義して、
これを
g[t_] = f[x_,y_,z_,t_]をx,y,zそれぞれ±無限大で積分
としてg[t_]をプロットしたいのですが、時間がかかりすぎて困ってます。
最適化するにはどうしたらいいんでしょうか？
また、g[t_]を計算した時点で保存して次回からいきなりPlot[g[t]]させる方法
ってあるんでしょうか？
734 名前：１３２人目の素数さん [05/01/30 22:35:23 ]: listplot
データを間引きすれば？
735 名前：733 mailto:sage [05/01/31 00:03:20 ]: >>734さん

レスサンクスです。

PlotPoints,PlotDivision,Maxbendをいじってもいるんですが、
速くすると精度は落ちるしってパターンでしかないので
ちょっと迷ってます。とりあえずPlotPointsをぐっと減らしてmaxbendもあげて概形を
見てから計算するようにして一晩放置して計算させたのに無駄だったってことだけは下げてます。

やっぱいいPC買うしかないんですかね今は、AthlonXP@1.7GHz orz..
736 名前：１３２人目の素数さん mailto:sage [05/01/31 08:35:26 ]: 数学で生きていくのであれば、商売道具に金かけた方が幸せでしょうね。
737 名前：１３２人目の素数さん mailto:sage [05/01/31 10:30:46 ]: >>735
いいPC買ったって2倍速にもならんだろ。

Integrateをそのまま使ってて、
記号積分に挑戦して時間がかかっているだけじゃないのか？
738 名前：735 mailto:sage [05/01/31 15:02:09 ]: NIntegrateにしてます。

理想的には不貞積分した時に式がどうなるか出してほしいのですが
それはあきらめました。
739 名前：１３２人目の素数さん mailto:sage [05/01/31 18:04:33 ]: >>738
不貞積分(w

では、Plot[Evaluate[NInte...],{t,...}]にして(それほど効果がないかもしれん)、
NIntegrateのオプションをいじる。

オプションをいじってもNIntegrate自体に時間がかかるんなら、
最適化も何も無理。NIntegrateを自分で書くぐらいしか手はない。
740 名前：１３２人目の素数さん mailto:sage [05/02/01 20:14:33 ]: 要素数５の0,1のみからなるリストを全種類作りたいのですが、
どのようにすればよいですか？
{0,0,0,0,0}、{0,0,0,0,1}、{0,0,0,1,0}、{0,0,0,1,1}・・・
（全部で32個）
741 名前：１３２人目の素数さん mailto:sage [05/02/01 21:06:43 ]: 要素が5個だけなら
たぶん手間的に32個書いたほうが楽だと思うけど…
742 名前：１３２人目の素数さん mailto:sage [05/02/01 22:23:19 ]: For文で要素をi,j,k,l,mとして、
for (int i=0;i<=1;i++){for (int j=0;j<=1;j++){　略
{i, j, k, l, m}
でどう？
743 名前：１３２人目の素数さん mailto:sage [05/02/01 22:50:59 ]: マセマティカで
FORはつかっちゃだめでしょう？
せめてDOにしなさい
744 名前：１３２人目の素数さん mailto:sage [05/02/01 22:56:08 ]: f[list_, elems_] := Map[Prepend[list, #] &, elems]
g[list_, elem_] := Flatten[Map[f[#, elem] &, list], 1]

Nest[g[#, {0, 1}] &, {{}}, 5]
でどう？　
745 名前：744 mailto:sage [05/02/02 00:37:41 ]: Moduleを使ってまとめれば、

genSeq[elms_, n_] :=
Module[
{f, g},
f[list_, elems_] := Map[Prepend[list, #] &, elems];
g[list_, elem_] := Flatten[Map[f[#, elem] &, list], 1];

genSeq[{0, 1}, 5]
746 名前：744 mailto:sage [05/02/02 00:40:52 ]: >>745
コピーを間違えた

genSeq[elms_, n_] :=
　 Module[
　　{f, g},
　　f[list_, elems_] := Map[Prepend[list, #] &, elems];
　　g[list_, elem_] := Flatten[Map[f[#, elem] &, list], 1];
　　Nest[g[#, elms] &, {{}}, n]
　 ]

で、

genSeq[{0, 1}, 5]
747 名前：１３２人目の素数さん [05/02/02 02:40:26 ]: Plot[{Sqrt[3 x^2 - 3], -Sqrt[3 x^2 - 3]}, {x, 1, 7}, AspectRatio -> Automatic,
PlotRange -> {{-7, 7}, {-7, 7}}]

では、問題なくグラフが画けるのに、
Plot[{Sqrt[3 x^2 - 3], -Sqrt[3 x^2 - 3]}, {x, 0, 7}, AspectRatio -> Automatic,
PlotRange -> {{-7, 7}, {-7, 7}}]
などと範囲を少し変えただけなのに「機械精度の実数ではありません」等と表示されてグラフが途切れて画かれてしまうのはなぜなんでしょうか？
５．０等新しいバージョンでは問題なく画けるのでしょうか？
748 名前：１３２人目の素数さん mailto:sage [05/02/02 03:39:44 ]: >>747
{x, 0, 7}とすると、0 <= x < 1で　3 x^2 - 3　< 0 だから、 Sqrt[3 x^2 - 3]が実数でない。
したがって、Mathematicaはエラーを出す。
749 名前：740 mailto:sage [05/02/02 10:23:57 ]: >>744
すごい！ありがとうございます！
早速使わせて頂きます。
ただ、Moduleを使うと、
{0, 0, 　, 0, 0}, {　, 0, 　, 0, 0}, {0, 　, 　, 0, 0}, ･･･
のように１が表示されなくなります・・・。
750 名前：１３２人目の素数さん mailto:sage [05/02/02 10:58:33 ]: >>749
>>746
751 名前：１３２人目の素数さん mailto:sage [05/02/02 12:14:45 ]: >>740
Outer使え。

Gen[n_] := Flatten[Outer[List, Sequence @@ Table[{0, 1}, {n}]], n - 1]

Gen[5]
752 名前：740 mailto:sage [05/02/02 13:21:19 ]: >>751
Outerをうまく使えないものか考えていました。
ありがとうございます。
すごくスッキリしました！

なかなかMathematicaプログラミング技術が上達しなくて、
日々奮闘しているので、こちらは大変勉強になります。
753 名前：１３２人目の素数さん mailto:sage [05/02/02 18:14:15 ]: >>752
そういうふうに、正しく(すっきりと)やろうとするのは大切だよ。
アドバイス:きれいなプログラムが書けるようになりたいなら、ForとかDoとかは使わないこと。
754 名前：１３２人目の素数さん [05/02/02 18:46:48 ]: >>748
初歩的な質問に親切に答えていただいてありがとうございます！
755 名前：１３２人目の素数さん mailto:sage [05/02/02 19:25:36 ]: 組み合わせを求める質問、あちこちで見たな。需要が多いみたい
だから、次のバージョンあたりで簡単にできる関数を入れてくれると
うれしい。
756 名前：１３２人目の素数さん [05/02/03 01:13:37 ]: ImplicitPlotって４．０では使えないのでしょうか？
<< Graphics`ImplicitPlot`
a := 1; b = 1.6; c := 1
p := 1; q := -1; r := -1
ImplicitPlot[ax^2 + 2bxy + cy^2 + px + qy + r == 0, {x, -3, 3}]
この様に入力しても「ImplicitPlotはプロテクトされています」と出てきてグラフィックが表示されないのですが…
757 名前：１３２人目の素数さん mailto:sage [05/02/03 15:37:03 ]: >>756
そのように入力したら変数名はaxやbxyなので、
エラーメッセージこそ違えど、動かないでしょうね。
で、本当はどう入力したんですか？
コピーして同じになるものを貼ってください。
758 名前：１３２人目の素数さん [05/02/03 19:00:17 ]: Mathematicaの質問とは少し的がはずれてますが、
VC#.Netを使用してMathLinkとの通信を行うソフトを開発しております。
しかし、結果の取得の部分でどうしてもうまくいかなく行き詰まっております。
その内容としては、オリジナルのパッケージを読み込み、評価式を送信すると、
グラフや数字の結果が得られるというものです。
その結果の中でも、Answer{int1,int2,...,intN}EndAnswer　で囲まれた部分を
受け取りたいのですがうまくいきません。
どうすればよいでしょうか？
よろしくお願いします。
759 名前：１３２人目の素数さん [05/02/03 21:39:49 ]: このプログラムをMathematicaで作動させたいのですがどう変形すればいいでしょうか？

またInvFFTは高速フーリエ変換です。
BEGIN
InitGauss(seed)
beta := 2*H + 1
FOR i := 0 TO N/2 - 1 DO
rad := power(i + 1, -beta/2)*Gauss()
phase := 2*3.14592*rand()/Arand
A[i] := rad*cos(phase)
B[i] := rad*sin(phase)
END FOR
InvFFT(A, B, X, N/2)
END

ブラウン運動をフーリエ変換で解析するのが目的です。
どうかお願いします。
760 名前：１３２人目の素数さん [05/02/04 00:15:51 ]: >>757
ありがとうございます！
例えば
<< Graphics`ImplicitPlot`
ImplicitPlot[x^3 + y^3-3xy== 0, {x, -2, 2},{y,-2,2},AxesOrigin->{0,0}]
この様な式が、本やネットに載っていたのでそのままかいても
does not have a single variable other than
とでてグラフ化できません。
この x^3 + y^3-3xy== 0をx^3 + y^3-3y== 0やx^3 + y^3-3x== 0にしたら出来るのですが…
なぜかxyの項が入るとプロットできません。
何かPackageがうまく組み込めていないとかあるのでしょうか？？
761 名前：１３２人目の素数さん mailto:sage [05/02/04 01:10:51 ]: >>760
だからxyはx yじゃないの。
axyはa x yじゃないの。
xyはx y == x * y とは別の変数なの。OK？
762 名前：１３２人目の素数さん mailto:sage [05/02/04 01:14:33 ]: xyじゃなくてx y(間にスペース)な。
763 名前：１３２人目の素数さん [05/02/04 01:55:26 ]: >>761
>>762
すみません。超初心者＆孤立無援なものでスペースがいる事を知れませんでした
ありがとうございました。
764 名前：１３２人目の素数さん mailto:sage [05/02/04 12:13:05 ]: >>763
マニュアルやオンラインヘルプくらい見ましょうよ…
765 名前：１３２人目の素数さん mailto:sage [05/02/04 23:55:52 ]: バグを発見してレポートを送ったらメールが来て、「Mathematica5.1では
正常に動作するようになりましたので、アップデートしてください」とな。
有料。ぉぃ。
766 名前：１３２人目の素数さん mailto:sage [05/02/05 00:51:28 ]: (´-`).。ｏＯ (そこでもう一押しして、無料でver upしてもらう！！)
767 名前：超初心者 [05/02/05 01:00:51 ]: 大学で触ったこともないのにいきなりレポートだされて困ってます。
いろいろ教えてください。お願いします。
まず一つ目です。
陰関数のグラフを書くにはどうしたらいいのでしょうか。
y=f(x）の形になっていれば書けるのですが、f(x
y)=0の形になってしまうとわかりません。
x^5+y^5-3xy=0を例にお願いします。
768 名前：１３２人目の素数さん mailto:sage [05/02/05 01:29:15 ]: >>767
<< Graphics`ImplicitPlot`

ImplicitPlot[x^5 + y^5 - 3x y == 0, {x, -2, 2}]
769 名前：超初心者 [05/02/05 01:45:00 ]: >>768ありがとうございます。
今日は遅いので続きは明日にききます。って明日メ切なんですけどね。
もしお暇でしたらよろしくお願いします。
770 名前：超初心者 [05/02/05 11:03:31 ]: 極座標r=acosθのグラフをかくにはどうしたらいいのでしょうか。
わかりません。お願いします。
771 名前：超初心者 [05/02/05 11:34:15 ]: 追加で767の問題で定数aがある場合はどうしたらいいのでしょうか？
例y2^-x^2(x+a)=0
772 名前：長所新車 [05/02/05 14:56:32 ]: 解決しました
773 名前：765 mailto:sage [05/02/05 19:45:13 ]: >>766
DSolve[...](...はある式)を評価するとMathKernelが落ちるっていう
バグなんだが。なんかバグフィックスに使われたような気がして
あまりいい気がしないコメントがResearch & Developmentって
部署のDirector名義で返ってきた。一押ししても無理やろうな。
基本的にWolframはケチっぽいし。
774 名前：１３２人目の素数さん mailto:sage 釣られてみるか [05/02/06 00:56:44 ]: >>771
Mathematicaはヘルプが充実しているからまずはそれを使え。
解説書もたくさん出ているから大きめの本屋で漁ってみて
手頃だと思うのを買ってみるといい。
自分の力で考えたり調べたりする前に安易に2chに頼るのは
止めとけ。
脳味噌腐るぞ。
775 名前：１３２人目の素数さん [05/02/07 18:53:11 ]: [Random[II]=(II+1)^(-3/2)*NormalDistribution,{II,1,100}]
ＩＩを１から１００まで繰り返したいのですが、うまくいきません。
どうすればいいでしょうか？
また正規分布の表し方はこれでいいのでしょうか？
776 名前：１３２人目の素数さん mailto:sage [05/02/07 19:13:54 ]: > ＩＩを１から１００まで繰り返したいのですが、うまくいきません。
変数を繰り返すとは意味がわからん。書いてある式も意味がわからん。

> また正規分布の表し方はこれでいいのでしょうか？
HelpのNormalDistributionの項を見よ。
777 名前：775 [05/02/07 20:06:20 ]: >すいいません。
ＩＩではなくiとして考えてください。
Ｃ言語でいうfor(i=1;i<=100;i++)がやりたいです。
このやり方でDoを使わずにできるときいたのですができませんでした。
778 名前：１３２人目の素数さん mailto:sage [05/02/07 21:10:02 ]: 文法も何も考えずにでたらめ書いて動くわけないだろ。
何を求めたいのかすらわからん。
779 名前：１３２人目の素数さん mailto:sage [05/02/07 22:20:10 ]: こういうことがしたいのか?かなりエスパー入ってるが。

<< Statistics`NormalDistribution`
R = Table[(II + 1)^(-3/2)*Random[NormalDistribution[0, 1]], {II, 1, 100}]
780 名前：１３２人目の素数さん mailto:sage [05/02/07 22:20:27 ]: >>777
あなたがやろうとしていることはC言語も知らないのに、
「printf関数とは何ですか」と質問しているに等しい
781 名前：777 [05/02/07 23:53:21 ]: >132人目の素数さん
やりたいことはまさにその通りです。
ありがとうございます。
もっと勉強してから出直します。
782 名前：779 mailto:sage [05/02/08 00:02:27 ]: どの名無しさんだか。
783 名前：１３２人目の素数さん [05/02/09 11:02:21 ]: 何のプログラミング言語でもいいから
２つの項を入力されたらmost general unifierを出力する
*(1,y) *(x,+(a,b)) -> {x->1,y->+(a,b)}
のようなものを作れ、という課題が出たのですが、mathematicaでもできるらしいと聞きました。
mathematicaではどうすればいいのかどなたか教えていただけないでしょうか。
784 名前：１３２人目の素数さん mailto:sage [05/02/09 12:43:37 ]: まともに動いてるプログラムに日本語のコメント入れたら構文エラー
で動かんくなった．(Ver. 4.0)SJISの処理でスキャナーがコケてるんだろう．
785 名前：１３２人目の素数さん mailto:sage [05/02/09 15:43:32 ]: >>783
"most general unifier"なんて知らないが、講義で求めるアルゴリズムを習ったんだろ？
そのまま実装してみろよ。
786 名前：１３２人目の素数さん mailto:sage [05/02/09 23:10:58 ]: >>783はﾏﾙﾁ
pc5.2ch.net/test/read.cgi/tech/1076418993/372
787 名前：１３２人目の素数さん mailto:sage [05/02/11 03:54:59 ]: 777にはリアルでウケたwww
788 名前：１３２人目の素数さん [05/02/11 09:55:11 ]: dataf=Table[f[m],{m,0,600}]

Do[WriteString[stmp, dataf[[m]]], "\n"], {m, 1, 601}]

でファイルに書き出しているのですが、出力が
{105.98}
{106.704}
{107.147}
{108.42}
{108.795}
みたいにカギ括弧がついてエクセルに持って行ったときに不便で仕方有りません。
かっこ無しで出力する方法って無いでしょうか？
789 名前：１３２人目の素数さん [05/02/11 15:05:59 ]: >>788
書き直しなさい、そのままでは実行できない。
790 名前：１３２人目の素数さん mailto:sage [05/02/11 17:51:47 ]: それっぽい関数がある気もするけど、こんな感じ?

(*PrintAsCSV[リスト(, 桁数)] : 何か書き出すぞ.*)
ClearAll[PrintAsCSV, hPrintAsCSV]
PrintAsCSV[x_, precision_:$MachinePrecision] := hPrintAsCSV[x, precision]
PrintAsCSV::notnumber = "変なもの`1`が含まれていますが何か?";
PrintAsCSV::nulllist = "空リストみっけ.";
hPrintAsCSV[x_, p_] /; ((! NumericQ[N[x]]) && (! ListQ[x])) || (Head[x] === Complex) := Module[{}, Message[PrintAsCSV::notnumber, "'" <> ToString[x] <> "'"]; "("<>ToString[x]<>")"]
hPrintAsCSV[{}, p_] := Module[{}, Message[PrintAsCSV::nulllist]; "(空リスト)"]
hPrintAsCSV[x_, p_] /; (! ListQ[x]) && (MachineNumberQ[x]) && (Head[x] =!= Complex) := ToString[InputForm[x, NumberMarks -> False]]
hPrintAsCSV[x_, p_] /; (! ListQ[x]) && (NumericQ[x]) && (Head[x] =!= Complex) := ToString[N[x, p]]
hPrintAsCSV[{x__}, p_] := hPrintAsCSV[x, p]
(*タブがデフォ.*)
hPrintAsCSV[x__, y__, z___, p_] := hPrintAsCSV[x, p] <> "\t" <> hPrintAsCSV[y, z, p]
(*んで, 改行条件の設定.コメントアウトで動作が変わるよん*)
hPrintAsCSV[{x__}, {y__}, z___, p_] := hPrintAsCSV[{x}, p] <> "\n" <> hPrintAsCSV[y, z, p]
hPrintAsCSV[x__, {y__}, z___, p_] := hPrintAsCSV[x, p] <> "\n" <> hPrintAsCSV[{y}, z, p]
hPrintAsCSV[{x__}, y__, z___, p_] := hPrintAsCSV[{x}, p] <> "\n" <> hPrintAsCSV[y, z, p]

(*デモ*)
{3, Sin[1], {x, {3}, {5/3, 2 + 3I, {0, 8}, 4.2866}, {{1, , 2}, {}, {4.68462134654532463543513465413246843213543546345}}}}
tmpstr = PrintAsCSV[%]
791 名前：790 mailto:sage [05/02/11 18:29:19 ]: 科学表記されて死ぬな。
9行目
○ToString[AccountingForm[N[x, p]]]
○ToString[FortranForm[N[x, p]]]
×ToString[N[x, p]]
我ながら暇人。
792 名前：１３２人目の素数さん mailto:sage [05/02/11 23:37:15 ]: リストの関数ってのは定義可能？たとえば

Risuto2[y_,z_]: =
CoefficientList[Normal[Series[Lc[x, y, z], {x, 0, 8}]]],x]

それとも
FirstCoeff[y_,z_]: =
CoefficientList[Normal[Series[Lc[x, y, z], {x, 0, 8}]]],x][[1]]
??
793 名前：１３２人目の素数さん mailto:sage [05/02/12 00:12:37 ]: >>792
リストを引数に取ることも戻り値にすることも出来る。
794 名前：788 mailto:sage [05/02/12 13:55:56 ]: 吸いません・・・　　すごいおちでした。

datafの定義の中の式を{}で囲んでました…
だめだもうむりぽ。
795 名前：１３２人目の素数さん [05/02/12 18:56:38 ]: 複数のリストをひとつのグラフに表示したいです。

l1 = Table[x, {x, 0, Pi, 0.1}];
l2 = Table[Sin[x], {x, 0, Pi, 0.1}];
ListPlot[{l1, l2}]

↑エラー。
どうすれば表示できますか？あほですいませんお願いします。
796 名前：１３２人目の素数さん mailto:sage [05/02/12 21:12:41 ]: >>795
<< Graphics`MultipleListPlot`
l1 = Table[x, {x, 0, Pi, 0.1}];
l2 = Table[Sin[x], {x, 0, Pi, 0.1}];
MultipleListPlot[{l1, l2}]
797 名前：796 mailto:sage [05/02/12 22:08:15 ]: >>795
もしパッケージを使いたくなかったら以下。DisplayFunctionのヘルプの
中の例にまんま載ってるし。Mathematicaはヘルプが最高の資料だぜぃ。
l1 = Table[x, {x, 0, Pi, 0.1}];
l2 = Table[Sin[x], {x, 0, Pi, 0.1}];
p1 = ListPlot[l1, DisplayFunction -> Identity];
p2 = ListPlot[l2, DisplayFunction -> Identity];
Show[{p1, p2}, DisplayFunction -> $DisplayFunction];
798 名前：１３２人目の素数さん mailto:sage [05/02/13 00:28:19 ]: >>Mathematicaはヘルプが最高の資料
同意
799 名前：１３２人目の素数さん [05/02/14 07:27:23 ]: << Graphics`MultipleListPlot`
lis={1,2,3};
MultipleListPlot[lis,lis/4]

を実行すると

MultipleListPlot[{1,2,3},{1/4,1/2,3/4}]

という結果しか帰って来ません。
どうすればグラフが表示されますか。
800 名前：１３２人目の素数さん mailto:sage [05/02/14 07:56:30 ]: >>799
自分の頭使えよ。ヘルプでMultipleListPlotやそれらしきものを探したか？
多分バージョンが古いんだろうな。
801 名前：１３２人目の素数さん [05/02/14 10:59:53 ]: マニュアル以外に
お勧めの本は無いのでしょうか？

テンプレも無いみたいでしたが
802 名前：名無しさん@Vim%Chalice [05/02/14 11:42:21 ]: >>801
プログラミングMATHEMATICA - バージョン3&4対応
www.amazon.co.jp/exec/obidos/ASIN/4894711613/qid=1108348626/sr=1-8/ref=sr_1_10_8/250-5779356-1145006
内容が言語の根幹にかかわることなので、5でも変わりない。
803 名前：１３２人目の素数さん [05/02/14 19:26:34 ]: 画像をImportして、その画像に対して中心を原点、四隅に対してx軸、y軸の値を設定したいんですが、
どうしたらいいんでしょうか。
804 名前：１３２人目の素数さん mailto:sage [05/02/15 20:55:20 ]: 初心者です．お教えください
フーリエ変換するときに，係数をデフォルト値と変えるために，毎回，FourieParametersを打ち込んでます．
こんなカンジ
FourierTransform[1, t, w, FourierParameters -> {1, -1}]
毎回　FourierParameters -> {1, -1}　を打ち込むのが面倒なんで，どこかに記憶することはできませんか？
記憶させた後は，
FourierTransform[1, t, w]
をやれば，上と同じ結果が出てくれるとうれしいです．
805 名前：１３２人目の素数さん mailto:sage [05/02/15 21:36:23 ]: >>804
できます。SetOptions。
起動時に読み込まれるファイルで設定することもできます。
806 名前：804 mailto:sage [05/02/15 22:27:22 ]: SetOptionsできたー！　>>805さんくす！
おこがましいですが、起動時に読み込まれるファイルで設定する方法って
どうするか教えていただけませんか？おねがいします。
807 名前：１３２人目の素数さん [05/02/17 05:24:12 ]: 1から10の間にある乱数を20個発生させて、
それらのうちで｢3以上8以下｣の範囲にある値の数を数えたいです。
前半は
Table[Random[Real, {1, 10}], {20}]
で出来たのですが、後半の数える部分が上手くいきません。
Countを使うのでしょうが、｢3以上8以下｣が設定できません・・。

どのように設定すればいいのか教えていただけませんか？
よろしくお願いします。
808 名前：１３２人目の素数さん [05/02/17 07:59:31 ]: >>807

risuto2[risuto1_]:=which[n<3,0,n<8,1,true,0];
Print[sum[risuto2]];

how about this?
809 名前：１３２人目の素数さん mailto:sage [05/02/17 15:52:01 ]: >>807
Countを使うなら:
Count[Table[Random[Real, {1, 10}], {20}], n_?(3 <= # <= 8 &)]

Length@Select[...]も悪くはない。
810 名前：１３２人目の素数さん mailto:sage [05/02/17 15:53:16 ]: nは要らなかった。_?(...)
811 名前：１３２人目の素数さん mailto:sage [05/02/17 16:19:54 ]: >>808
risutoはねぇだろ
812 名前：807 mailto:sage [05/02/17 17:51:15 ]: >>808-811
ありがとうございました。809さんので無事にできました。
813 名前：１３２人目の素数さん mailto:sage [05/02/17 19:52:30 ]: >>811

榊原さんの教科書風にしたわ
814 名前：１３２人目の素数さん mailto:sage [05/02/18 12:56:39 ]: 軽蔑すべき輩だと思う
815 名前：１３２人目の素数さん [05/02/18 19:15:37 ]: わざとだろ(w
突っ込まれるの待ってたんじゃね?
816 名前：１３２人目の素数さん mailto:sage [05/02/19 16:26:06 ]: 質問なんですが
ttp://documents.wolfram.com/v5/Built-inFunctions/NumericalComputation/EquationSolving/FurtherExamples/NBMLImages/NDSolve/NDSolve_79.ja.gif
この画像の3行目の頭にある
u(1,0)[1,t]
はInputFormではどう書くのですか？
817 名前：１３２人目の素数さん mailto:sage [05/02/19 16:47:28 ]: Derivative[1, 0][u][1, t]
でした。スマソorz..
818 名前：１３２人目の素数さん [05/02/19 20:18:09 ]: Mathematicaからまともな表を表示する方法はありますか？
ちなみにLinkForExcelは可変長リストを返す関数に関して糞なのでダメ。
TeX系はインタラクティブじゃないのでイヤです。
819 名前：１３２人目の素数さん [05/02/19 21:42:06 ]: Do[Plot…],{i,0,100}]
って感じでパラメーター変えながら100枚ほどプロットしたのですが、
それぞれのグラフのタイトルとして i=1 とかって表示させるにはどうしたらいいですか？
820 名前：１３２人目の素数さん mailto:sage [05/02/19 22:46:20 ]: >>818
これじゃだめかい？だめだろうねぇ。
t = Table[Subscript[a, i, j], {i, 5}, {j, 5}];
DisplayForm[FrameBox[GridBox[t, RowLines -> True, ColumnLines -> True]]]

>>819
Table[Plot[Sin[i x], {x, 0, 2π}, PlotLabel -> "i=" <> ToString[i]], {i, 3}]
821 名前：１３２人目の素数さん mailto:sage [05/02/19 23:43:40 ]: >>818
Javaで自作。

Link for Excelが標準的な解なわけだから、
> 可変長リストを返す関数に関して糞
とやらをなんとかして使うのがいいと思うが。
822 名前：819 mailto:sage [05/02/19 23:43:50 ]: >>820
多謝
823 名前：１３２人目の素数さん mailto:sage [05/02/20 08:58:22 ]: こんど、PC新しくするんだけど、Mathematicaには
Pentium4とAthlon64のどちらが向いていますか？
今のところ考えているのはPen4 3.4GHz(EM64T)か
Athlon64 3500+くらいのやつです。
824 名前：１３２人目の素数さん mailto:sage [05/02/20 11:08:39 ]: Linuxなら64bitのMathematicaがあるよ。Windowsは知らん。
825 名前：818 mailto:sage [05/02/20 18:08:45 ]: >>820-821
ありがとうございます。
とりあえずはLinkForExcel+Excel上のグチャグチャで凌ごうと思います。
826 名前：１３２人目の素数さん mailto:sage [05/02/21 10:41:09 ]: >>824
なるほど64bit版ですか。今のところ32bit Windows版を考えてるんですが、安定してるなら
64bit Linux版に乗り換えたいところですね。

どなたかPen4とAthlon64の両方で使っているという方がいましたら、感想などよろしくお願いします。
827 名前：１３２人目の素数さん [05/02/22 12:56:28 ]: << Statistics'Continuous Distributions'

で入れると、どうやっても

Get::noopen: Statistics'Continuousを開くことができません．

でエラー返されちゃうんですが。
バージョンは4です。
828 名前：１３２人目の素数さん mailto:sage [05/02/22 13:00:44 ]: >>827
<<Statistics`ContinuousDistributions`

じゃないの？　ContinuousとDistributionsの間に空白が入るのはおかしい。
829 名前：827 mailto:sage [05/02/22 15:01:59 ]: >>828
解決しました。
「'」と「`」が間違っていました。
ありがとうございました。
830 名前：１３２人目の素数さん mailto:sage [05/02/22 16:25:36 ]: >>829
バックダッシュは間違えやすいよね
831 名前：１３２人目の素数さん mailto:sage [05/02/22 19:15:36 ]: >>830
バッククウォートと言います。

ダッシュってのは横線。'をダッシュというのは日本の中高校数学だけ。
832 名前：１３２人目の素数さん mailto:sage [05/02/23 02:32:07 ]: ちなみに'の正式名称は"prime"ね。
833 名前：１３２人目の素数さん [05/02/27 19:33:25 ]: 上の方と同じような質問なのですが

<< Graphics`ImplicitPlot`
ImplicitPlot[x y - y^2 == 3, {x, -2, 2}]

と打つと
Set::"write": \
"\!$Options[ImplicitPlot]$のタグ\!$ImplicitPlot$\
はプロテクトされています．"

とでるのですが何がいけなかったでしょうか･･
834 名前：１３２人目の素数さん mailto:sage [05/02/27 22:21:26 ]: >>833
俺んとこではでないぞ@Mathematica5
835 名前：１３２人目の素数さん mailto:sage [05/02/28 22:28:13 ]: Mathematica link for Excel使っている人いる？
836 名前：１３２人目の素数さん mailto:sage [05/03/02 21:39:17 ]: すまへん。LOG[1+a x]にフェルミ分布関数の微分をかけたような
被積分関数を０から無限大まで数値積分したい。
そのとき、Methodはどれを選ぶと早いのでしょうか？実際は
もっと複雑な被積分関数なので、まあまあのせいどで
現実的な時間に終わってほしい。台形？やっぱり？
リカーションとか分割とかナイスな設定あったら教えて。
837 名前：１３２人目の素数さん mailto:sage [05/03/03 05:56:09 ]: >>835
持ってるけど使ってないぞ。
838 名前：１３２人目の素数さん [05/03/03 11:00:48 ]: age
839 名前：１３２人目の素数さん [05/03/03 15:32:07 ]: やっぱり計算が遅い。対数関数と指数関数の積をインテグランドに
したものはやっぱり計算が遅い。前者の対数関数を
ｘの多項式で近似し、後者をガウス関数で近似すればええのんかな？
今のところ、積分区間をゼロから無限大まで引っ張っているけど
どうせ重要な積分区間は限られているし、その狭い区間だけ
計算させるようにすれば格段に短縮されるのかな？もう本当
つらい。
840 名前：１３２人目の素数さん [05/03/03 16:59:32 ]: Do[式, {i, min_i, max_i, step_i}]
→ step_i おきに i を min_i から max_i まで変えながら式を繰り返す．

てなこっちゃけど、結果を対数表示させたいので、min_i*step_i^n
ってループさせる方法は？
841 名前：１３２人目の素数さん mailto:sage [05/03/03 19:46:19 ]: >>840
式をそう書いたら言いやん。
842 名前：１３２人目の素数さん [05/03/11 15:57:38 ]: ＞＞376

全く下らん問題やな。数値積分にして、メソッドのところを
モンテカルロにしたら？次にインテグレートで出力結果をえる。
そして二つを比較。
843 名前：１３２人目の素数さん mailto:sage [05/03/12 01:59:26 ]: >>842
230日前のレスにレスする君に脱帽。
844 名前：１３２人目の素数さん mailto:sage [05/03/12 15:33:27 ]: 5.1.1のうpでーたが来たよう
845 名前：１３２人目の素数さん [05/03/13 07:05:41 ]: Mathematicaってどこの数学科でも授業で使い方教えてくれるんですか？
846 名前：１３２人目の素数さん [05/03/18 22:57:03 ]: 円周率の平方根は有理数だそうです。

In[1] := Limit[x - Sqrt[Pi],
x -> 19151686109235253480650512155425821931577/
10805181810206786582137323831443648757793]
Out[1] = 0
847 名前：１３２人目の素数さん mailto:sage [05/03/18 23:26:28 ]: >>846
Block[
{$MaxExtraPrecision=100},
Limit[x-Sqrt[Pi], x-> 19151686109235253480650512155425821931577/
10805181810206786582137323831443648757793]
]
848 名前：１３２人目の素数さん mailto:sage [05/03/19 03:48:57 ]: >>846
4.2だとそうはならんな。
849 名前：１３２人目の素数さん mailto:sage [05/03/19 17:04:04 ]: 昨年卒業したので、Ｓｔｕｄｅｎｔ版からＰｒｏ版へのアップグレード
期限がもうすぐなんだけど、卒業証明書って必須なの？実家に置いてきちゃったよ・・・。
850 名前：１３２人目の素数さん mailto:sage [05/03/19 23:37:44 ]: >>849
バーシティーウェーブに問い合わせしたら？
851 名前：１３２人目の素数さん [2005/03/24(木) 17:06:51 ]: 保守上げ
852 名前：１３２人目の素数さん mailto:sage [2005/03/24(木) 22:16:56 ]: >>849
後輩か教授に金渡して買ってもらえ。
853 名前：１３２人目の素数さん [2005/03/28(月) 07:56:51 ]: interpolation function
良い感じ。
しかし
NIntegrateInterpolatingFunction　意味不明
どうやったら積分収束がはやくなるの
854 名前：１３２人目の素数さん [2005/03/29(火) 03:50:46 ]: MATLAB とかMAPLEをサイトライセンスで入れていると、
MATHEMATICAのサイトライセンスの料金を高くするというような
公正な競争を阻害するライセンスの販売システムになっている。
大学ではみんなボイコットして、MAPLEに乗り換えよう。
855 名前：１３２人目の素数さん mailto:sage [2005/03/30(水) 22:12:03 ]: Mathematicaが直感的に操作できて（ﾟдﾟ）ｳﾏｰ　
856 名前：１３２人目の素数さん mailto:sage [2005/03/31(木) 01:37:06 ]: Mathematica 高すぎ
857 名前：１３２人目の素数さん mailto:sage [2005/03/31(木) 14:09:32 ]: Export[ , ,"GIF"]
が動きません。
こんなエラーが出ます。

Export::fmterr: 無効なPNM形式です

"GIF" 以外はJPEGでもBMPでも動くようです。

documents.wolfram.com/v4-ja/RefGuide/Export.html

一体なにが問題なのか分るかたお教えください。
858 名前：つづき mailto:sage [2005/03/31(木) 14:14:15 ]: documents.wolfram.com/mathematica/functions/Export

ここを見ても、正しいことをしているようです。
バージョン5.1です
859 名前：１３２人目の素数さん [皇紀2665/04/01(金) 14:51:23 ]: >>818

昔作った関数でこういうのがあった。

TF = TableForm[PaddedForm[#, {6, 5}] & /@ # & /@ #,
Alignments -> Right, TableSpacing -> {0, 0}] &;

Table[…]//TF

みたいな使い方をしてる。

もう忘れたから、
Mathematicaに貼り付けて
「選択項目のヘルプ」でしらべてちょ。
860 名前：１３２人目の素数さん [2005/04/02(土) 21:13:23 ]: 二分木なデータ構造をMathematicaでやろうとしたらどう実装すればいいですか？
Lisp風なリスト？
861 名前：１３２人目の素数さん mailto:sage [2005/04/02(土) 23:30:44 ]: Mathematica自体がLisp風な訳だが。
862 名前：１３２人目の素数さん [2005/04/03(日) 11:01:02 ]: >>861
でも、リストの深いところを破壊しようとすると全部書き換える事になって遅くなりませんか？
863 名前：862 mailto:sage [2005/04/03(日) 11:24:11 ]: つうかリストの破壊が基本的にできない気がしてきました。書き換えしかない。
864 名前：１３２人目の素数さん mailto:sage [2005/04/03(日) 12:19:28 ]: >>862
二分木のデータを変更すると、そこに至るための上の部分の木も再生成されるが、log(n)だ。
だから、普通Schemeで二分木使う時はset!でデータ書き換えたりしないと思うのだが。

Mathematicaでは「書き換え」はシンボルに対してしかできないだろう。
一部を書き換えると全部が「コピー」または「再生成」される。
どうしても書き換えたいなら、シンボルによる二分木を作っておいて、あとでシンボルに値を
アサインするようにすればいいんじゃないか。
865 名前：862 mailto:sage [2005/04/03(日) 22:19:37 ]: >>864
そうですよね、「コピー」または「再生成」しかできない。
（ちょっとこっちな言語の用語に詳しくないもんで、すいません。）
Unique[]ってのがあるみたいなんで、それでノードを乱打して繋げてみます。
ありがとうございます。

#今日リストの大コピー合戦で実装してみたら、
#Length[Flatten[リスト]] = 10000あたりでもう超ウルトラ遅い。
#ノード数は200ぐらいなのにorz
866 名前：１３２人目の素数さん mailto:sage [2005/04/06(水) 23:34:24 ]: 教えてください
Simplify[Abs[a + b*I], {a ∈ Reals, b ∈ Reals}]
とすると，
Abs[a + I b]
という答えが返ってきて困っています．aもbも実数と言ってるんだから，もっと簡単になるだろ！
Sqrt[a^2 + b^2]
という解を得るには，どうしたらよいのでしょうか？
867 名前：１３２人目の素数さん mailto:sage [2005/04/07(木) 00:35:10 ]: >>866
機械にとっては、次のどちらが簡単だと考える方が自然かね？
Abs[Plus[a, Times[Complex[0, 1], b]]]
Power[Plus[Power[a, 2], Power[b, 2]], Rational[1, 2]]
868 名前：866 mailto:sage [2005/04/07(木) 01:36:28 ]: 自己解決しました．867さん，ご意見ありがとう
ComplexExpand[Abs[a + b*I]]
でやってくれました
869 名前：１３２人目の素数さん [2005/04/08(金) 13:02:41 ]: >>868
> 自己解決

自己解凍みたいw
870 名前：博士 ◆OBpuRHyGa. mailto:もっちーＡＧＥ [2005/04/08(金) 14:58:47 ]: lim f(x)=sin∫f(x)となるようなf(x)をヒルベルト空間を用いて
x→∞
求めよ。バナハ空間でもよい。
871 名前：１３２人目の素数さん [2005/04/08(金) 15:17:58 ]: >>870
積分の範囲は？
872 名前：博士 ◆OBpuRHyGa. mailto:もっちーＡＧＥ [2005/04/08(金) 16:12:06 ]: 0→∞
873 名前：１３２人目の素数さん mailto:sage [2005/04/08(金) 20:41:10 ]: >>870

なんでヒルベルト空間なんているんだよ。
lim[n→∞]f(x)=0
∫[0,∞]f(x)dx=0
にすればいいだけじゃん。
たとえば
f[x]=(-1)^[x]/[1+x/2]
874 名前：873 mailto:sage [2005/04/08(金) 20:45:48 ]: >>873
lim[n→∞]f(x)=0
↓
lim[x→∞]f(x)=0
875 名前：１３２人目の素数さん [2005/04/14(木) 15:41:02 ]: >>868

MyRe[H_] :=
1/2*(ComplexExpand[H, TargetFunctions -> {Re, Im}] +
ComplexExpand[Conjugate[H], TargetFunctions -> {Re, Im}])
876 名前：１３２人目の素数さん [2005/04/24(日) 10:50:44 ]: 保守のため、上昇
877 名前：１３２人目の素数さん [2005/04/28(木) 00:33:04 ]: Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals]で
Conjugate[a + I b]となり, a-I bにならないのですがどうしてでしょう？

ちなみに、
Simplify[Conjugate[I b], a ∈ Reals && b ∈ Reals]は
-i b　になるし、
Simplify[Conjugate[a+I ], a ∈ Reals && b ∈ Reals]は
-i+a　になります。
878 名前：１３２人目の素数さん mailto:sage [2005/04/28(木) 19:14:20 ]: >>877
Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals]でa - I bになるが。
@Mathematica 5.0
879 名前：１３２人目の素数さん [2005/04/28(木) 19:55:03 ]: >>878

Mathematica4.0.1 mac版
では

Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals]は

Conjugate[a + I b]

です。878さんのはwindows版ですか？
880 名前：878 mailto:sage [2005/04/29(金) 19:29:24 ]: >>879
Win版だけど、バージョンの違いだと思われ。OS変えるだけで
わざわざカーネルのソース書き換えないだろう。
881 名前：１３２人目の素数さん mailto:sage [2005/05/01(日) 10:16:29 ]: Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals, ComplexityFunction -> LeafCount]

V4系にはConjugate[..]のほうが単純に見えるんだろう。
882 名前：１３２人目の素数さん [2005/05/03(火) 10:00:50 ]: age
883 名前：１３２人目の素数さん [2005/05/10(火) 20:25:40 ]: hos
884 名前：１３２人目の素数さん [2005/05/11(水) 19:13:01 ]: ちょっと質問なのでが
Mathematicaを使ってアイパターンを描画させたいのですがMathematicaでできるでしょうか？
gnuplot使っていたのですがアイパターンを描けずMathematicaに移行したいと思っています
885 名前：１３２人目の素数さん mailto:sage [2005/05/11(水) 21:38:09 ]: gnuplotでできなかったのは、誰にも聞かなかったからで、
Mathematicaでできたとしたら、それは誰かに教えてもらえたからでしょうね。
886 名前：１３２人目の素数さん [2005/05/12(木) 17:56:57 ]: プログラム板でレスなかったのでこちらで質問させてください。

CUIで使ってるときに、1行に表示される文字数を変更する方法はありませんか？
例えば、

In[1]:= N[Pi,100]
Out[1]= 3.1415926535897932384626433832795028841971693993751058209749445923078\
> 16406286208998628034825342117068

としたときに、結果を改行せずに表示してほしいのですが。
887 名前：１３２人目の素数さん mailto:sage [2005/05/12(木) 18:51:26 ]: Mathematica 5.01のMac版とSolaris版で試したのですが、
In[22]:= 1. 10^(-309) +0.0
Out[22]= 0.
上のようになるのはバグでなく機械精度の問題ですか？

また、この結果を
1. 10^(-309)
となるようにする方法ってありますか？
888 名前：887 mailto:sage [2005/05/12(木) 19:18:22 ]: すいません
In[22]:= 1. 10^(-309) +0.0
Out[22]= 0.
こうなるのは、有効桁数のこと考えたら当たり前ですね。逝ってきまつ
889 名前：１３２人目の素数さん mailto:sage [2005/05/12(木) 21:03:13 ]: >>886
無理
890 名前：１３２人目の素数さん [2005/05/15(日) 10:27:27 ]: さっきランダムウォークの
コードを書いた。はろうか
891 名前：１３２人目の素数さん [2005/05/15(日) 10:57:05 ]: 今起きたばかりで
　コンドームと一瞬読めてしまった俺ｶﾞｲﾙ
892 名前：１３２人目の素数さん mailto:sage [2005/05/15(日) 16:55:26 ]: 迷ってるのなら、貼れ
893 名前：１３２人目の素数さん [2005/05/19(木) 22:01:32 ]: >>892

今日また、コードを書いた。
894 名前：１３２人目の素数さん [2005/05/20(金) 18:44:08 ]: Ver.5.1　９ラック版ある？
895 名前：１３２人目の素数さん mailto:sage [2005/05/20(金) 20:06:58 ]: 　　　　　　〃〃　ｶﾞｯ
　　　＿__〃　 ∧＿∧
　　　|　　|　　（・∀・ ,,）
　　 .=|　　|==○○=　ヽ
　　从|＿_|∩　ヽ　　､(⌒)っ
　（　｀Д´）/　　ヽ＿）
　　　　　　/ ←>>894
896 名前：１３２人目の素数さん mailto:sage [2005/05/21(土) 21:57:03 ]: >>894
刑務所行って探して来い
897 名前：１３２人目の素数さん [2005/05/25(水) 18:10:48 ]: おまいらMathematicaMark 2004の結果教えてくり
やり方は
<<Utilities`Benchmark`
BenchmarkReport[]

おいらはﾁﾝｺﾊﾟｯﾄﾞ T42 PenM 1.7GHz 512MB Mathematica5.1.1で1.59
898 名前：１３２人目の素数さん mailto:sage [2005/05/31(火) 12:33:44 ]: ３点O(0,0,0)、A(1,0,0)、B(0,1,1)を通る平面上に、
Aを中心とする半径１の円をｙ軸の周りに１回転してできる回転体を
Mathematicaを使って書きたいんですけど
どんなコマンドにしたらいいですか？
899 名前：ChaosicSoul ◆/yaJbLAHGw [2005/05/31(火) 12:50:54 ]: Re:>>898
パラメータ表示で曲面を描くコマンドがあったはずだ。一個のパラメータで円を作って、もう一個のパラメータは回転用に使えばいいだろう。
円のパラメータ表示は、(1+cos(t),√(2)sin(t)/2,√(2)sin(t)/2)でいいだろう。
これに回転用の行列を掛ければ完成だ。パラメータ表示を縦ベクトルにして、
((cos(s),0,-sin(s)),(0,1,0),(sin(s),0,cos(s)))^T を左から掛ければいいだろう。
900 名前：ChaosicSoul ◆/yaJbLAHGw [2005/05/31(火) 12:51:19 ]: ちなみに正しいスペルは chaotic のようで。
901 名前：８９８ mailto:sage [2005/05/31(火) 12:55:36 ]: コピペでMathematicaに貼り付けたらできちゃうように
コマンドを作ってもらってもいいでしょうか…
902 名前：ChaosicSoul ◆/yaJbLAHGw [2005/05/31(火) 13:02:29 ]: Re:>>901 今Mathematicaが手元に無いから何とかして自分で作ってくれ。
903 名前：８９８ mailto:sage [2005/05/31(火) 13:07:28 ]: >902
ありがとうございました。

自分はあんまりMathematicaは使いこなせないんで、
ChaosicSoulさんがMathematicaを使えるようになったり、
他の方でも、できるようならコマンドを書いていただけるとうれしいです。
皆さんお願いします。
904 名前：１３２人目の素数さん mailto:sage [2005/05/31(火) 15:26:04 ]: そこまで人に頼るのかよ
905 名前：Zeno ◆5nZQbNmQPs [2005/05/31(火) 17:56:11 ]: Reply-to:>>903
　　　　　　　　　　　　　君は将来性がない人間だな。
　　　　　　　　　そんな簡単なグラフは自分で書くようにしろ。

________________________________________________________
セオドア・カジンスキー「産業社会とその未来」 (aka. マニフェスト) の抜粋
www.cnn.com/US/9509/unabomber/09-19/am/excerpt.html
906 名前：１３２人目の素数さん mailto:sage [2005/05/31(火) 18:25:23 ]: >>903
何故かこっち(Methematica 5.0)では途中で固まるが、これであってるはず
ParametricPlot3D[Transpose[{{Cos[s], 0, -Sin[s]}, {0, 1, 0}, {Sin[s], 0, Cos[s]}}].{1 + Cos[t], Sin[t]/Sqrt[2], Sin[t]/Sqrt[2]}, {t, 0, 2π}, {s, 0, 2π}];
907 名前：８９８ mailto:sage [2005/05/31(火) 21:37:46 ]: >>905
英語ばっかりだ～
読めん…

>>906
ありがとー！
やってみます。
908 名前：１３２人目の素数さん mailto:sage [2005/06/01(水) 00:12:28 ]: 学校の宿題を代わりにやってくれるスレはここですか
909 名前：ねえねえ [2005/06/01(水) 09:57:48 ]: science3.2ch.net/test/read.cgi/math/1051275230/l50
>>297
誰か私にちからを…
910 名前：ChaosicSoul ◆/yaJbLAHGw [2005/06/01(水) 10:37:09 ]: Re:>>909 Cのファイル入出力でバイナリモードを使えばできる。何が改行コードかに注意。
911 名前：ねえねえ [2005/06/01(水) 11:14:48 ]: 承知!!
912 名前：ねえねえ [2005/06/01(水) 11:27:30 ]: >>910
ありがとー
出来そうです。
913 名前：１３２人目の素数さん [2005/06/02(木) 13:40:42 ]: u, v を実数, k を定数としたとき、
-5u^{4}v + 10u^{2}v^{3}-v^{5} = k
のグラフをUV平面に書きたいのですが，
方法が分かりません。
とりあえず k=1 の時を書こうとしましたが、
Plot を使うには v について解く必要がありそうですし，
Solve とか試してみましたがエラーで全然動きません。
どなたか教えて下さい。
914 名前：913 mailto:sage [2005/06/02(木) 14:31:55 ]: Plot[{{Root[1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 1]}, {Root[
1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 2]}, {Root[
1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 3]}, {Root[
1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 4]}, {Root[
1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 5]}}, {u, -10, 10}]

とかしてみたらそれらしいのが出ました(k=1)。
何か違う気もしますが，もう少しいじってみます。
915 名前：１３２人目の素数さん mailto:sage [2005/06/02(木) 15:06:21 ]: >>913
ImplicitPlotで。

UVP[k_] :=
ImplicitPlot[ - 5 u^4 v + 10 u^2 v^3 - v^5 == k, {u, -10, 10}, {v, -10, 10}]
UVP[1]
916 名前：915 mailto:sage [2005/06/02(木) 15:07:04 ]: <<Graphics`ImplicitPlot`
ね。
917 名前：１３２人目の素数さん mailto:sage [2005/06/02(木) 22:11:37 ]: >>915
うまいね、UVPという関数を作る感じか。
918 名前：918 [2005/06/02(木) 23:00:40 ]: 9=1+8
9-1=8
919 名前：919 mailto:sage [2005/06/02(木) 23:02:02 ]: ｸｲｯｸ　ワイパー
920 名前：１３２人目の素数さん [2005/06/07(火) 11:13:05 ]: 今日、大学生協行ったら5.2へのフリーアップグレードつきのパッケージが出てた。
俺が2週間前に買ったときは無かったのに
　　　　　　　　　　　　　　　　　┌─┐
　　　　　　　　　　　　　　　　　｜も.｜
　　　　　　　　　　　　　　　　　｜う　|
　　　　　　　　　　　　　　　　　│来│
　　　　　　　　　　　　　　　　　│ね│
　　　　　　　　　　　　　　　　　│え .|
　　　　　　　　　　　　　　　　　│よ .|
　　　　　　　　　　　　　　　　　│ !!.│
　　　　　　ｼｮﾎﾞｰﾝ　　　ﾊﾞｶ　└─┤　　　　ﾌﾟﾝﾌﾟﾝ
　　　　　ヽ(´･ω･`) ﾉヽ(`Д´)ﾉ　　(`Д´)ﾉ　　　（　`Д）
　　　　　　|￣￣￣|─|￣￣￣|─|￣￣￣|─□（ヽ┐U
　　～～　￣◎￣　　.￣◎￣　　￣◎￣　　 ◎－>┘◎
921 名前：１３２人目の素数さん [2005/06/07(火) 16:51:27 ]: ここは大学で出た宿題がわかんないんだけど聞いてもいいのかな?
922 名前：１３２人目の素数さん mailto:sage [2005/06/07(火) 17:21:59 ]: 5.2ってまだ出てないよね
いつ出るんだろ
923 名前：１３２人目の素数さん mailto:sage [2005/06/07(火) 18:05:57 ]: >>921
いいんジャマイカ？
924 名前：９２１ [2005/06/07(火) 23:07:59 ]: やったね＾＾ありがと
定点A（a,0）と楕円(x^2/2)+y^2=1上の動点Pとの距離APの最小値ｂを求めよ
また、aを横軸、ｂを縦軸にとって、aとｂの関係をグラフに示せ
っていう問題なんですが、条件を置いてからのそこから先の展開ができません。
925 名前：１３２人目の素数さん mailto:sage [2005/06/08(水) 02:15:02 ]: スレ違いだから簡単に
Pを(√2 cosθ,sinθ)と置いて距離APの2乗を導き、微分=0を調べればいいんでは
つか高校レベルじゃない？arccosとか出てくるから大学なのか
926 名前：１３２人目の素数さん mailto:sage [2005/06/09(木) 18:57:38 ]: 5.1.1がリリースされたみたい。　5.1自体つい半年前に出たばかり。
ただいま、某サイトからダウソ中。
927 名前：１３２人目の素数さん mailto:sage [2005/06/09(木) 20:22:38 ]: 5.1.1は前からあるよ。
つかワレ厨は氏ね！
928 名前：１３２人目の素数さん mailto:sage [2005/06/09(木) 20:47:05 ]: ウルフラムのホムペからも金払えばダウソ出来るよ。
929 名前：１３２人目の素数さん mailto:sage [2005/06/10(金) 02:35:20 ]: 既知
930 名前：１３２人目の素数さん [2005/06/15(水) 18:34:10 ]: m
931 名前：１３２人目の素数さん [2005/06/17(金) 02:36:27 ]: LiveGraphics3Dで3Dグラフィクスを作ってるんだが、Mathematicaでちゃんと
PlotRange使ってプロット範囲指定してるのにWriteLiveFormで変換してHTMLで見ると
PlotRangeで指定したことがまるっきり無視されるんだが、どうすればいいんだ？
932 名前：１３２人目の素数さん mailto:sage [2005/06/17(金) 08:12:56 ]: jisaku.u-satellite.net/uploader/files/869.gif
933 名前：１３２人目の素数さん [2005/06/17(金) 12:26:16 ]: ｲﾝﾃﾙｵﾜｯﾄﾙ
934 名前：１３２人目の素数さん mailto:sage [2005/06/17(金) 20:17:49 ]: ｿﾝﾅｺﾄﾜｶｯﾄﾙ
935 名前：１３２人目の素数さん [2005/06/18(土) 12:28:34 ]: matlabにしてもmathematicaにしても、ここまで差が出るとは。
理系の研究室にopteronが売れまくるのも当然だな
936 名前：１３２人目の素数さん mailto:sage [2005/06/18(土) 16:28:56 ]: すいません，質問なんですが，
ListPlot[{{x1,y1},{x2,y2},…]のように，
3次元のデータ点{x1,y1,z1},{x2,y2,z2},…を
3Dプロットする方法はないですか？
937 名前：１３２人目の素数さん mailto:sage [2005/06/18(土) 22:21:02 ]: >>936
Graphics`Graphics3D`のScatterPlot3D
938 名前：１３２人目の素数さん [2005/06/21(火) 20:01:05 ]: x≧0かつy≧0かつx+y≦3のとき
z=x-2-2xy+2y^2-2x-2y+6　のグラフを書きたいんですがどうしたらよいでしょうか

もう一つ同じく不等式系とグラフで 2x^2+3xy+2y^2≦7のとき正の定数aを用いて
z=(x+a)(y+a)
のグラフも教えていただけると幸いです

お願いします
939 名前：939 mailto:sage [2005/06/21(火) 20:24:37 ]: 939 ｿｹｯﾄ
940 名前：１３２人目の素数さん mailto:sage [2005/06/22(水) 02:31:11 ]: >>938
その不等式くらい解け。
Plot3D[x - 2 - 2x y + 2y^2 - 2x - 2y + 6, {x, 0, 3}, {y, 0, 3}];
941 名前：１３２人目の素数さん [2005/06/22(水) 19:54:15 ]: 938

2x^2の不等式の奴はむりだとおもわれ
942 名前：１３２人目の素数さん mailto:sage [2005/06/22(水) 22:02:53 ]: >>938
region[x_, y_] /; 2 x^2 + 3 x y + 2 y^2 <= 7 = 1
region[x_, y_] = Infinity
Plot3D[region[x, y] (x + a) (y + a) /. a -> 3, {x, -3, 3}, {y, -3 , 3},
PlotPoints -> 80]

手抜きだが、エラーメッセージは無視してくれ。
943 名前：１３２人目の素数さん [2005/06/23(木) 16:04:48 ]: mathematicaって途中の式を表示する機能ってありますか？
944 名前：１３２人目の素数さん mailto:sage [2005/06/23(木) 16:55:07 ]: >>943
Traceは？
945 名前：１３２人目の素数さん mailto:sage [2005/06/23(木) 19:31:21 ]: >>944
これなら数学の勉強にも使えそうでつね、ありがとん。
946 名前：１３２人目の素数さん [2005/06/29(水) 21:50:40 ]: ３次元のデータ(x,y点での値)
を2次元平面状に色であらわしたいのですが、どうすればいいでしょうか？
たとえば、
f(x,y) = 1000 なら赤色
f(x,y) = 100 なら青　
みたいに大きいほど赤　すくないほど　青　みたいに表したいです。
947 名前：１３２人目の素数さん [2005/06/29(水) 22:09:29 ]: 安定な特異点と不安定な特異点を各１つ以上もつ非線形なベクトル場を考え、Mathematica でプロットせよ。
ってどうやればいいんですか？
948 名前：１３２人目の素数さん mailto:sage [2005/06/29(水) 22:31:50 ]: >>946
ContourPlot と ColorFunction
949 名前：１３２人目の素数さん [2005/06/29(水) 22:50:05 ]: age
950 名前：１３２人目の素数さん mailto:sage [2005/06/29(水) 22:53:15 ]: >>947
Graphics`PlotField3D`
PlotVectorField3D
951 名前：１３２人目の素数さん [2005/06/29(水) 23:25:32 ]: >>950
詳しく教えてもらえたら嬉しいです！
952 名前：１３２人目の素数さん mailto:sage [2005/06/30(木) 10:04:30 ]: >>951
> 安定な特異点と不安定な特異点を各１つ以上もつ非線形なベクトル場
はもう考えてあるの？
953 名前：１３２人目の素数さん [2005/06/30(木) 11:47:01 ]: >>951
それすらわかりません。模範解答みたいなの教えてもらえませんか？
954 名前：１３２人目の素数さん [2005/06/30(木) 11:52:38 ]: ↑>>952でした
955 名前：１３２人目の素数さん [2005/06/30(木) 15:27:59 ]: age
956 名前：１３２人目の素数さん mailto:sage [2005/07/01(金) 01:22:30 ]: V5の日本語版の尻おしえてくれぇぇ
957 名前：１３２人目の素数さん mailto:sage [2005/07/01(金) 02:21:51 ]: 　　　　　　　　　　　　　＿＿＿＿
　　　　　　　　　　　.／, - ､, - ､　　￣ヽ
　　　　　　　　　　./-┤　｡|｡　 |――-､　ヽ
　　　　　　　　　|　ヽ｀- ○- ´ ／　　ヽ　 |
　　　　　　　　　|　-　　 |　　　―　　 |　 |
　　　　　　　　　|　 ´　　|　　　｀ヽ　 .　|　 |ヽ
　　∩　　　　　　人｀､ _ 　|　　　＿.- ´ |　.|　　＼
　　|　　⌒ヽ　／　＼　￣￣￣　　　ノノ　　　＼
　　|　　　 |´　　　　|￣―--―― ´ヽ　　＿　　／⌒＼
　　＼_ 　 _/-―――.|　（Ｔ）　　　　　　｀l 　　 Τ（　　　）
　　　￣　　　　　　 |　　￣　　　　　　　}　　　 |　＼＿／
　　　　　　　　　　　　 | ､--―　￣|　　　/
ＦＵＣＫ　ＹＯＵぶち殺すぞ・・のび太・・・・！お前は、大きく見誤っている。この世の実態が見えていない。
まるで３歳か４歳の幼児のように､この世を自分中心・・
求めれば､回りは右往左往して世話を焼いてくれる､そんなふうに、まだ考えてやがるんだ・・。臆面もなく・・・・！

甘えを捨てろ。お前の甘え・・その最もたるは、毎日のように叫んでいる､その、「ドラえもぉ～ん！」だ。
泣きつけば道具を出してくれるのが当たり前か・・・・？なぜそんなふうに考える・・・・？バカがっ・・・・！
とんでもない誤解だ。世間というものはとどのつまり、誰も他人の助けになど応えたりはしない。
飢餓で苦しむアフリカの子供たち・・戦争で住む家を失った者たち・・
食べ物を持て余し､馬鹿げたことに金をつぎこむ日本人は、彼らの求めに応えてやったか？応えちゃいないだろうが・・！
これは身近にないからだとか、そういうことじゃあない。目前にあってもそうなのだ。
何か得られるものでもない限り､他人を救ったりしない。それが基本だ。
その基本をはき違えているから、わざわざ２２世紀からオレが来る羽目になったんだ・・・・！！

無論中には、助けるものもいる。しかしそれは自分にとって都合がいいからであって・・つまりは・・のせられてるってことだ。
なぜ・・それがわからない・・？なぜ・・それに気付かない・・・・？
958 名前：１３２人目の素数さん mailto:sage [2005/07/01(金) 02:31:03 ]: >>956
素直に、英語版使いな。。　偉そうな顔をしている自称苦楽屋がいるけど、連中に
そんな腕はない。
959 名前：１３２人目の素数さん mailto:sage [2005/07/01(金) 08:56:34 ]: 　　　　　　　　_....._{{ 〃
　　　　　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　｀ヽ',　 ',-<
　　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）)
　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ.
　　　',　',　|Vｧ=､ﾞ､｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　ﾞ',
　　　ヽ､', l:!Kノ}. 　　 f:_.)iﾞi: ﾘ !　l　ル　
　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉﾞ ,ﾘ l ﾘ'´ 　　　
.　　　　 ',|!!､　　r‐┐ 　｀ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣　
　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　＜　 >>956氏ね、ヴォケ！
　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　　　　＼＿＿＿＿＿＿＿＿＿＿
　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､　
　　/＼＼　　　　',　　 }　　 //｀ヽ
960 名前：１３２人目の素数さん mailto:sage [2005/07/01(金) 23:53:14 ]: 英語版のシリアルを入力し、英語になってしまったんで、
シリアルを変えたいんですが、どうやったらできるんですか？
961 名前：１３２人目の素数さん mailto:sage [2005/07/02(土) 00:19:16 ]: 5.1の場合、

C:\Documents and Settings\All Users\Application Data\Mathematica\Licensing\mathpass

を消す。　それ以前の場合は、

　Wolfram Research\Configuration\Licensing\mathpass

を消す。　mathpassってファイルにパスワードが書いてあるから、そのファイルを関係有りそうな
フォルダから検索すれば見つかる。
962 名前：１３２人目の素数さん mailto:sage [2005/07/02(土) 01:09:19 ]: 1+2
ってかいて、SHIFT+ENTERおしたら、ウィンドウズのエラーが表示されて計算しないんだけど・・・
なんでかわかりませんか？
963 名前：１３２人目の素数さん mailto:sage [2005/07/02(土) 07:45:47 ]: エラーの内容ぐらい書けよ
つか、呪われてるんとちゃうｗ
964 名前：１３２人目の素数さん mailto:sage [2005/07/03(日) 00:05:34 ]: >>962
kernelから壊れてるんじゃないのか？
965 名前：１３２人目の素数さん mailto:sage [2005/07/05(火) 10:34:04 ]: 5.2ﾏﾀﾞｰ?
ﾁﾝﾁﾝAA略
966 名前：１３２人目の素数さん [2005/07/06(水) 19:36:00 ]: age
967 名前：１３２人目の素数さん [2005/07/08(金) 21:43:37 ]: ある複素関数の実数部を取り出して、それをPlot3Dしたいのですが、
うまくいきません。
例えば(x+yI)^3の実数部の関数をxとy、それぞれ0から5までの範囲で
3次元プロットしたい場合、次のようにやってみたのですがエラーが出て
しまいます。

Plot3D[Re[(x+yI)^3],{x,0,5},{y,0,5}]

どのようにすればいいでしょうか？
968 名前：１３２人目の素数さん mailto:sage [2005/07/08(金) 22:06:10 ]: yI != y I。下のをコピペしてみよう。

Plot3D[Re[(x + y I)^3], {x,0,5}, {y,0,5}]
969 名前：１３２人目の素数さん [2005/07/08(金) 22:14:16 ]: >>968
あー、できました！ありがとうございます！
初心者で基本も分かっておらず、すいませんでした。
970 名前：犬笠銀次郎 mailto:http://www6.speednet.ne.jp/～91nj1r0 [2005/07/08(金) 22:16:48 ]: Mathematica は高いらむ～～～～～。
971 名前：１３２人目の素数さん [2005/07/09(土) 23:01:45 ]: 計算して出した関数の逆関数のグラフをプロットさせてみたいとき、
どうすればいいでしょうか？

例えば
Plot[InverseFunction[x^2], {x, 0, 10}]

ではダメでした。これは逆関数が一意に決定しないからですか？
しかし、

Plot[InverseFunction[x], {x, 0, 10}]

でもダメでした。
どうかアドバイスをお願いします。
972 名前：１３２人目の素数さん mailto:sage [2005/07/10(日) 01:59:06 ]: Plot[InverseFunction[# &][x], {x, -3, 3}]
973 名前：１３２人目の素数さん mailto:sage [2005/07/10(日) 03:26:59 ]: >>971
ImplicitPlot
974 名前：１３２人目の素数さん mailto:sage [2005/07/14(木) 19:52:09 ]: ver5.1が出てからだいぶ経ちましたが、ver5.0以前には無い機能でこれはイイという機能はありますか？
975 名前：１３２人目の素数さん [2005/07/14(木) 21:36:02 ]: Ver5.2 発表　下げ
976 名前：１３２人目の素数さん [2005/07/14(木) 21:44:29 ]: A,B,C,D,E,Fの６個の玉に糸を通して首飾り
を作るとき何通りの作り方があるか？
（⇒回転以外に裏返して同じになるもの
も全部同じ首飾りとして考える。）

円順列ではないようなので、
どうやって解けば良いのかわかりません。
教えてください。
977 名前：１３２人目の素数さん mailto:sage [2005/07/14(木) 21:49:50 ]: 高校生に聞け
978 名前：１３２人目の素数さん [2005/07/14(木) 21:50:19 ]: コンビ屋の超難問です
979 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 00:40:41 ]: www.wolfram.com/products/mathematica/newin52/
Athlon 64 X2が欲しくなってきた
980 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 03:12:13 ]: >>976
(6-1)!/2
981 名前：１３２人目の素数さん [2005/07/15(金) 10:39:29 ]: (d^2x)/(dt^2)+μ(x^2-1)(dx/dt)+x=a+εcos(bt)の数値計算プログラム作成は
どうやればいいのでしょう？
982 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 15:01:09 ]: ﾙﾝｹﾞ･ｸｯﾀ!
983 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 16:38:14 ]: 次ｽﾚ
science3.2ch.net/test/read.cgi/math/1121413040/
984 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:17:11 ]
985 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:17:18 ]
986 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:17:39 ]
987 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:17:57 ]
988 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:18:42 ]
989 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:18:45 ]
990 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:19:15 ]
991 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:19:17 ]
992 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:19:37 ]
993 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:19:38 ]
994 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:19:59 ]
995 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:20:00 ]
996 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:20:19 ]
997 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:20:21 ]
998 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:20:39 ]
999 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:20:59 ]
1000 名前：１３２人目の素数さん mailto:sage [2005/07/15(金) 18:21:02 ]: 〓 Mathematica 3 〓
science3.2ch.net/test/read.cgi/math/1121413040/
1001 名前：１００１ [Over 1000 Thread]: このスレッドは１０００を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef