mathematicaによるシミュレーション教えてください

**名無しさん＠１周年** · NG

mathematicaを使って時間を変数としてもつ関数を
求めたのですが，この関数の例えば0.1秒刻みの値を
出力させるにはどうしたらよいのでしょうか？
できれば，Exelに出力させたいのですが。
ご存知の方，教えていただけませんか。

**めるとも** · NG

マセマティカのメーリングリストを探したほうがいいのでは？
『マセマティカによる数値計算』という本が出てました。

ところでマセマティカでグラフが容易に書けるのに、
またはプログラミングが容易にできるのに、わざわざ
プログラミング言語でルンゲクッタのプログラムを書いて
微分方程式のパスを求めたりする必要があるのでしょうか？
（無いのではないでしょうか？）
また、シミュレーションの雑誌に書いてあることは、
より早く計算できるプログラムについてだと勝手に想像
しているのですが（読んでないから）、遅くてもいいや
と思う人にとって、早く走るプログラムを勉強することは
無駄だと思うのですが、いかがなものでしょうか？
誰か質問に答えていただけないですか？

**名無しさん＠１周年** · NG

Blackboxでやるよりもソース書いてやる方が，
自分で色々修正できるし，後々便利だよ．
やっぱり元が全部分かってるのが一番安心できます．

**とくめー** · NG

名無しさん＠１周年さんのいう通りだと思います。
以前、モンテカルロ法を使ってシミュレーションするのにmathematicaと
Cの２通りで実験してみたのですが、mathematicaの場合はどんなに工夫しても
Cよりは実行時間がかかるようです。

**名無しさん＠１周年** · NG

ところでmathematicaはどこのサイトからダウンロードできるのでしょうか？

**名無しさん＠１周年** · NG

フリーじゃないってば＞マセマティ子

5 · NG

>6
フリーじゃないのは知っています。
で，どこのサイトからダウンロードできるのでしょうか？

**名無しさん＠１周年** · NG

＞５
ネタかと思ってた、、、
www.wolfram.comへどうぞ。ダウンロードは
サポートしてないと思うけど、、、

5 · NG

>8
www.wolfram.com も知っています。
もう一度だけ伺います。どこのサイトからダウンロードできるのでしょうか？

**とくめー** · NG

mathematicaはフリーなの？
studentバージョンなら確か2万円くらいで売っていたような
気がするけど.

5 · NG

参りました。質問を変えます。
私は普通の会社員です。最も安くmathematicaを購入した場合
価格はいくらぐらいになるのでしょうか？

学割を利用するなら放送大学に入学するという技も知っています。
この場合は入学金に18,000円かかりますが、授業を１つもとらなければ
それ以上の費用はかかりません。
OSはwindows2000 or linux（redhat6.1）です。linux版なら安そうですね。

M · NG

つーか、student version ならば２万円ぐらい。ただし学生onlyです。
大学の生協で売ってますけど、放送大学の場合は？

それ以外だと、代金１０万以上すると思う。
何故mathematicaにこだわるの？

**とくめー** · NG

mathematicaにこだわらなくてもCで書いたほうが計算は
絶対早い。グラフィックスはきれいだけどね。

**名無しさん＠１周年** · NG

なんだ。５はwarez野郎か

**実習生さん** · NG

論文でMathematicaが一番初心者向きとあった。
それができればマトラブ、それができればC等
マセマティカは計算速度が遅いらしい。
これも言語だから習熟するのが大変のようだ。
ちょっとしたことなのにプログラムを組まなくちゃいけない。
単純なことで使いたいなら別のソフトを探したほうがいい。

**実習生さん** · NG

ところでシュミレーションにマセマティカを
使うにはどういう本を勉強すれば言いのでしょうか？

**実習生さん** · NG

『マセマティカによる数値計算』
はマセマティカのコマンドでやっていること
（ブラックボックス）の中身を説明しています

マセマティカを使ったシュミレーションの本ではなかったと思います

**名無しさん＠１周年** · NG

マスマティカブックっていう１万円くらいの分厚い本がありますが、
あれはヘルプと同じなのでしょうか？
ビニールがかぶっているので立ち読みできず、困っています。
ヘルプと同じなら買う必要もないだろうしなぁ・・・

**名無しさん＠１周年** · NG

日本語バージョンのマセマティカには、ヘルプも
日本語、また日本語訳のマセマティカブックがついてきます。

日本語追加キット（名前忘れ）を英語バージョンマセマティカ
に入れると、ヘルプが日本語になります。また日本語訳の
マセマティカブックがついてきます。

オンラインヘルプにはマセマティカブックが入っています。
しかしヘルプは他のものも入っています。
勉強するのに、オンラインヘルプだと使いにくいです。
その１部（か１章）は必読ということになっています。
損はないので、日本語訳を買いましょう

**便乗して、誰か教えて** · NG

Mathematicaで、不等式の領域を書かせるのって、
どうすれば、よいのでしょうか?
たとえば、x>0, y>0, y<-0.5x+1, y<-2x-2
とかっていうのを書きたいと思っているのですが、
ご存知の方、宜しくお願い致します。

**実習生さん** · NG

ちょっとまっててね。明日調べて報告するから。

**名無しさん＠１周年** · NG

明日？とっくにすぎたぞ。

**実習生さん** · NG

ごめんなさい。しばらくまってください

**名無しさん＠実習生** · NG

>>16
まずは、語学からやるべし！
simulation
シュミレーション　　　×
シミュレーション　　　○

恥っ、ミャー

**実習生さん** · NG

やっぱりできないです。ＭＬできいてくださいな

**実習生さん** · NG

(*次のものをコピーしてペーストすれば良いです。
なにが語学からだ。恥じなのは、これくらいのことを自分で考え
ようとしないヘタレ根性ではないか！！ちなみに、20であげている
例は良くないですね。最後のところはContourPlotでなくてDensityPlotでも行けます
やっとできたフゥ。ポイントはいくつかの条件をIfを使って、
すべての条件を満たすときのみ、１、どれかを満たさないとき０と
いう指示関数を作ったことです*)
If[x <= 0, 0,
If[y <= 0, 0,
If[y >= -0.5x + 1, 0,
If[y >= -2x - 2, 1,
0
]
]
]
];
ContourPlot[f[x, y], {x, -1, 1}, {y, -1, 1}, PlotPoints -> 100]

**実習生さん** · NG

自分を誉めてあげたい　（有森作）

**実習生さん** · NG

うまく答えられるかわからないし、時間がかかるかもしれないけど
、前に考えて勉強になったので、どうぞ質問してくださいな

**実習生さん** · NG

というか、質問があったらメーリングリストに投稿してみては？
mathematicaユーザー会のメーリングリストです。

**実習生さん** · NG

どなたか、マセマティカでの
ニューラルネットワークのプログラムを教えてください。
とくに、競合学習における、サモンマップの導出など
わかりません。よろしくお願いいたします。

**あぼーん** · NG

あぼーん

**実習生さん** · NG

僕の勉強の報告

Mathematica数値数式プログラミング上坂著
を読んで、今は
例題で学ぶMathematica[基礎プログラム編]
のプログラミング例を読んでいる。
Wolframの本は第２章の途中まで読んでいる。
関数よりも、プログラミングのアルゴリズム
を理解することのほうがより大切だと気づきました。

**実習生さん** · NG

今後の予定は、MaederのMathematicaProgrammerとか
Maple/Mathematicaによるシュミレーション工学とか
を読んでみる予定。

ところで、マセマティカのプログラミングに役立つ本を
教えていただけませんか？（いちおうニューメリカルレシピインＣ
は買って持ってますが）

**名無しさん＠１周年** · NG

詳しい事知りませんけど
『Mathematica ニューラルネットワークと最適化』玄光男・井田憲一（共立出版）
というのはどうですか

**名無しさん＠１周年** · NG

>>34
その本、Mathematicaとあんまり関係ないです。

NG

Ｍａｔｈｅｍａｔｉｃａの学生ライセンスは
フルタイムスチューデント（学業専念型で本業を持たない
こと）対象だったと思うので、
社会人学生などは除外されていると思ったが。

**名無しさん＠１周年** · NG

エクセルのデータをMathematicaに取り込んでグラフ化したいのですが、
どのようにやればいいのか教えて下さい。お願いします。

**名無しさん＠１周年** · NG

age

恵子 · NG

複素数を標準系に直した場合ABSを無くす方法ってある？例えば
ComplexExpand[1/(a+b*I)]をすると，
a/Abs[a+bI]^2-Ib/Abs[a+bI]^2になってしまうんだけど，
Abs[a+bI]^2をa^2+b^2にしたい。
フェラでも何でもするから，どなたか教えてください。

**名無しさん＠１周年** · NG

>>39
aまたはbが複素数の可能性があるから、
それ以上は展開できないと思う。

a^2+b^2とできるのはa,bが実数とした場合だから。

恵子 · NG

>>40さんへ
aやｂを実数とするにはどうすればよろしいですか？
ComplexExpandはI以外すべて実数と見なして展開
するとしてありました..
早く成果ださないと教授と寝なきゃいけない、とほほ・・・
お願い助けてください。

恵子 · NG

もう結構ですわ。変換規則をりようしたらできあましたわ。
フェラがほしかったら，もっと早くおしえてくださることねぇ
それではごきげんよう

**名無しさん＠１周年** · NG

ExpandのポリシーとしてAbs[a + I b]^2を展開したりは
しないのかな。
TargetFunctionsを適当に変えてやるとイケる。

恵子 · NG

>>43
まあ、ありがとう、さきっぽだけよ

**名無しさん＠１周年** · NG

寝るのと頑張るのだったら、寝る方がらくじゃん？寝れば？

**名無しさん＠１周年** · NG

上げ。。
　　どうか教えてください。
mathematicaで原始太陽系のシミュレーションしようと思ってるんですが
わからないことがあります。

NDSolveで得られた解を一定の時間間隔ごとにグラフに表示したいんですが
InterpolatingFunctionということなので、できないです。なにか方法はないでしょうか？

**あぼーん** · NG

あぼーん

**Mathematica愛好家** · NG

>>46
Interpolatingした関数でも
Table[{t,f[t]},{t,1,10,1}]
f[t]:Interpolateした関数
{t,1,10,1}:１から１０まで１きざみでTable作る。
とし、ListPlotすればいいです。

**名無しさん＠１周年** · NG

f={{1},{1,2},{1,2,3},{1,2}}
のようなリストからf[[1]]から一つを取り出し、f[[2]]から一つを取り出し、
f[[3]]から一つを取り出し、f[[4]]から一つを取り出すことを考えるとき
のすべての組み合わせを出すにはどうすればよいですか？
教えてください！！

**あぼーん** · NG

あぼーん

**Mathematica愛好家** · NG

>>49
Product[Length[f[[i]]],{i,1,Length[f]}]
とすればいいのでは？
だめ？意味が違う？
質問の答えが１２でよいのなら、僕のやりかたででます。

49 · NG

>>51
ありがとうございます。
すべての組み合わせの個数を出すのではなく、すべての組み合わせを表示させたいのです。
どうすればよいですか？
{{1,1,1,1},{1,1,1,2},{1,1,2,1},{1,1,2,2},{1,2,1,1},{1,2,1,2},
{1,2,2,1},{1,2,2,2},{1,3,1,1},{1,3,1,2},{1,3,2,1},{1,3,2,2}]
のような解がほしいのです。
よろしくおねがいします

**名無しさん＠１周年** · NG

Outer[List, {1}, {1, 2}, {1, 2, 3}, {1, 2}]

49 · NG

>>53
ありがとうございます
それだと、n番目の要素がi個の場合だとできませんよね
{{1}, {1, 2}, {1, 2, 3}, {1, 2}}
の最初と最後の{をとるひつようがありますよね
わからないんですよ。そこが。
教えていただけませんでしょうか。
よろしくおねがいします

**あぼーん** · NG

あぼーん

**Mathematica愛好家** · NG

>>49
In:=
f = {{1}, {1, 2}, {1, 2, 3}, {1, 2}};
myOuter[l_, {s__}] := Outer[l, s];
myOuter[List, f]

Out:=
{{{{{1, 1, 1, 1}, {1, 1, 1, 2}}, {{1, 1, 2, 1}, {1, 1, 2, 2}}, {{1, 1, 3,
1}, {1, 1, 3, 2}}}, {{{1, 2, 1, 1}, {1, 2, 1, 2}}, {{1, 2, 2,
1}, {1, 2, 2, 2}}, {{1, 2, 3, 1}, {1, 2, 3, 2}}}}}
これででます。Mathematicaのユーザー会のメーリングリスト参照

49 · NG

>>56
ありがとうございます。
さんこうになりました。
mathematicaユーザーの会のメーリングリストは調べてみたんですけど、
有料なんですよね？
たまにしか使わないので・・・。

f={{{a, b}}, {{a, b}, {a, c}}, {{a, d}, {a, e}}, {{e, f}}};
のとき、
out:={{{a,b},{a,b},{a,d},{e,f}},{{a,b},{a,b},{a,e},{e,f}},
{{a,b},{a,c},{a,d},{e,f}},{{a,b},{a,c},{a,e},{e,f}}}
のような結果を得るのはどのようにすればいいのでしょうか？
>>56さんの示していただいた例でやってみましたがすべての深さまでの外積
をとってしまいますよね。むずかしいです・・・。
どうすればいいでしょうか？

49 · NG

４９です
よくメーリングリストの入会案内を見たら無料でした。
すいません。
で、早速登録しました。
過去のメールは参照できるのでしょうか？

**名無しさん＠１周年** · NG

質問なんですけど、
siki=i～nまでの四則演算
で式からi～nの文字を読み出して、下の式に代入して
それぞれの文字についてテーブルを作るとかできませんか？
Table[siki,{i,10},{j,100},,,{n,5}]
よろしくお願いします。どなたか教えてください

**Mathematica愛好家** · NG

>>49
mug-ctl@jp-mathusers.org
にタイトルなしで、本文に
get (番号)
をいれると過去の記事がみれます。
本文には英語と数字しか入れてはいけません。
本文に
summary
といれて送ると過去の記事の番号の表が送られてきます。
記事の内容と番号を照らし併せてget文でメールすれば過去の記事が手に入ります。
ちなみに
get 1082-1088
を本文に入れてメールを送れば49さんの質問の内容がのっています。

**Mathematica愛好家** · NG

49さんへ。

In:=
g = {{{a, b}}, {{a, b}, {a, c}}, {{a, d}, {a, e}}, {{e, f}}};
Outer[List, Apply[Sequence, g], 1]

Out:=
{{{{{{a, b}, {a, b}, {a, d}, {e, f}}}, {{{a, b}, {a, b}, {a, e}, {e, f}}}},
{{{{a, b}, {a, c}, {a, d}, {e, f}}}, {{{a, b}, {a, c}, {a, e}, {e, f}}}}}}
になります。Outerの深さを引数として指定すればよいのです。

**Mathematica愛好家** · NG

>>59
siki=i~nまでの四則計算ってのは具体的にどういった計算ですか？
よくわからなくて.......
+,-,*,/がランダムなのですか？

**名無しさん＠１周年** · NG

Apply[List, Unevaluated[{1 + 3 + 4, 2 + 2}], 1]
とすると

{{1, 3, 4}, {2, 2}}と出るので、あとはOuterでリストは
作れるでしょう

**名無しさん＠１周年** · NG

mathematicaにはexcelのテキストボックスのようなものは
ついていないのでしょうか？

49 · NG

>>Mathematica愛好家さん
ありがとうございます。こうすればよいんですね。
今月このプログラムに苦しんでいたんですよ。

59 · NG

siki={a*b*c,d/e+f-g,s*t/u}
でsiki[[1]]からsikiの長さまで
それぞれの変数(siki[[1]]ならa,b,c)についてその変数をさがして
a,b,cそれぞれについて1からnまで代入したすべての組み合わせを出したいのです。
すなわちn=3としてsiki[[1]]は
{1*1*1,1*1*2,1*1*3,1,2*1,1*2*2,1*2*3,1*3*1,,,3*3*2,3*3*3}={1,2,3,2,,18,27}
のような結果を得たいのです
よろしくお願いします

**mathematica愛好家** · NG

>>49
よかったです。

**mathematica愛好家** · NG

>>59さん
仕事の合間に少しずつ考えてるんですけど、スマートな方法が思いつきません。
もうちょっと、考えてみます。

**mathematica愛好家** · NG

>>64さん
Mathematicaで関数すくって、Exelで図にしたいのですか？
それだったら、たとえば任意の関数f[x]に対して、
xが0から10まで0.1きざみで図を書きたい場合、

list=Table[{x,f[x]},{x,0,10,0.1}];
Export["/(好きなDirectory)/aaa.txt", list, "Table"];

でasciiテキストにしてExelで読み込めばいいです。
Directoryの階層についてはsystemによって違うのでhelpでみてください。

**名無しさん＠１周年** · NG

違います。
エクセルでテキストボックスを作ると自由にレイアウトができます
同じようにマセマティカで文章を書いてレイアウトをうまくやるには
どうすればいいかという質問です。

**mathematica愛好家** · NG

>>70
テキストボックスはMathematicaにはありません。
ただ、Mathematicaでは座標の位置やフォント情報を指定すれば、graph上
で軸に名前を付けて表示したりできますが難しくて結構、ストレスがたまります。
私は、Mathematicaで関数など処理したデータをきれいに印刷したりレイアウトを行う際は、
テキストとして出力してExcel,kaleida graph等で改めてプロットし直しています。
テキストでx,yカラムで出力する方法は69でかきましたが、この方法が一番きれいにプロット
できるとおもいます。
お役に立てずにすみません。

**mathematica愛好家** · NG

>>70
あっ、図のレイアウトじゃなくて文章のレイアウトですか？
ワープロソフトではないのでテキストボックスのようにレイアウトするのは
かなり困難でしょう。
様々な人が.nbのNote Bookを公開してますがレイアウトのきれいさはあれくらい
が限界でしょう。
Mathematicaにはいろいろテンプレートとしてスタイルシートがあります。
フォントも変えたりしてレイアウトできますがテキストボックスのような
便利なものはありません。多分。

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

>>72さん。おひさしぶり。
ウルフラムサイトでarXivの形式で文章を作るスタイルシートや
サンプルがダウンロードできます。サンプルをみると数式が
traditionalよりも数式らしく書かれています。4.1じゃないと
完全対応じゃないみたいですが。

4.1になって、4.0とレイアウトがより高度になったのですか？

**mathematica愛好家** · NG

>>74さん。こちらこそ、お久しぶりです。
先ほどwolframのホームページをチェックしてみました。
arXivは論文投稿用のstyle sheetなのかな？
Publiconも投稿用のものですが使ったことがありません。
私はMacOSXを使用してますがOS9でしかうごかないのでまだ試してませんが、デモをみると確かにふつうのスタイルシートより”論文らしさ”がましているような気がします。
ただ、数式のパレットとスタイルシートを組み合わせれば同程度のものはできるとおもいます。
この数式パレットは4.0でも4.1でもあるのでかなりふつうの数式として出力できるので試してみてください。
確かに4.0より4.1のほうがレイアウトがきれいになり、より見やすくなりました。
74さんはMathematicaのノートブックそのままで論文投稿するような仕事をしている人ですか？
私は、学会発表などでもPowerPointに出力（グラフなど）を貼り付けてやっているのであまりスタイルシートも活用していません。
参考にならなくてごめんなさい。

59 · NG

やっぱり、上記の問題が解けません。
引数として[siki,n]としたいのです。
教えていただけないでしょうか？
どんな方法を使ってもかまいません。
よろしくお願いいたします

**名無しさん＠１周年** · NG

もとの問題と、教えて貰ったものをまとめて書いて教えを請わないと
わからないよ。>>76

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

スレ違いかもしれませんが、
ＶＣ＋＋でMathematicaを使うことはできますか？
できればプロットした絵なんかを表示したいんですが。

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

だれかMathematicaを一から勉強するのに最適な本教えて。

**mathematica愛好家** · NG

>>59さんへ
あんまり期待した答えかもしれないけどかいときます。
In[ ]:=
n = 3;
data = Table[i, {i, 1, n}];
list1 = Partition[Flatten[Outer[List, data, data, data]], 3];
list2 = Partition[Flatten[Outer[List, data, data, data, data]], 4];
list3 = Partition[Flatten[Outer[List, data, data, data]], 3];

Clear[a, b, c, d, e, f, g, s, t, u];
Flatten[Table[{siki[[1]] //. {a -> list1[[i, 1]], b -> list1[[i, 2]],
c -> list1[[i, 3]]}}, {i, 1, Length[list]}]]
Flatten[Table[{siki[[2]] //. {d -> list2[[i, 1]], e -> list2[[i, 2]],
f -> list2[[i, 3]], g -> list2[[i, 4]]}}, {i, 1, Length[list]}]]
Flatten[Table[{siki[[3]] //. {s -> list3[[i, 1]], t -> list3[[i, 2]],
u -> list3[[i, 3]]}}, {i, 1, Length[list]}]]

Out[ ]:=
{1., 2., 3., 2., 4., 6., 3., 6., 9., 2., 4., 6., 4., 8., 12., 6., 12., 18., \
3., 6., 9., 6., 12., 18., 9., 18., 27.}

{1., 0., -1., 2., 1., 0., 3., 2., 1., 0.5, -0.5, -1.5, 1.5, 0.5, -0.5, 2.5, \
1.5, 0.5, 0.333333, -0.666667, -1.66667, 1.33333, 0.333333, -0.666667, \
2.33333, 1.33333, 0.333333}

{1., 0.5, 0.333333, 2., 1., 0.666667, 3., 1.5, 1., 2., 1., 0.666667, 4., 2., \
1.33333, 6., 3., 2., 3., 1.5, 1., 6., 3., 2., 9., 4.5, 3.}

です。nを変えれば、3だけではなく4でも5でもいいです。
でも、list1,2,3は変数の数によるので注意が必要です。
ちなみに、//Nはレスを書き込む際にコピーペーストがうまくいかないので
小数ででるようにするために入れたものです。実際に計算するときははずすと分数表示になってみやすいです。
時間がだいぶかかってごめんなさいね。

**mathematica愛好家** · NG

>>79さん
Mathematicaの上でC++を使いたいのですか？
それとも、C++の上でMathematicaを使いたいのですか？
MatheLinkとかlibraryをつかえばC++の上でlibararyとしてMathematicaの組み込み関数などをつかえますよ。

**mathematica愛好家** · NG

>>81さん
Mathematicaを学ぶ方法としてもっとも適していると思うことは、
自分がしたい計算（たとえば微分積分など）についてMathematica Book
（アプリケーションの中にもあるし、製本したものもある、高いけど）
を丁寧によんでいろいろ試してみることだと思います。
いろいろな人がMathematicaの本を書いてるけど、基礎ができてからでないと挫折してしまうおそれがあります。
Mathematica Bookが一番のテキストだと思います。目次をみておもしろそうなのを自分で実行してみて勉強してみましょう！

59 · NG

>>mathematica愛好家さん
ありがとうございます
いまから、試してみます
この計算結果がどうしても学会に使いたいのです

81 · NG

>>mathematica愛好家さん
親切にありがとうございます。
自分が使っているのはstudent版だからMathematica Bookの本は持ってないので
アプリケーションの中にあるのでコツコツやっていこうと思います。

**あぼーん** · NG

あぼーん

**mathematica愛好家** · NG

>>59さん
よかったです。お役にたてて。
学会頑張ってください！

**mathematica愛好家** · NG

>>81さん
Mathematicaいいやつなんで、頑張ってみてください。
ある程度わかってきたら各分野でMathematicaの本が出版されているので
専門書でさらに磨きをかけたり、Mailing Listに登録して（無料）
みんなにいろいろ質問してみてください。
もちろんこのスレでみんなに質問することもいいことだと思いますよ。

59 · NG

>>mathematica愛好家さん
やっと、さきほど自分の研究用のものに教えていただいたプログラムを
改良することができました。
結果は学会の要旨に載せられそうです。

Flatten[Table[{siki[[1]] //. {a -> list1[[i, 1]], b -> list1[[i, 2]],
c -> list1[[i, 3]]}}, {i, 1, Length[list1]}]]

こうすれば、よかったのですね
siki[[i]]の変数がn個の時にも使用できるようにがんばって改良してみます。

**Mathematica愛好家** · NG

>>59さん
Length[list]って以前は書いてしまったけどどこにもlistは作ってませんでしたね。
うっかりレスを書いているときにまちがえてしまいました。
Length[list1]等にしないといけませんね。間違いに気付いてくれて良かったです。
要素がn個の時はなんだかややこしそうですね。
要旨がんばってください！
ちなみに、
In[]:=
Clear[list1, list2, list3, siki];
siki = {e1[[1]]*e1[[2]]*e1[[3]], e2[[1]]/e2[[2]] + e2[[3]] - e2[[4]], e3[[1]]*e3[[2]]/e3[[3]]};
siki[[1]]
siki[[2]]
siki[[3]]
sub1 = Table[e1[[j]] -> list1[[i, j]], {j, 1, 3}];
sub2 = Table[e2[[j]] -> list2[[i, j]], {j, 1, 4}];
sub3 = Table[e3[[j]] -> list3[[i, j]], {j, 1, 3}];
n = 3;
data = Table[i, {i, 1, n}];
list1 = Partition[Flatten[Outer[List, data, data, data]], 3];
list2 = Partition[Flatten[Outer[List, data, data, data, data]], 4];
list3 = Partition[Flatten[Outer[List, data, data, data]], 3];
Flatten[Table[{siki[[1]] //. sub1}, {i, 1, Length[list1]}]]
Flatten[Table[{siki[[2]] //. sub2}, {i, 1, Length[list2]}]]
Flatten[Table[{siki[[3]] //. sub3}, {i, 1, Length[list3]}]]
として置換リストsub1-3までつくり、要素をe1[[j]]などで一般型でかくといいかもしれません。
結果はsikiが前と同じ定義なので同じ結果が出力されるはずです。

**あぼーん** · NG

あぼーん

**あぼーん** · NG

あぼーん

**Mathematica愛好家** · NG

>>59さん
追加です。
e[[j]]として要素を定義しますが、演算子は
siki={Sum[e[[j]],{j,1,n}],Product[e[[j]],{j,1,n}],...}
等にすればいいのでしょうか？
またうまくいけば教えてください。

59 · NG

>>Mathematica愛好家さん
ご教授ありがとうございます

式はほかのソフトウエアから取り込みます。
式の要素はn個で、式の数もi個です。
それを区別するには区切り子として*,/,+,-などをとらえる必要があると思っています
もう少し考えて見ます。

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

MAXIMAやMATHLABは比較的盛況なのに、
Mathematicaはなぜこんなに閑古鳥が泣いているのだろうか？

**名無しさん＠１周年** · NG

JLink使ってる人はいませんか？

**名無しさん＠１周年** · NG

スレの趣旨とは全然関係ないですが、スチューデント版って注文から納品まで
どのくらいかかるものなんでしょう？

最初、注文したときは「一週間くらい」といわれ、２週間たったのに連絡して
こないから聞きにいったら「二、三週間くらい」といわれたのですが、たかが
ソフトの注文にそんなにかかるもんなんでしょうか？
単に生協のおばちゃんが忘れてただけじゃないかと・・・

**名無しさん＠１周年** · NG

>>99
多分、生協のおばちゃんがすぐに卸業者に電話しなかったためでしょう。
ふつう、注文したら一週間ちょっとで、きます。
Mathematica 4.1を注文したときはそれくらいでした。
卸業者の在庫がなくて、という話なら別ですが…

**名無しさん＠１周年** · NG

堺屋太一が恐ろしいシミュレーションをしています。
http://www.geocities.co.jp/WallStreet-Stock/7553/

中国の需要拡大に伴い、世界的な資源危機が発生、
エネルギー資源、希少金属、食糧の国際価格が跳ね上がる。

また、日本は少子高齢化によって若年労働力が激減するため輸出競争力が低下し、
中国の輸出に押される。貿易赤字が急拡大し、2018年には1ドル300円にまで円安が進む。

国内の平均物価は現在の3倍に急騰、消費税は20％に近づく。その上、燃料税も増税されるため
ガソリン市販価格はリッター千円になる。・・・ということらしいのですが。

**λφ³** · NG

>>99
わたしもバージョンアップを頼んでますが、Mathematica4.2日本語版は
WolframResearchの出荷の遅れがあって出荷は１０月になるって言って
ました。

**名無しさん＠１周年** · NG

となると、私の注文も１０月まで待たされることになる可能性が高いって
ことでしょうか。それならそうと生協のおばちゃんもちゃんと報告して欲
しい、、、つーかそもそも約束の期日を越えても連絡してこないってのも
普通の社会なら許されないことでしょうけど。

**λφ³** · NG

>>103
生協の人も怒ってました。WolframResearchに抗議するとか言ってたけど
音沙汰ないなぁ。どうなってんだか。

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

>MAXIMAやMATHLABは比較的盛況なのに、
MAXIMAって別スレッドにありますけどね。
やっぱり世界的に見て人気ある数式処理はMathematica/Maple/Mathcad
数値演算はMatlab
MATHLABって何？

97 · NG

>>106
Matlabの間違いです。すみません。。。。。。

**あぼーん** · NG

あぼーん

**あぼーん** · NG

あぼーん

**あぼーん** · NG

あぼーん

**あぼーん** · NG

あぼーん

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

mgのメーリングリスト（英語版マセマティカメーリングリスト）

でどんなことが話題になっているか、簡単（触りだけ）に紹介しましょう。
登録していない人は
http://smc.vnet.net/mathgroup.html
を開けて、そこのImportant MathGroup Email Addressesの節の中の
Subscribe to MathGroup をクリックすることをお薦めします。

**名無しさん＠１周年** · NG

教えてください．mathematicaで不等式ってできますか？
どうやればいいのですか？初心者なんでお願いします！

**名無しさん＠１周年** · NG

<<Algebra`InequalitySolve`
InequalitySolve[x (x^2 - 2) (x^2 - 3) > 0, x]

とすると答えがでます。その他いろんなことは、ヘルプで
Master Index のボタンをおして、inequalityとしてみてはいかが？

**114** · NG

１１５さんありがとうございました．ヘルプ見ましたが，お教えいただいたように入力したのですが，動きません．
パッケージが足りないとかあるんですかね？

**あぼーん** · NG

あぼーん

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

`InequalitySolve`はごらんのとおり前と後ろにバックスラッシ「`」がついて
います。そのままコピーペーストするとうごきます

**あぼーん** · NG

あぼーん

**114** · NG

119さんありがとうございます．
ところで，ヘルプを参考にしつつ，
実行してみました．
その手順としては，
<<Algebra`InequalitySolve` （enter）
（反応なし）
InequalitySolve[x (x^2 - 2) (x^2 - 3) > 0, x]（enter）
ヘルプではこれで解がでるはずなんですが，次のようにアウトプットされます．
InequalitySolve[x (x^2 - 2) (x^2 - 3) > 0, x]
つまり，解が出てきません．何か足りないのでしょうか？
ちなみにmathematica3.0です．

**114** · NG

119さん，バージョン4.0だったら，不等式解けました．
ありがとうございました．

**名無しさん＠１周年** · NG

age

**名無しさん＠１周年** · NG

mathematica日本語キットをインストールしても日本語になりません。
なんで？

**名無しさん＠１周年** · NG

>>124
いや、日本語キットをいれて立ち上げるとすぐに日本語になりますよ。
ちなみにわたしはMac使ってますけど、OS9では問題なく日本語キット
が動きますよ。OSXでは、OSの環境を日本語にしていれば、
本体をインストールするだけで日本語が表示されます。
ただ、工場から出荷する際、間違えて英語版を日本語版のパッケージに
いれたというミスを聞いたことがありますが、
インストールした際に日本語が表示されないときは英語版なので
交換してくれるそうです。（当たり前ですが）
Windowsにインストールしたのであれば全て再インストールしたほうがいいかも知れませんね。
ただ、再インストールすると、Math IDが変ったりしてPasswordの再発行が必要になるかも知れませんが…

**名無しさん＠１周年** · NG

Windowsで再インストールしても日本語化できない。

**名無しさん＠１周年** · NG

>>126
Mathematicaで使用する言語をJapaneseにしてる？

**名無しさん＠１周年** · NG

オプションインスペクターで日本語にしましょう

**名無しさん＠１周年** · NG

>>128 それってどこの機能？

**名無しさん＠１周年** · NG

age

**名無しさん＠１周年** · NG

ファイル->パレット->言語設定
でも駄目か？

**名無しさん＠１周年** · NG

ファイル->パレット->
に言語設定っていうのがない！
パレットには７つ選択肢があって、
言語設定って英語でなんていいます？
Language Preference とかじゃないですか？
そういうのはないですが。

**名無しさん＠１周年** · NG

age

**名無しさん＠１周年** · NG

>>132
ver4.1 では、GengoSettei になってる。

**名無しさん＠１周年** · NG

mathmaticaでＹの字のように２方向から進行する波が
一つに収束するようなものをアニメーションするにはどうるればいいんでしょうか？

**名無しさん＠１周年** · NG

>>134 GengoSetteiなんてない

**名無しさん＠１周年** · NG

>>134 GengoSetteiなんてない。ver4.0です。

**名無しさん＠１周年** · NG

あげ

**名無しさん＠１周年** · NG

mathematicaで有限要素法電磁界解析って無謀ですか？
Cでやるよりかなり時間かかりそうなんですが。

**名無しさん＠１周年** · NG

>139
なんの目的でそういう事をするのか？

１）結果だけ欲しく、物設計用
仮に Mathematica を持っている、としても
無謀です。実行スピード、メモリサイズで
割に合いません。 Mathematica が動くスパコンを
持っているとしても、構築の手間が膨大です。
教科書コードの丸移植でもコンパイラでしょう。

２）離散化手法自体をあれこれ試す
数学が得意ならば Mathematica でもよいかも。。。
まあそういう人が作ってくれたら神でしょうがね。

**名無しさん＠１周年** · NG

x y A
のようなデータの等高線を描かせて、囲まれた面積（たとえばA=20の面）
の面積を計算させることは可能でしょうか？
等高線を描かせたのはいいが、そこで行き詰っています。。。

**名無しさん＠１周年** · NG

>>141
x,yが均等な間隔なら、A>=20となる点の数を数えれば良いかな。
そうでないなら、ちょっとめんどい。

**名無しさん＠１周年** · NG

>>141
なんかの三次元関数z[x,y]とz=Aの切片がなす閉曲線の面積？
閉曲面の表面の面積？
何をもとめたいの？

**名無しさん＠１周年** · NG

4.2やっと購入しました。OSXの上で動かしています。
何となく全ての動作がはやくなりました。すごいです。
4.0->4.1より4.1->4.2のほうがパワフルです。

**名無しさん＠１周年** · NG

でももう二年たつんですね。このスレ…

**名無しさん＠１周年** · NG

前までmathematicaつかっていたけど、やっぱｃやフォートランのほうが
なにぶんやりやすいぞ。

**名無しさん＠１周年** · NG

>>146
確かに計算が早いけどMathematicaは組み込み関数が便利なので
ずっとMathematica使っています。データのフィッティングも簡単やし。

**名無しさん＠１周年** · NG

Mathematica4.2になってから非線形の最適化プログラムが強化されましたが
そのなかで、NMinimizeというコマンド（非線形simplex法）がありますが実際使ったことのある人、
使い心地を教えてくださいませんか？また、Nelder-Meadもオプションについたみたいなんですけど
これも使ったことのある人、感想おしえてください。

**148** · NG

NMinimizeは変数が二個までしかフィッティングできないので
多量変数（三個以上）の場合はFindMinimumで極小値を求めないとだめっぽい。
NMinimizeはそういうアルゴリズムなんかなぁ？

**名無しさん＠１周年** · NG

NMnimizeのMethodオプションにランダムサーチやNelder-Mead
がついてるし、NMinimizeは非線形simplex法を採用していてかなり
すごい！でも、Nelder-Mead自体は1965年に発表された方法だから
採用するのに時間かかりすぎ？
また、Analysis of Variance (ANOVA)でfittingの精度をきちんと
見積もることができるようになったのはうれしい。

**名無しさん＠１周年** · NG

Mathematica のホームページを見るとNMinimizeでは
「遺伝的プログラミング」を使用できると書いてあるが、
本当か？だったら本当に面白いが、「遺伝的アルゴリズム」
の間違いじゃないだろうか。

・・・と思って英語版レビューのページみたら、これはやっぱり
遺伝的アルゴリズムじゃないか！

**名無しさん＠１周年** · NG

>>151
遺伝的アルゴリズムってどういうイメージ？
よくしらないんで教えてもらえるでしょうか？

**名無しさん＠１周年** · NG

>>152
わかりやすいのがあんまり無いな。
http://florida.mes.titech.ac.jp/ga-j.html

**名無しさん＠１周年** · NG

test

**152** · NG

>>153
ありがとうございました。
さっそく勉強してみます。

**あぼーん** · NG

あぼーん

**名無しさん＠１周年** · NG

須藤元気がもっと強くなりますように。

**名無しさん＠１周年** · NG

mathematicaで多重積分（ｎ重責分）を表現するにはどうすればよろしいですか？
これができないものかと悩んでいます。。。

**名無しさん＠１周年** · NG

Integrate[f[x,y],{x,0,1},{y,0,1}]
じゃだめ？

**158** · NG

2重積分ならわかるんですけどね。
一般にｎ重となるとどうしていいかわからなくて。。。。
∫…∫ｄｘ１ｄｘ２・・・・ｄｘｎ
っていうのをmathmaticaを使って計算させてみたいんです。

**名無しさん＠１周年** · NG

>>160
具体的にどういう計算をしたいか良くわからんが、
Ver4.1以降ならBoole関数使うと上手く行くかも。

**名無しさん＠１周年** · NG

>160
つまり、引数が可変ということ？
だとすると、こんなのはどうでしょう（頭で考えただけで未検証）
１）積分変数とその区間の組を List で作る
２）そのリストに被積分関数を Prepend
３）Apply で List を Integrate

**162** · NG

もしも Integrate[f[r_vector], Domain] なんてのがご希望でしたら
すいません。

大人しくリーマン和かモンテカルロを構成するのかなぁ。。。

**名無しさん＠１周年** · NG

あのお、マスマティカで経路積分がしたいんですが、よｒ

**名無しさん＠１周年** · NG

>164

時間分割タイプの式をちらっとしか見たことありませんが
ちょっとしたソフトウエアプロジェクト並になると思います。
ググったら、２～３の米国大学でシラバスのトピックで
見かけた内容ですね。

領域が変なタイプの場合の対応として
積分領域を含む単純タイプの新領域 +
　　領域内外区別用の Heviside x 被積分関数
にしてみたり

有限分割で式を組み立てて、 Limit もしくは
分割数大で数値的な収束を見る

C とかでライブラリがありそうな気がしますので、
そのソースを解読し、 Mathematica ならではの
リスト処理をするのが早いと思います。

**名無しさん＠１周年** · NG

>>164
無限重積分なんてできっこないだろ！

**名無しさん＠１周年** · NG

単純にxをn回積分するのも簡単にできないのかなぁ？

**名無しさん＠１周年** · NG

n個のデータをn個のファイルに保存するにはどうすればよいですか？

**名無しさん＠１周年** · NG

ひとつひとつ、Exportするのはだめ？
多分、任意のn個のdataを自動的にファイルの名前にaaa_###.txt
としてセーブしたいんですよね？
めちゃめんどうです。

**169** · NG

パターンの認識は他のプログラム言語（perl、awk etc）
に比べてかなり面倒です。
data[[n]]としてわりあてて
Table[Export[data[[n]].txt,data[[n]],"Table"],{n,1,100}]
として発生しておき、実行するしかないでしょう。
やったことないのでちゃんと動くか分かりませんが…

**名無しさん＠３周年** · NG

誤差逆伝搬法のアルゴリズムを教えてください．
各層3層でのさくっと書いてくださる方，お願いします．

BP · NG

>>173
バックプロパゲーションならネットでも簡単にひっかかるだろ
検索してみな

**名無しさん＠３周年** · NG

「ますます出来るMathematica」
http://www.catkk.co.jp/mathcd/index.html

これ使うと勉強効率あがるかしら？

**名無しさん＠３周年** · NG

VBAからmathematica動かしたいんですが、
詳しく書かれてるＨＰとかありませんかね？
wolframのＨＰだけだと、よくわかりません、、、

**名無しさん＠３周年** · NG

age

**200** · NG

200get
ついでにage

_ · NG

http://homepage.mac.com/hiroyuki44/hankaku03.html

**超初心者** · NG

図を描くのに必要なのですが、
ImplicitPlotはどうやったら使えるようになるのでしょうか?
ヘルプにある以下のプログラムを打ち込んでも
下の行の式をそのまま返してくるのですが･･･

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

きちんと各行shift+enterで入力決定しております。
よろしくお願いします。

**名無しさん＠３周年** · NG

>>178
Packageがうまく組み込めていないかもと思われ。
Input menu->Get File Pathで、<<のあとにパスを入れてミソ
きれいな楕円だね？

**超初心者** · NG

>>179
お返事ありがとうございます。
申し訳ありませんが、具体的なプログラムを書いていただけますでしょうか？
<<C:\Program Files\Wolfram Research\Mathematica\5.0\SystemFiles\Graphics\Packages\FileGraphics.m Graphics`ImplicitPlot`
ImplicitPlot[x^2 + 2 y^2 == 3, {x, -2, 2}]
ではダメでした。

**名無しさん＠３周年** · NG

<<"C:\\Program Files\\Wolfram
Research\\Mathematica\\3.0\\AddOns\\StandardPackages\\Graphics\\ImplicitPlot.m"
で、shift+enter
ImplicitPlot[x^2 + 2 y^2 == 3, {x, -2, 2}]
で、また、shift+enter

**超初心者** · NG

>>181
大変ありがとうございます。
実は、何をやったわけではないのですが、何度か新規作成を繰り返していたら
なぜか出来るようになってしまいました。不思議です。

**名無しさん＠３周年** · NG

Solveで連立線形方程式を解こうとしてるのですが、
解きたい変数の数が何十個とあります。

仕方ないので

vars = {v1,v2,…全部の変数を入力してる…,v30}
Solve[ eqns, vars ]

とやってるんですが、これではあんまり大変なので
なにか上手い方法はないでしょうか。

**名無しさん＠３周年** · NG

>>183
var=Table[v[[i]],{i,1,30,1}];
にした少しは楽になるのでは？

>>182
WindowsでMathematica使うとよくそんなことあるの？（パッケージがうまくよめないとか）
Mac使ってますが、一度もそんなことないなぁ

**名無しさん＠３周年** · NG

>>183
Solveじゃなくて行列使って計算するのじゃだめなのかな？

**名無しさん＠３周年** · NG

mathematicaの入門とか、初心者向けの本、教えてください。

**名無しさん＠３周年** · NG

>182,184
横レスですが、Mathematicaは操作の途中で生成した、
オブジェクト変数名を起動時からずっと記憶しますよね。
途中で重複した変数名が残ってしまって悪さをしたのでは？
それら途中履歴をクリアする、Clear文っていうのもあります。

>186
Mahtematica関係の本は入門書からとんでもない高度な話題のまで
いろいろありますが、自分のオススメは、

牧野書店　Mathematica 数値数式プログラミング
ISBN4-7952-0113-7 C3055

**名無しさん＠３周年** · NG

初心者です。shift+enterを押しても、計算結果が出力されません。１＋１を入力してもエラーがでます。
何が原因なんでしょ。

**名無しさん＠３周年** · NG

>>188
1+1を１＋１っていれているからじゃないか？

**名無しさん＠３周年** · NG

>>188
まずは使っているOS, Mathematicaのバージョンなどを書いてください

**名無しさん＠３周年** · NG

マスマティカで　きれいに需要曲線の変化を描きたい
たとえば　ケインジアン交差図と投資関数から導かれるIS曲線
Y=C(Y-T)+I(r)+G　=　a+b(Y-T) + c-dr + G
⇔Y=(1/1-b){a+c+G-bT-dr}
流動性選好理論から導かれるLM曲線
Ld(r.Y)=eY-fr≡Ld(M/P)
⇔r=(e/f) - (1/f)(M/P)
これらを　(r.p.y)空間上で交差させて
この交点から導かれる総需要曲線
Y=(z/1-b)[a+c+G-bT-{d/f-fb)}(M/P)}]
z=f/[f+de/(1-b)] をその接点を切り取る形で表示して
これに　政府がマネーサプライを増加させたときや
租税を増加させたときどのような変化がおこるかわかりやすく図示
またピグー効果が起こった時　C(Y-T,M/P)　=　a+b(Y-T) +ν(M/P)
と消費関数がなるわけですが
このとき　総需要関数がどのように変化するかもわかりやすく
示したいんです
でもmatematica で　数理的に説明するよりも
自分で勝手に三次元グラフを作ってしまったほうが手っ取り早いんですね
Y=(z/1-b)[a+c+G-bT-{d/f-fb)}(M/P)}] この式から見て　Mの増加や
C(Y-{T+kY})　=　a+b(Y-{T+kY}) となる累進課税時に曲線の傾きがどうなるかは簡単にわかりますが
実際にグラフで視覚的に体験して見たいのです
Mの増加は　LM曲線を左方にシフトさせますが
それが需要曲線の変動とどうつながるか
三次元グラフを使って立体的に示したいんです
新着レス 2004/06/25(金) 21:57
88 名前：デフォルトの名無しさん投稿日： 04/06/25 21:57
例えば　コブダグラス型の生産関数です
z = 60 x^0.6 y^0.4
g = Plot3D[z, {x, 0, 10}, {y, 0, 10}, ViewPoint -> {-1, -3, 0.7}]
ContourPlot[z, {x, 0, 10}, {y, 0, 10}, PlotPoints -> 30]

**名無しさん＠３周年** · NG

これを応用して

r = 4 - 0.5y (is曲線)
g = Plot3D[r, {x, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 0.5, 0.7}]
ContourPlot[r, {x, 0, 10}, {y, 0, 10}, PlotPoints -> 30]
164 名前：１３２人目の素数さん投稿日： 04/06/25 11:35
s = 10y - p^2 (lm曲線 1/pではエラーなので-p^2にした)
g = Plot3D[r, {p, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 1, 0.7}]
ContourPlot[r, {p, 0, 10}, {y, 0, 10}, PlotPoints -> 30]
t = Show[r, s]
とすると
Show::"gcomb": "\!\(Show[\(\(2, \
y\)\)]\)のグラフィックオブジェクトの合成処理においてエラーが発生しました．"

と
なってしまう

**名無しさん＠３周年** · NG

解決しました
http://www.geocities.jp/toritsumacro2004/

**名無しさん＠３周年** · NG

パターン・マッチを使って関数を定義する際に，
変数とパラメータの区別ってどうしたらできますか．

例えば，x_^y_ にマッチさせて，
x と y のどれが変数なのかによって，
返す結果を変えたいということなんですけど．

**名無しさん＠３周年** · NG

パターン・マッチを使って関数を定義する際に，
変数とパラメータの区別ってどうしたらできますか．

例えば，x_^y_ にマッチさせて，
x と y のどれが変数なのかによって，
返す結果を変えたいということなんですけど．

**名無しさん＠３周年** · NG

１シートあたりのデータ数の最高数（エクセルだと６５０００くらい）はいくつですか？

**名無しさん＠３周年** · NG

最近来なくなっちゃったみたいだけど、mathematica愛好家って良い香具師だな

**名無しさん＠３周年** · NG

楕円のシュミレーションなんですけど、楕円の回転のさせ方おしえてください。

**名無しさん＠３周年** · NG

バージョン5.1って出てたのか・・・

なんか、面白そうな機能がついてるな。
http://www.wolfram.com/products/mathematica/newin51/equationtrekker.html

**名無しさん＠３周年** · NG

　　　　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　・３・） (　　＾＾）＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ