- 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/
- 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 ]
- ここまで投げやりなのも珍しいな？ 低凶？
|
 |
