- 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/
- 153 名前：１３２人目の素数さん mailto:sage [04/06/24 00:18]
- >>152
パレットのSymbolizeをクリック。
ノートブックにフォーカスを写して
expr_<Ctrl+->x<Shift+Enter>
これで定義される。
ただし、<A+B>はキーAを押しながらキーBを押すことをあらわす。
a<Ctrl+->x // FullForm
でa_Subscript_xとなることを確認。
- 154 名前：152 [04/06/24 22:56]
- >>153
ありがとうございます。できました！
ただ、152の計算で定義したあと一度目のShift+Enterでは
入力ミスではないですか。q_Subscript_ｙはp_Subscript_ｙに似ています
・
・
・
と似たようなエラーがいくつか出たんですが、消してもう一度Shift+Enter
を押したらちゃんとできました。
よく分かりませんが、使えるよう勉強させてもらいます。
- 155 名前：１３２人目の素数さん mailto:sage [04/06/24 23:01]
- >>154
> 入力ミスではないですか。q_Subscript_ｙはp_Subscript_ｙに似ています
入力ミスじゃないんだから無視しろ。
- 156 名前：１３２人目の素数さん [04/06/25 00:29]
- どうでもいいよ。
- 157 名前：155 mailto:sage [04/06/25 01:16]
- わかりにくかったか。
>>152
そのメッセージはwarningであってerrorじゃない。親切に入力ミスの可能性を
指摘してくれているだけで、変数の定義自体は正常に行われているので、
入力ミスでない今回の場合は無視すれば良い。
- 158 名前：152 [04/06/25 01:51]
- >>155
アドバイスありがとうございます。
そんな親切な機能があったのですね、分かりました。
- 159 名前：１３２人目の素数さん [04/06/25 02:29]
- mathmatica で 需要関数 Y= いろいろな定数/P の分析で
マネーサプライが増えた時と ぴグー効果が起こった時 の状態を描きたいんですが
まったくわかりません
とにかく
plot してみたんだけど グラフが ショボイっていうか軸により過ぎててわかんない
どうやって調整するんでしょうか
あと それぞれのグラフを色分けしたいんですが タグがまったくわからん
ひとつでも間違うとちんぷんかんぷん
ISLM分析 にPを付け加えた 3次元分析でも それぞれの定数が変化した時に起こる
需要曲線の傾き変動をきれいに江が北印ですが
だれか
ここにタグ貼ってください よろしくお願いします
- 160 名前：１３２人目の素数さん [04/06/25 02:39]
- たとえば plot[Y=1/p,2/p,1/(2p),(p,0,10)]
という基本タグ
をどのように拡張して行けば それぞれの線をきれいに色分けして きれいなグラフで
シフトがわかるようになるのかということです
- 161 名前：１３２人目の素数さん [04/06/25 02:40]
- あと IS-LM モデルのpを加えた三次元表示で P軸できった 総需要曲線が
どのようにシフトするか分析するために
三次現プロットしなければ生けないんですがこのタグも出来たらわかるかたいませんかね
- 162 名前：１３２人目の素数さん [04/06/25 10:50]
- こんなの見つけてきました
例えば コブダグラス型の生産関数です
ほんとは k=1/3 l=2/3 なんですが
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]
- 163 名前：１３２人目の素数さん [04/06/25 11:22]
- Lmもどき
r = 4 - 0.5y
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]
- r = 10y - 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]
isもどき 1/pなんた゜か゜ なんだかエラーになってしまうので変えた
- 165 名前：１３２人目の素数さん [04/06/25 11:40]
- でこれを重ねようとすると エラーになる
- 166 名前：１３２人目の素数さん [04/06/25 12:35]
- ちょっと聞いてくれ
夜露死苦って素数なんだぜ！
分かりにくく答えが出しにくい青春時代を良く表してるよね！
- 167 名前：１３２人目の素数さん [04/06/25 18:44]
- 「ますます出来るMathematica」
www.catkk.co.jp/mathcd/index.html
これ使うと勉強効率あがるかしら？
- 168 名前：１３２人目の素数さん mailto:sage [04/06/25 19:03]
- 付属のヘルプを読む方がはるかにためになるだろうな。
- 169 名前：１３２人目の素数さん [04/06/25 19:31]
- ようやく合成表示がわかったのに エラー なんだけど
r = 2
g = Plot3D[r, {x, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 0.5, 0.7}]
s = y
g = Plot3D[s, {p, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 1, 0.7}]
gjy = Show[r, s]
Show::"gcomb": "\!\(Show[\(\(2, \
y\)\)]\)のグラフィックオブジェクトの合成処理においてエラーが発生しました．"
- 170 名前：１３２人目の素数さん [04/06/25 19:36]
- いやー難しい
ミクロのコブダグラス関数しか表示してるサイト無い
ISLM分析から接線をADであらわしてるサイトなんて無い
だれか うまくシミュレートできるタグわからないかな
- 171 名前：１３２人目の素数さん [04/06/25 19:39]
- タグ覚えるのがこれまた大変なソフトだな
もう使いたく無い
Flash MX よりかは扱いやすいが
ヘルプは醜くて拉致があかない
ミクロ分析の解説本はで照るみたいだけど
マクロも出るの解説本でてないかな
- 172 名前：１３２人目の素数さん mailto:sage [04/06/25 19:40]
- >>169
r = 2, s = y なんだから Show[r, s] = Show[2, y]だね。
Show[2,y]でなにが表示されると？
ここまでの書き込みを読んでいても思ったのだが、君もうちょっと頭使わないと。
- 173 名前：KingOfKingMathematician ◆H06dyzvgzA [04/06/25 20:18]
- お勧めトリップ集
KingOfKingMathematicianの後に付けるのがおしゃれ。
H06dWILLhA : #/{\@%YwX
H06djy9xBA : #SgHdO'H%
H06dYXOYLA : #*｢A@?NVF
H06dhKnt9A : #[Aｼsudｾl
H06dWifa1A : #{SfbN(6ｦ
H06dyzvgzA : #QAiEｼEp- ←使用中
- 174 名前：１３２人目の素数さん [04/06/25 20:48]
- >>172
できませんよ
あなたこそちゃんとタグの使い方マスターしてくださいね
Show::"gcomb": "\!\(Show[\(\(2, \
y\)\)]\)のグラフィックオブジェクトの合成処理においてエラーが発生しました．"
- 175 名前：１３２人目の素数さん [04/06/25 21:11]
- 簡単に行ってしまうと
マスマティカで きれいに需要曲線の変化を描きたいんですよ
たとえば
ケインジアン交差図と投資関数から導かれる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.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 で 数理的に説明するよりも
自分で勝手に三次元グラフを作ってしまったほうが手っ取り早いんですね
- 176 名前：１３２人目の素数さん [04/06/25 21:13]
- 簡単に行ってしまうと
マスマティカで きれいに需要曲線の変化を描きたいんですよ
たとえば
ケインジアン交差図と投資関数から導かれる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 で 数理的に説明するよりも
自分で勝手に三次元グラフを作ってしまったほうが手っ取り早いんですね
- 177 名前：１３２人目の素数さん [04/06/25 21:17]
- Y=(z/1-b)[a+c+G-bT-{d/f-fb)}(M/P)}]
この式から見て Mの増加や
C(Y-{T+kY}) = a+b(Y-{T+kY}) となる累進課税時に
曲線の傾きがどうなるかは簡単にわかりますが
実際にグラフで視覚的に体験して見たいのです
- 178 名前：１３２人目の素数さん [04/06/25 21:19]
- Mの増加は LM曲線を左方にシフトさせますが
それが需要曲線の変動とどうつながるか
三次元グラフを使って立体的に示したいんです
- 179 名前：１３２人目の素数さん [04/06/25 21:20]
- 「マスマティカ」じゃなくて「マセマティカ」だから、あしからず。
- 180 名前：KingOfKingMathematician ◆H06dC8bpwA [04/06/25 22:13]
- Mathematicaは普通の電卓の代わりに使うのは少し難しい。
- 181 名前：KingOfKingMathematician ◆H06dC8bpwA [04/06/25 22:14]
- 一般人が受け入れられるもの、それは関数電卓でもない。
- 182 名前：KingOfKingMathematician ◆H06dC8bpwA [04/06/25 22:15]
- だから高い。
- 183 名前：KingOfKingMathematician ◆H06dC8bpwA [04/06/25 22:15]
- 10年後はどうだろう？
- 184 名前：KingOfKingMathematician ◆H06dyzvgzA [04/06/25 22:16]
- 回し。
H06dWILLhA : #/{\@%YwX
H06djy9xBA : #SgHdO'H%
H06dYXOYLA : #*｢A@?NVF
H06dhKnt9A : #[Aｼsudｾl
H06dWifa1A : #{SfbN(6ｦ
H06dyzvgzA : #QAiEｼEp- ←使用中
- 185 名前：１３２人目の素数さん [04/06/25 22:16]
- 結局わかる人はここにいないんだね
所詮 似ちゃんなんて その程度だね はははははは
- 186 名前：１３２人目の素数さん [04/06/25 22:17]
- KingOfKingMathematician
あんた持って無いの？ ny にいっぱい(ry
- 187 名前：KingOfKingMathematician ◆H06dC8bpwA [04/06/25 22:39]
- Re:>186 吾はmaxima使い。FourierTransformが無いのはちょっと痛いけどね。
- 188 名前：１３２人目の素数さん mailto:sage [04/06/25 22:40]
- Mathematica ４ と Mathematica ５ の違いって何？
- 189 名前：172 mailto:sage [04/06/25 23:42]
- >>174
なるほど>>169の意味も>>172の文意も全く分からないわけね。
頭悪すぎるから諦めた方がいいんじゃないかな。
タグとか言ってるし…
答えてもらえなくて残念だったね。頭が悪くても
人にものを尋ねる態度ぐらいは学習できるかもしれないから、がんばってね。
- 190 名前：１３２人目の素数さん mailto:sage [04/06/26 00:39]
- >>167
1冊何か、なら、マニュアルでもあるWolframの分厚い本ですかね。
あれに書いてないことは「出来ない」ということでしょうから。
初心者用には以下の Mathematica数学の道具箱〈上〉・〈下〉 をすすめます。
www.amazon.co.jp/exec/obidos/ASIN/489242174X/qid=1088177783/sr=1-52/ref=sr_1_2_52/250-9710320-2294634
www.amazon.co.jp/exec/obidos/ASIN/4892421731/qid=1088177783/sr=1-54/ref=sr_1_2_54/250-9710320-2294634
大学１年生くらいでこういう風にグラフを書けることや計算ができることを発見出来ると、人生変わる気がする。
もっとはやく知っとけばよかったと思います。
>>188
このスレに書かれていることに関しては違いが分からない程度。
- 191 名前：KingOfKingMathematician ◆H06dyzvgzA [04/06/26 08:40]
- お勧めトリップ集
KingOfKingMathematicianの後に付けるのがおしゃれ。
H06dWILLhA : #/{\@%YwX
H06djy9xBA : #SgHdO'H%
H06dYXOYLA : #*｢A@?NVF
H06dhKnt9A : #[Aｼsudｾl
H06dWifa1A : #{SfbN(6ｦ
H06dyzvgzA : #QAiEｼEp- ←使用中
- 192 名前：KingOfKingMathematician ◆H06dC8bpwA [04/06/26 10:22]
- 誰か201まで回しといてくれ。
- 193 名前：１３２人目の素数さん [04/06/26 23:00]
- ３DPLOTをするとします。そのためにプログラミングを
ノートブックに書きますよね。それを実行したとき、グラフが
プログラムの直ぐ下に表示されるのですが、グラフだけ別のノートブックに
表示させたいのです。どうすればよろしいでしょうか？
- 194 名前：１３２人目の素数さん [04/06/26 23:49]
- Option inspector を見ましょう
"ctr"+"shift"+"o"
で出てきます
- 195 名前：１３２人目の素数さん [04/06/27 18:04]
- 見ましたが、グラフを別のノートブックに表示させるには
どうすればよいのかわかりません。アドバイスをお願いします。
因みに文字やグラフの大きさは調整できました。助かりました。
- 196 名前：１３２人目の素数さん [04/06/27 19:41]
- <<195
ただ単にコピーandペーストすればよいと思います。
- 197 名前：１３２人目の素数さん [04/06/27 21:03]
- >>189
あのー ここに頼らないで自前で完成したよ
結局ここって使えないねー なにも教えてくん無いし
ほらここ
www.geocities.jp/toritsumacro2004/
- 198 名前：１３２人目の素数さん [04/06/27 21:06]
- mathematicaよりもGNU機能が高いソフトを使ったら直感的に簡単に出来た。
このソフトはC言語書く、つまりプログラミングに慣れていないひとには
ちょっと酷過ぎると思いました。
まっただ視覚的にシミュレートしたかっただけなんだけど
- 199 名前：１３２人目の素数さん [04/06/27 21:14]
- 自慢しすぎ〜アホじゃない？
- 200 名前：１３２人目の素数さん mailto:sage [04/06/27 21:16]
- C言語というかプログラミング一般の文法かな？forとかifとか…
あれは知っていたら便利だし、手計算なんかよりも書く文章量が減って(ﾟдﾟ)ｳﾏｰ なんだから、
そこまで酷評するとこでもないよ。
- 201 名前：１３２人目の素数さん [04/06/27 21:18]
- >>198
GNU機能これ何、教えて
- 202 名前：１３２人目の素数さん [04/06/27 21:23]
- GUI が充実しているの間違いだった
- 203 名前：１３２人目の素数さん [04/06/27 21:24]
- Graphical User Interface の略です
- 204 名前：１３２人目の素数さん [04/06/27 21:31]
- Mathematica って使い安いですか？
なんかコマンド覚えるのが苦痛で。
次期バージョンではＣＵＩからＧＵＩへ、
使い勝手の劇的なる向上を 期待してます。
- 205 名前：１３２人目の素数さん [04/06/27 21:33]
- mathematicaでplotした3Dをドロップして拡大縮小したり マウスだけでいろいろな面で
切断できるようにしてほしい
そしてその切断面でアニメを再生できるように
- 206 名前：１３２人目の素数さん [04/06/27 21:35]
- 例えば サインカーブの変動を３次元でいろいろな面に切って観察したいと思いませんか？
この操作がマウスだけでできるようにして欲しいのです
あとマウスで作った立体図形などに 最も近時出来る関数を一発表示してくれるとか
- 207 名前：１３２人目の素数さん [04/06/27 21:38]
- あと面できらずにコサインカーブ面で切ってその表面でのサインカーブの
表面を観察したりして見たい。
マスウだけで
- 208 名前：１３２人目の素数さん [04/06/27 21:42]
- そしてその表面上にある点を設置して
サインカーブの変動と共に移動する点の座標をリアルタイムで観察したい
この面はコサインで切ってあるので消えたり現れたりする
- 209 名前：１３２人目の素数さん [04/06/27 21:44]
- これができると例えば 三次方程式の解の座標を
その係数項の変化と共に観察できるようになる
- 210 名前：１３２人目の素数さん [04/06/27 22:36]
- >>196 「ただ単にコピーandペーストすればよいと思います。」
それでは意味がありません。計算結果をグラフで出力する場所を
別のノートブックに設定したいんです。いちいちコピペしていたのでは
面倒ですし、コピペだとプログラミングの直後にグラフが表示されてしまうことが
解決できません。
どなたかまっとうな答えをご存知ないでしょうか？
評価するたびに別枠でノートブックが自動的に開いて
そこにグラフのみ表示させるようにしたいのです。
MATLABだと簡単にできるのですが、マセマティカではどうするのでしょうか？
- 211 名前：１３２人目の素数さん [04/06/27 22:59]
- mapleでは うまくできますよ
{>がつくところがインプットするところなので
--
|>
|>
|>
|>
|>
--
としてインプットする場所を増やせば下に一括表示されます
ただ1800ドルするのと、販売地域がカナダとアメリカだけなところが痛いです
- 212 名前：１３２人目の素数さん mailto:sage [04/06/27 23:00]
- マセマティカでも
--
|
|
|
|
|
--
とすればいいだけだったんだ 解答になって無いか
- 213 名前：１３２人目の素数さん mailto:sage [04/06/27 23:15]
- https://check.cybernet.co.jp/maple/request/request.shtml
日本語出てるのか
- 214 名前：１３２人目の素数さん mailto:sage [04/06/27 23:22]
- 金持ちの使えるソフトだからマセマティカとは格が違うな
大卒初任給より高い
- 215 名前：１３２人目の素数さん mailto:sage [04/06/27 23:30]
- 今9だけどあそこには英語の7がかろうじてあった
- 216 名前：１３２人目の素数さん mailto:sage [04/06/27 23:32]
- 高校生のとき買っとけば良かった 大学以上だと17万かよ
このライセンスは、学校教育法によって定められた、高等学校・高等専門学校向けに提供されます。
■ シングルライセンス
ライセンス料金
1ユーザ 47,250円（本体価格45,000円）
- 217 名前：１３２人目の素数さん [04/06/27 23:42]
- >>208
できるね
www.cybernet.co.jp/maple/hiroba/nami/nami.html
- 218 名前：１３２人目の素数さん [04/06/28 00:13]
- www.cybernet.co.jp/maple/hiroba/digital_filter/digital_filter.html
なんかよくわかんないけど ものすごいソフトだな これって
- 219 名前：１３２人目の素数さん mailto:sage [04/06/28 01:09]
- >>217-218
この程度ならMathematicaでも簡単に作れるが・・・
インターラクティブにしたいというなら、面倒だな。LiveGraphics3Dなんかと
組み合わせないと。
- 220 名前：１３２人目の素数さん mailto:sage [04/06/28 16:44]
- >>210
こんなものかな？
動きが気に入らなかったらHelp引いて直してくれ。
newNbExec[f_, args___] := Module[{new = NotebookCreate[]},
NotebookWrite[new,
ToString[HoldForm[f]] <> "[" <> ToString[HoldForm[args]] <> "]", All];
SelectionEvaluate[new];]
newNbExec[Plot, Sin[x], {x, 0, 1}]
- 221 名前：１３２人目の素数さん mailto:sage [04/06/28 16:49]
- SetAttributes[newNbExec, HoldAll];
newNbExec[f_] := Module[{new = NotebookCreate[]},
NotebookWrite[new, ToString[HoldForm[f]], All];
SelectionEvaluate[new];]
newNbExec[Plot[Sin[x], {x, 0, 1}]]
のほうが綺麗か。
- 222 名前：１３２人目の素数さん [04/06/28 21:08]
- mathematica4.0とmaple8(5,6)持ってるけど、個人的にはmapleの方が気に入ってる。
mathematicaはもうアップグレードできなくなったし。
- 223 名前：152 [04/06/28 21:28]
- ax,ay,az,bx,by,bz（ｘｙｚは下付けの添字）というシンボルを変数にし
Solve関数を使って多元連立方程式を解いたのですが、解がax,ay,az,bx,by,bzの
順番で出てこないのを、順番通りにしたい場合どうすればいいのでしょうか？
Sort関数などを使ってもダメですよね。
- 224 名前：KingOfKingMathematician ◆H06dC8bpwA [04/06/28 21:44]
- Re:>223 List関数とPart関数を組み合わせて無理やり作るのはだめ？
- 225 名前：UltraMagic ◆NzF73DOPHc [04/06/28 21:45]
- NO!
また名前間違えた。
- 226 名前：１３２人目の素数さん mailto:sage [04/06/28 22:14]
- >>223
> Sort関数などを使ってもダメですよね。
ダメではない。試してないか、Sortの使いかたがわかってないか。
ans = Solve[...]
の答えは {{...}}だから、内側の{...}の中身をソートするためには、
Sort /@ ans
としなくてはならない。
- 227 名前：223です [04/06/29 17:34]
- >>226
ありがとうございます。
どうもまだ、Sort関数の使い方が分かってないみたいです。
Sort[Solve[連立方程式,{変数}]]
や
Solve[連立方程式,{変数}];
Sort[%]
などはやってみたんですがこれではダメなわけですよね。
さらには、
Solve[連立方程式,Solt[{変数}]]
なんていう変なのまで試みたんですが当然ダメでした。
どうもSolveのansの内側の{}をソートするってのがうまくできないんですが
よろしければ、教えていただけないでしょうか。
- 228 名前：１３２人目の素数さん mailto:sage [04/06/29 18:12]
- >>227
？？？>>226に書いてあるように、
ans = Solve[<連立方程式,{変数}>];
Sort /@ ans
と入力すればいいはずだが？
それでできないんなら実際にどういう式をいれたか見せてくれ。
- 229 名前：223です [04/06/29 18:42]
- >>228
何度も申し訳ありませんでした。できました、ありがとうございます！
- 230 名前：１３２人目の素数さん [04/06/29 20:04]
- どなたか教えてください。
p1={a, b, c, d}
p2={e, f, g}
p3={h, i, j, k}
p4={l, m}
の要素数の違うリストp1からp4まで要素をひとつずつ取って、
{a, e, h, l}, {a, e, h, m}, {a, e, i, l} …
のような感じで重複なく全種類の要素数４のリスト作りたいのです。
おねがいします。
- 231 名前：UltraMagic ◆NzF73DOPHc [04/06/29 21:55]
- Re:>230 List関数をPart関数とDo関数を使え。
- 232 名前：１３２人目の素数さん mailto:sage [04/06/30 03:42]
- >>231 泥くさそうだな。
>>230
BindListM[lst_List, f_] := Join @@ Map[f, lst];
ReturnListM[x_] := {x};
MyOuter[] := {{}}
MyOuter[l_List, ls___List] :=
BindListM[l,
Function[x,
BindListM[MyOuter[ls], Function[tl, ReturnListM[Prepend[tl, x]]]]]]
MyOuter[p1,p2,p3]の様に使う。
# このコードが良く分からなければHaskellという言語を
# 通じてMonadを勉強するといい。シンタックスシュガーによって
# Haskellでは上のコードがこう書ける。
# myouter [] = [[]]
# myouter (l:ls) = [(x:tl) | x <- l, tl <- myouter ls]
しかし実は便利な組込み関数Outerがあるので、
Flatten[Outer[List,p1,p2,p3,p4],3]
(後ろの3はpjの個数 - 1)
でよい。
- 233 名前：１３２人目の素数さん mailto:sage [04/06/30 03:51]
- 実用するなら、MyOuter[ls]が何回も計算されないように:
MyOuter[] := {{}}
MyOuter[l_List, ls___List] :=
Module[{xs = MyOuter[ls]},
BindListM[l,
Function[x, BindListM[xs, Function[tl, ReturnListM[Prepend[tl, x]]]]]]]
- 234 名前：230 mailto:sage [04/06/30 15:06]
- ありがとうございます。
Outerを使いたかったので、
Flatten[Outer[List,p1,p2,p3,p4],3]
が大変わかりやすかったです！
- 235 名前：１３２人目の素数さん [04/06/30 18:38]
- なんか昔からMathematica4.0の起動が不安定で
普通に使えるときとか、ダイアログも何も無しに終わっちゃうときとか
パラパラあったんですが、最近どうも起動失敗ばっかりで、
Wolframのサイトのヘルプも試してみたんですがどうもダメで、
他にこんな症状の人います？
- 236 名前：１３２人目の素数さん mailto:sage [04/06/30 18:46]
- 誰か、自分の気に入ってる面白いコマンド教えてくれ
- 237 名前：１３２人目の素数さん [04/06/30 18:50]
- Limit[
Limit[
Integrate[a x/((x - b)^2 + a^2), {x, -L, L}],
{L -> ∞}],
{a -> +0}]
これを現在実行してるんですけど、１５分くらい経つけど答えでないんですよね〜。
バグですかね〜。
普通に問題とけるぞ！って方は教えて頂けるとありがたいです。
- 238 名前：１３２人目の素数さん mailto:sage [04/06/30 19:00]
- >>237
置換で普通にできないべか
- 239 名前：１３２人目の素数さん [04/06/30 19:02]
- >>238
具体的に教えていただけないでしょうか？
- 240 名前：１３２人目の素数さん mailto:sage [04/06/30 19:03]
- >>239
おいおい、わからない質問スレへ行けよ
- 241 名前：１３２人目の素数さん [04/06/30 19:04]
- >>240
わかりました。
- 242 名前：１３２人目の素数さん mailto:sage [04/06/30 19:31]
- >>237
バージョンはいくつ？
4.2.1だと一瞬でb \[Pi]とでるぞ。5以上ならサポートに報告すれ。
Limit[Limit[Integrate[a x/((x - b)^2 + a^2), {x, -L, L}], L -> Infinity],
a -> 0]
- 243 名前：１３２人目の素数さん mailto:sage [04/06/30 19:43]
- ４．０でもｂπ出たぞ。
- 244 名前：１３２人目の素数さん mailto:sage [04/06/30 20:48]
- >>242
5.01だと無限ループに入るみたいだね。
これに限らず、ver5以上は式の簡略化などでフリーズする
確率が高いような気がする。
- 245 名前：１３２人目の素数さん mailto:sage [04/06/30 22:51]
- それはフリーズではなく計算中
- 246 名前：１３２人目の素数さん mailto:sage [04/06/30 23:06]
- >>245
あー、すまん。フリーズじゃなくてカーネルエラーで落ちるの間違い。
前スレにも出したけど、こんな式を評価すると確実に落ちる。
7910*Csc[(4*Pi)/15]*Sqrt[((-(1 -
Cos[82000/797267]^2))*Sin[(4*Pi)/15]^2 +Sqrt[(1 -
Cos[82000/797267]^2)^2*Sin[(4*Pi)/15]^4 - (-1 +
Cos[82000/797267]^2)*Sin[(4*Pi)/15]^2*(Cos[82000/797267]^2
- Sin[(4*Pi)/15]^2)])/(Cos[82000/797267]^2 -
Sin[(4*Pi)/15]^2)]//FullSimplify
- 247 名前：１３２人目の素数さん [04/07/01 22:27]
- いやー、mathematicaは楽しい言語だねー。はまっちゃったよ。
- 248 名前：１３２人目の素数さん mailto:sage [04/07/02 08:27]
- 言語か？
- 249 名前：１３２人目の素数さん [04/07/02 09:20]
- >>248 おいおい、247の言うようにmathematicaはプログラミング言語だよ。
単なる優れた電卓と思ったら大間違いだよ。
「Mathematica 方法と応用」 がお薦め。
- 250 名前：１３２人目の素数さん [04/07/03 20:19]
- ちょいヒマな人教えてくれないかな？Mathmaticaで、次の極限値出せる？
lim[n→∞]∫[0→π/2]｜sin2nx*cosx｜dx
コマンド教えてくれないか？
- 251 名前：１３２人目の素数さん [04/07/03 20:48]
- >>250
Limit[Integrate[Sin[2 n x] Cos[x], {x, 0, Pi/2}], {n -> ∞}]
答えは0でした。
- 252 名前：１３２人目の素数さん mailto:sage [04/07/03 20:48]
- >>251
絶対値は？
- 253 名前：１３２人目の素数さん [04/07/03 20:54]
- >>251
手計算でやったら２/πになったんだけど・・・間違いかな
手計算でやってもかなり莫大な計算量になったから、
Mathematicaのコマンド教えてもらおうと思ったんだけど・・・
|
 |
