〓 Mathematica 四 〓
レス数が950を超えています。1000を超えると書き込みができなくなります。
0001132人目の素数さん
2007/12/10(月) 16:00:01, - ' ,..、、.ヾ{{フ'⌒`ヽ、
/ ,:', -‐‐` ´ '´⌒ヽ ヾ:、
. ,' ,'´ ,ィ ,ィ ,' , `ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,))
! | ! .,'-.{ ! !|; |`、.}゙!.! |. ! ヽ.
', ', |Vァ=、゙、 `゙、!-_:ト,リ', l ! | ゙',
ヽ、', l:!Kノ}. f:_.)i゙i: リ ! l ル
| l!iヾ- ' , .!__:ノ ゙ ,リ l リ'´
. ',|!!、 r‐┐ ` ノ'. /,イ / ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄
'i!゙、ヽ、 ゙ー' _, ィ,:',:''´ < Mathematica に関する話題はここに書いてね!
゙:、ィ、jヾー::: 'iヘ .ノ',リ. \___________________
,、- '´ ヽ、゙、 { `>"、
/\\ ', } //`ヽ
過去スレ
>>2
0002132人目の素数さん
2007/12/10(月) 16:03:02http://science2.2ch.net/test/read.cgi/math/985023298/
〓 Mathematica 2 〓
http://science3.2ch.net/test/read.cgi/math/1078534285/
〓 Mathematica 3 〓
http://science6.2ch.net/test/read.cgi/math/1121413040/
0003にょにょ ◆yxpks8XH5Y
2007/12/10(月) 16:04:410004さくら姫
2007/12/10(月) 16:18:56,イ /_ -‐――‐‐- _ ヽ ヽ
. / ,イ/ / ` ヘ ヘ
/ ./ / / / / ∧ ∧ ハ
. / ./ / / / / ハ ハ l
i ./ / / / /,' { l ∧ ハ |
|/ / ,'. ! / /. ! ハ. lト、 ,ハ | ハ l |
| i i | /| .ハ | | { | i |l ト、 }l || |
| | | | l {Hl {ナT爪 | V | | 爿 辷} 八 リ l l
| | | ,ハ. ト |ィfテ示kヽ| ヽl }/ 」/rzッ/_|/ ,' l
ヽ ト、ト、 ∧ l! /|l{、__,ハ| ´/{ }心/}. / ,'
. 〉、 ∧l l刈刋リ {州爪リ }/,'/ /
l \ ヽ. ヽ--'' ヽrzン / /
|| 八 { Xxx ! xxX / /
ヽト、 \ ー ‐ /./ /
ヽ\ ヽ、 /,.' イ|
| ヽハ、> 、_ ,. ィ/ //|i!
| `¨ヽk}> 、 __,. r<{r-く- ''" リ >>1 スレ立て乙です
l / ̄\ / ̄ヽ ,'
l / `ー‐''´ ∧ /
ヽ _,r〈 /⌒V ハ-、_/
,.r<\ | / l ̄l V / /\
,. -<. | / ̄. . ̄ト、 / `丶、
/ lト-< └┐ .┌┘ ヽンl 丶、
/{ l |l └‐┘ || l ヽ
. / .| l || || l ト、
0005132人目の素数さん
2007/12/10(月) 16:21:540006132人目の素数さん
2007/12/10(月) 18:54:43y''=3とかの微分方程式って
DSolve[y''[x]==3,y[x],x]
て、入力してもとけないんですけど、どうやったらとけますか?
すいませんが、教えてください。
0007132人目の素数さん
2007/12/10(月) 19:13:53解けるけど、、、。
{{y[x] -> (3 x^2)/2 + C[1] + x C[2]}}
0008132人目の素数さん
2007/12/10(月) 20:17:520009132人目の素数さん
2007/12/10(月) 23:58:14「
:関数取得( := )(Input入力)
:変数取得( = )(Input入力)
:関数定義( := )(Module定義)
:関数定義( := )(Module定義)
:定義関数呼び出してのメイン演算、条件文込み
」
てのが一行ずつのインタプリンタでは上手く実行されてるとして、これを
通常のプログラムみたいに何かしらの一個の読み込みコマンドを入れたら
全部包括的にやってくれるようにまとめ直したいんだけど、
どうすればいいんですか?
あと、MessageDialogue、Input、ChoidceDialogの実行に明確な順列をつけるには
どうすればいいんですか?
抽象的な質問で非常に申し訳ないんだけど、参照するコマンドとかだけでもわかればありがたい。
0010132人目の素数さん
2007/12/13(木) 04:20:0012394, 12367, 12384, 12425, 12394, 12356, 12418, 12398, 23455,
34892, 12375, 12390, 12288, 12393, 12358, 12379, 20170, 24180,
12418, 19968, 20154, 12398, 12463, 12522, 12473, 12510, 12473,
20057, 65367, 65367, 65367};
FromCharacterCode[str]
0011132人目の素数さん
2007/12/13(木) 14:58:00なくだらないもの実
行して どうせ今年
も一人のクリスマス
乙www
0012132人目の素数さん
2007/12/14(金) 03:10:10109, 105, 116, 115, 117, 46, 99, 111, 109, 47, 103, 97, 109, 101,
47, 99, 111, 109, 105, 110, 103, 47, 50, 48, 48, 55, 47, 48, 55,
47, 48, 51, 47, 104, 45, 49, 48, 52, 95, 55, 52, 51, 56, 52, 95,
114, 52, 53, 46, 106, 112, 103, 46, 106, 112, 103};
Import[FromCharacterCode[str4]]
0013132人目の素数さん
2007/12/15(土) 08:55:59{x -> u ,y -> v}
と求めた後、u,vをそれぞれx,yの値とするにはどうしたらいい?
いまはいちいち手入力で
x=u;
y=v;
とか打ってるんだけど、面倒で。
あほーな質問でスマソ。
0014132人目の素数さん
2007/12/15(土) 11:38:22{x0, y0} = {x, y} /. {x -> u, y -> v}
という感じで解をx0,y0に代入してる
x,yに解を代入しないのは、もう一度x,yで方程式記述するとき
すでに値が入っていると思うように動かないので、それを避けるため
0015132人目の素数さん
2007/12/15(土) 16:19:35{x, y} = {x, y} /. %[[1]]
さらに>>14方式なら
{x0, y0} = {x, y} /. %[[1]]
でおk
0017132人目の素数さん
2007/12/16(日) 02:14:45ttp://www.wolfram.com/products/player/download.cgi
ttp://www.wolfram.co.jp/products/player/download.cgi
を開いても真っ白のままで何も起こらないです
Windows XP+IE6とVista+IE7で両方同じ現象。
0018132人目の素数さん
2007/12/16(日) 23:14:12XP+IE7で確認しました。文字コードの影響みたい
右クリック -> エンコード -> その他 -> Unicode(UTF-8)
0019132人目の素数さん
2007/12/17(月) 09:52:38せっかくなので登録画面のご意見欄にこの問題のこと書いておきます
0020132人目の素数さん
2007/12/19(水) 21:31:54業務の関係でmathematicaを購入することになりました
どこでどのようにして購入すればよいのでしょうか?
やはりウルフラムから購入するのでしょうか・・・
0021132人目の素数さん
2007/12/19(水) 21:32:560022132人目の素数さん
2007/12/19(水) 23:05:210023132人目の素数さん
2007/12/20(木) 12:07:24はい。
大きなソフト売り場を持つ店舗なら、店頭でも買えますよ。
0024132人目の素数さん
2007/12/20(木) 23:52:32積分や面倒な計算も結構すぐやってくれる。記号計算が出来る使いやすいソフトは
やはりMathematicaだと思う。昔はMapleの方が性能が良いと言われていたけど今は
どうなのかな?Mapleは使い慣れないせいか、ほとんど使ってない。
0025132人目の素数さん
2007/12/22(土) 01:03:45新スレおめでとうございます
, _ ノ)
γ∞γ~ \
| / 从从) )
ヽ | | l l |〃/VVVVVVVVVVVVVVVVVVVVVVVVN\
`从ハ~ ワノ∩ ハニャウェ━━━━━━━━━━━━ブ >εε=ヽ( `Д´)ノ ウワァァァァン
⊃ .VVVVVVVVVVVVVVVVVVVVVVVVVN/
, _ ノ)
γ∞γ~ \
| / 从从) ) 人人人人人
ヽ | | l l |〃 ニコニコ ≪ヽ( `Д´)ノ ≫ ビリビリ
`从ハ~ ワノ) ^Y^Y^Y^Y^Y
, _ ノ)
γ∞γ~ \
| / 从从) )
ヽ | | l l |〃 ハニャーン! ヽ(・∀・ )ノ ハニャーン!
`从ハ~ ワノ∩
0026132人目の素数さん
2007/12/23(日) 16:22:38わざわざウルフラムから買わなくても
日本に代理店とかあったんじゃなかったっけ?
0027132人目の素数さん
2007/12/25(火) 00:16:520028132人目の素数さん
2008/01/07(月) 12:35:06Mathematicaを再インストールするときは、パスワードを再発行してもらう必要はありますか?
MathIDって変わってしまうんでしょうか?
0029132人目の素数さん
2008/01/08(火) 13:30:26MathIDってことはスチューデント版だと思うんですが、
再インストールすると確かに変わりますので再発行が
必要です。
0030132人目の素数さん
2008/01/09(水) 13:37:35ありがとうございます。
再発行してもらいます。
0031132人目の素数さん
2008/01/11(金) 22:51:24リモートデスクトップからわざわざ立ち上げるという手間がわずらわしい。。。
誰か教えてください。
mathmatica5でコンピュータはXPです。
0032132人目の素数さん
2008/01/31(木) 02:37:090033132人目の素数さん
2008/01/31(木) 02:38:070034132人目の素数さん
2008/01/31(木) 02:47:03カーネルは超能力者です
0035132人目の素数さん
2008/02/01(金) 00:35:17Mathematica6.0Student@Vistaで使ってます。
Export関数でFLV形式がサポートされた、ということで出力しようと思ったのですが、出力されるFLVファイルが再生できません。
だいたいファイルサイズからして7.4KBとビデオの大きさではありません。
入力した式は
Export["MovingPower.flv", Manipulate[Plot[x^a, {x, 0, 5}, PlotRange -> {{-1, 5}, {-1, 5}}], {a, -2, 2}]]
です。
再生環境はffdshow+MediaPlayerClassicです。
別のFLVは再生できるのでコーデックの問題ではないと思います。
「式がちげぇよトンチキ」「私もできない」などありましたら教えてください。よろしくお願いします。
あと何でさくらタンなんですか。CCS板住人としてはここで彼女に会えるとは思っていなかったハァハァ。
0036132人目の素数さん
2008/02/01(金) 05:56:570037132人目の素数さん
2008/02/01(金) 08:52:17できればkwskお願いします。
0038132人目の素数さん
2008/02/09(土) 21:57:400039援護保守
2008/02/10(日) 05:44:27119, 46, 108, 117, 99, 107, 121, 45, 99, 104, 46, 99, 111, 109, 47,
105, 110, 102, 111, 47, 105, 109, 103, 47, 105, 110, 102, 111, 95,
101, 118, 48, 55, 49, 50, 48, 50, 95, 114, 101, 112, 111, 114,
116, 95, 116, 105, 116, 46, 106, 112, 103}]]
0040132人目の素数さん
2008/02/11(月) 20:05:320041132人目の素数さん
2008/02/24(日) 20:34:490042132人目の素数さん
2008/02/24(日) 22:38:580043132人目の素数さん
2008/02/24(日) 22:52:05きみは外資金融にでも内定が決まっているのかい?
0044132人目の素数さん
2008/02/25(月) 00:16:30英語版ならバイナリと鍵genが至る所に出回っているよ。
まぁ、しかし実際便利な道具だと思う。面倒な計算は大抵やってくれるし、いちいち本をひっくり返す必要がない。
0045132人目の素数さん
2008/02/25(月) 21:09:26関数型プログラミングできるし簡易スクリプト言語として便利そう
0046132人目の素数さん
2008/02/25(月) 22:36:060047132人目の素数さん
2008/02/26(火) 00:16:280048132人目の素数さん
2008/02/26(火) 12:42:10つStudent
0049132人目の素数さん
2008/02/26(火) 13:55:370050132人目の素数さん
2008/02/26(火) 23:16:1300511stVirtue ◆.NHnubyYck
2008/02/26(火) 23:20:35ところで、私に通貨をください。
0052132人目の素数さん
2008/02/27(水) 16:30:08マチマチカにもクローンの類はあるの?
0053132人目の素数さん
2008/02/27(水) 16:44:5000541stVirtue ◆.NHnubyYck
2008/02/27(水) 17:33:36しかし、 アメリカ合衆国の公用語に通じていないと厳しいかもしれぬ。
0055132人目の素数さん
2008/02/28(木) 13:35:35マニュアル?
豊富な例を自分の手でいじってると
大抵のことは分かる気もする
0056132人目の素数さん
2008/02/29(金) 10:01:250057132人目の素数さん
2008/03/01(土) 11:35:08こんにちは、
下記の式、y1とy2があります。
y1を評価させて、y2の通り表示させるには、どうしたら良いでしょうか?
y1 = Expand[(a1*p1 + a2*p2 + a3*p3 + b*mc)^2]
y2 = p1^2*a1^2 + p2^2*a2^2 + p3^2*a3^2 + mc^2*b^2 + (a1*a2 + a2*a1)*p1*p2 + (a2*a3 + a3*a2)*p2*p3 + (a3*a1 + a1*a3)*p3*p1 + (a1*b + b*a1)*p1*mc + (a2*b + b*a2)*p2*mc + (a3*b + b*a3)*p3*mc
ちなみに、y1をそのまま評価しますと、下記になり、y2の通り表示しません。
p1^2*a1^2 + 2*p1*p2*a1*a2 + p2^2*a2^2 + 2*p1*p3*a1*a3 + 2*p2*p3*a2*a3 + p3^2*a3^2 + 2*mc*p1*a1*b + 2*mc*p2*a2*b + 2*mc*p3*a3*b + mc^2*b^2
0058132人目の素数さん
2008/03/02(日) 17:47:41これは、はまるべきでないパターン。
単に論文や、教科書が正しいかどうか見たいだけなら
In[18]:= y1=Expand[(a1*p1+a2*p2+a3*p3+b*mc)^2]
Out[18]= b^2 mc^2+2 a1 b mc p1+a1^2 p1^2+2 a2 b mc p2+2 a1 a2 p1 p2+a2^2 p2^2+2 a3 b mc p3+2 a1 a3 p1 p3+2 a2 a3 p2 p3+a3^2 p3^2
In[17]:= y2=p1^2*a1^2+p2^2*a2^2+p3^2*a3^2+mc^2*b^2+(a1*a2+a2*a1)*p1*p2+(a2*a3+a3*a2)*p2*p3+(a3*a1+a1*a3)*p3*p1+(a1*b+b*a1)*p1*mc+(a2*b+b*a2)*p2*mc+(a3*b+b*a3)*p3*mc
Out[17]= b^2 mc^2+2 a1 b mc p1+a1^2 p1^2+2 a2 b mc p2+2 a1 a2 p1 p2+a2^2 p2^2+2 a3 b mc p3+2 a1 a3 p1 p3+2 a2 a3 p2 p3+a3^2 p3^2
In[19]:= y1==y2
Out[19]= True
以上でy1=y2より十分。
どうしてもこの形に意味があるならだけど、このy2の形はあくまで
「筆者の人間側の都合でこうなった」というのを忘れてはいかん。
Collectとかで係数をまとめればそれらしき形にはなってくるが 規則がわからんので
完全な一致をMathematicaの関数のみで達成するのは難しいと思うぞ。
0059132人目の素数さん
2008/03/07(金) 06:04:35またアップデート自体はどのように行いましたか?
0060132人目の素数さん
2008/03/07(金) 17:19:51こんな形してるんだって思った。
でもプレステのガンダムseed destinyのほうがすごくない?
動きは早いしまるで3次元空間にいるかのようだ。
plot3Dのf(x,y)曲面よりガンダムのほうが複雑な形してるじゃん。
ガンダムとマテマテカはどっちのほうがコードの量が多いんでしょうか?
0061132人目の素数さん
2008/03/07(金) 18:46:380062132人目の素数さん
2008/03/07(金) 19:06:17お前がそういうコード書けばいいんじゃ
0063132人目の素数さん
2008/03/09(日) 01:47:43ニューラルネットワークの単純なBPなら10行以内でコーディングできる。
対してガンダムは10進BASICだからかなり長くなる。
0064132人目の素数さん
2008/03/12(水) 12:01:18一つのグラフのyデータは100〜200程度で、もう一つは0.1〜1.5のオーダーですが、
見やすく一枚で表示したいです
0065132人目の素数さん
2008/03/16(日) 09:50:28フリーで良いのが出たら乗り替えたいけど、ここまで整備するのは大変そう。
0066132人目の素数さん
2008/03/16(日) 18:10:010067132人目の素数さん
2008/03/24(月) 20:52:35例えば、
s = Play[Sin[2 \[Pi] 1000 t] Exp[-9 t], {t, 0, 0.5}]
これで音を再生しても歪みは特にないのですが、
これを
Export["aaa.wav", s]
とし、
Import["aaa.wav"]
すると、えらく歪んでいる。どうなってるのでしょう?私だけ?
0068132人目の素数さん
2008/03/28(金) 01:15:22もうExp^(-)じゃない変なピークになってる
0069mahworld13
2008/03/30(日) 11:25:02"Think of Laura" by Christopher Cross
http://www.youtube.com/watch?v=NKuzZ7wIShw 02:42
http://www.youtube.com/watch?v=nYEWeBarLww 03:02 Live 1998
http://www.youtube.com/watch?v=0OTFLgW_9Qg 04:16 Live @Vancouver 2007/08/29
(1983 hit dedicated to the late teenager Laura)
0070132人目の素数さん
2008/04/03(木) 00:10:09もsutudent版だった
正真正銘な がくセイ! sutudent の 雛 だから じゃ 姉ぇーの
(暫く ロハに近いお値段で 酷使 しやがル)
0071132人目の素数さん
2008/04/03(木) 00:23:11との こと。
ったくアメリカ合衆国の公用語に通じていないし
且つ 日本語も覚束ない
の条件達を 完璧に満たす にも 拘らず
Mathematica を 学セイの分際で 酷使する
ヤツ_______,________,________ を 俎上に
0072132人目の素数さん
2008/04/03(木) 10:30:03mathematicaっていつから暗号作成機能ついたの?
00731stVirtue ◆.NHnubyYck
2008/04/03(木) 12:33:450074132人目の素数さん
2008/04/03(木) 12:54:14007572
2008/04/03(木) 20:59:31たうえのたきぶんしきょうがもきぬじばたけぬきしたきてまぬすきたよ
0076132人目の素数さん
2008/04/04(金) 02:16:32In[55]:= key = ToCharacterCode["ルイズかわいいよルイズ"]
Out[55]= {12523, 12452, 12474, 12363, 12431, 12356, 12356, 12424, \
12523, 12452, 12474}
(*一応鍵の長さを確認*)
In[66]:= Length[key]
Out[66]= 11
(*よって2バイト文字11文字までに対して情報量的安全である、続いて平文*)
plain = ToCharacterCode["ルイズの姉も可愛いよぉ"]
Out[58]= {12523, 12452, 12474, 12398, 22985, 12418, 21487, 24859, \
12356, 12424, 12361}
(*一応平文の長さを確認*)
In[64]:= Length[plain]
Out[64]= 11
(*おk*)
0077132人目の素数さん
2008/04/04(金) 02:17:47In[60]:= cipher = BitXor[plain, key ]
Out[60]= {0, 0, 0, 37, 26950, 198, 25515, 20883, 175, 44, 243}
In[61]:= FromCharacterCode[cipher]
Out[61]= "\.00\.00\.00%楆Æ掫冓¯,ó"
(*一応わけのわからない文字列になった*)
(*復号化*)
In[67]:= plain2 = BitXor[cipher, key ]
Out[67]= {12523, 12452, 12474, 12398, 22985, 12418, 21487, 24859, \
12356, 12424, 12361}
In[68]:= FromCharacterCode[plain2]
Out[68]= "ルイズの姉も可愛いよぉ"
(*一応復活しているのがわかる*)
実装から見ればかなりムチャクチャだけど あくまで一例ということで
0078132人目の素数さん
2008/04/05(土) 23:39:39でもwolframのオンラインショップでは2500ドルくらいで買える見たいじゃないですか。
なんだか日本人は日本の代理店経由でないと買えないみたいですが、独占禁止法に違反してませんか?
0079132人目の素数さん
2008/04/06(日) 09:41:13なお、認証その他に不都合が発生したとしても、当社としては責任を負いかねます」
0080132人目の素数さん
2008/04/06(日) 19:43:180081132人目の素数さん
2008/04/07(月) 02:21:37鍵genはver4のがそのまま使えるw
0082132人目の素数さん
2008/04/07(月) 09:05:194000−2500=1500ドルもサポート料と称してとるわけ?
ボッタクリバーみたいだね。
0083132人目の素数さん
2008/04/07(月) 18:41:35btってブータンのことですか?
0084132人目の素数さん
2008/04/09(水) 03:35:09Resolve[ForAll[p, p > 0,
Exists[q, q > 0,
ForAll[x, x >= 0 && 0 < Abs[x - a] < q,
Abs[Sqrt[x] - A] < p]]], {a, A}, Reals]
Resolve[ForAll[p, p > 0,
Exists[q, q > 0,
ForAll[a, a > 0,
ForAll[x, x > 0 && Abs[x - a] < q, Abs[1/x - 1/a] < p]]]]]
0085132人目の素数さん
2008/05/04(日) 03:20:52データをベクトルとして認識させるのも厄介だし。
Newton-Raphson法も簡単にプログラムが書けない。
EM algorithm にいたっては、てんで弱い。
プログラム形式でalgorithmを書くなら、Matlabの方が向いているのか?
0086132人目の素数さん
2008/05/04(日) 23:01:25それはそうと
データーをベクトルとして認識させるのが面倒って?
それにベクトルの記法にしたって{{},・・・・・}たかがこれだけじゃない。 Importは知ってる?
Newton-Raphson法は書かかなくくてもFindRootで実装されている。
単にアルゴリズム書きたいならcかFortranでいいんじゃないすか?
前提がフリー環境だから文献も多そうだし
0087132人目の素数さん
2008/05/05(月) 00:26:54厄介きわまりないってところかな。
テイラー展開と微分とが交互になると、FindRootでは解に行かれないんだわ。
0088132人目の素数さん
2008/05/05(月) 13:46:060089132人目の素数さん
2008/05/11(日) 07:57:49なんで普及しないんや?
値段とアクチが原因か
0090132人目の素数さん
2008/05/11(日) 09:49:260091132人目の素数さん
2008/05/13(火) 22:31:160092132人目の素数さん
2008/05/17(土) 19:29:370093132人目の素数さん
2008/05/26(月) 03:27:14これを表示させないためにはどうすればいいでしょうか?
0094132人目の素数さん
2008/05/27(火) 17:26:20調べても、式を評価しますと書いてあるだけで評価するとは何なのかまったくわかりません
具体的に教えていただければ幸いです
0095132人目の素数さん
2008/05/27(火) 17:56:31Holdを外す
0096132人目の素数さん
2008/05/29(木) 22:31:20初心者です。
i*jの行列リストの各行の任意の成分{i1,i2}をつかって
2次元グラフをプロットをしたいのですが、方法がわかりません。
ListPlotで作図したいですが教えてもらえないでしょうか。
0097132人目の素数さん
2008/05/30(金) 03:22:580098132人目の素数さん
2008/05/30(金) 15:44:53使用環境くらいちゃんと書け
0099132人目の素数さん
2008/05/31(土) 20:06:55行列の任意の要素を抽出
http://reference.wolfram.com/mathematica/ref/Part.html
平面のプロット(ListPlot3D)
http://reference.wolfram.com/mathematica/ref/ListPlot3D.html?q=ListPlot3D
0100132人目の素数さん
2008/06/02(月) 21:00:480101132人目の素数さん
2008/06/03(火) 00:03:310102132人目の素数さん
2008/06/03(火) 11:43:45写真とるか絵師が天才的に描くかしかない以上
数式処理ソフトでアレを再現したらノーベル生理学賞は間違いないと思う
0103132人目の素数さん
2008/06/05(木) 16:04:50;
listcoutourpot:ntri ; データが矛盾する三角形分割を生成しています.データを摂動して有効にすることができます.
;
というエラーが出てプロットで来ません.
plotrange->{full,full,automatic},
datarange->automatic
でプロットしています.
どうしたらいいでしょうか.
0104132人目の素数さん
2008/06/05(木) 19:16:02n=1とか。
0105132人目の素数さん
2008/06/05(木) 19:21:51http://amatias.cocolog-nifty.com/diary/2003/12/post_14.html
0106103
2008/06/06(金) 11:34:171つの列の要素の末尾すべてにバッククォーテーションがついています。
これが問題なのでしょうか。
バージョンは6.0です
5.0で実行すると[2*2以上の長方形行列ではありません」とでます
0107132人目の素数さん
2008/06/06(金) 18:23:230108132人目の素数さん
2008/06/07(土) 21:20:19Maxima以外に有名なソフトはありますか?
0109132人目の素数さん
2008/06/07(土) 22:20:34http://science6.2ch.net/test/read.cgi/math/1051275230/
無料ならMapleかなー。
0110132人目の素数さん
2008/06/07(土) 22:55:58MAPLEは30万以上するだろ
0111132人目の素数さん
2008/06/07(土) 23:54:56違法なの覚悟で言えば、
Mathematicaのクラックバージョンが最強だと思う。
0112132人目の素数さん
2008/06/08(日) 00:01:38なんてこのスレで言ってみる。
0113132人目の素数さん
2008/06/08(日) 00:47:540114132人目の素数さん
2008/06/08(日) 02:39:46因みに、
stackmatrix は行列 A、Bを縦に積んで新たな行列
A
B
を、
concatは行列A,Bを横に積んで新たな行列
[AB]
を作ります
0115132人目の素数さん
2008/06/08(日) 07:04:53<<LinearAlgebra`MatrixManipulation`
AppendColumns[A,B]
AppendRows[A,B]
0116132人目の素数さん
2008/06/08(日) 19:57:07あるいは、下でも可。
(http://reference.wolfram.com/mathematica/ref/Join.html)
行追加:Join[A,B]
とか、
列追加:Transpose[Join[Transpose[A],Transpose[B]]]
深さレベル指定すれば、Joinだけでもできる?試してない。
0117132人目の素数さん
2008/06/08(日) 23:54:08マジか…
買おうと思ってたんだが止めようかな
0118132人目の素数さん
2008/06/10(火) 14:49:002000年頃、MuPADが学生や研究者には無料で使えて、
簡易版Mathematicaという感じで使ってた。
が、今は有料化へ。30日間の無料試用版があるとはいえ、何だか悲しい。
唯一の希望はMaxima、ライセンス上、常に無料であることと、
LISP(Common Lisp)・Emacs、そしてTEXによる美しい表示が可能なので、
使いこなせれば楽しいのではないでしょうか。
0119132人目の素数さん
2008/06/10(火) 22:13:44GPLなら出力物は著作権表示とかなしで自由に使えるんだよね?
bisonとかいう例外は除けば
0120132人目の素数さん
2008/06/11(水) 18:59:590121132人目の素数さん
2008/06/14(土) 14:54:50学校の課題で以下の問題をMathmaticaでやらなければいけないんですが、
やりかったがわからないので教えて下さい
2次方程式ax^2+bx+cの近似解を、解の公式を用いて求めよ
ただし、ひとつの入力セルのみ用いるので入力セルの1〜3行目で各係数a、b、cへの代入を行い、4行目以降から解の公式を用いて結果を出力するように作成する
各係数は自由に設定してよい。NSolveやFactorを用いてはいけない。
0122132人目の素数さん
2008/06/14(土) 18:27:460123132人目の素数さん
2008/06/14(土) 21:49:530124132人目の素数さん
2008/06/14(土) 22:35:13これ30秒かかるか?・・・・瞬殺だろう
0125132人目の素数さん
2008/06/15(日) 03:19:23MathematicaのHelpは使いやすいので、使い方を覚えるといい。
まず、何のコマンドを使うかを知るのが、一番難しい。
ver5.2? ver6.0? 学校なのでまだ5.2かな、5.2ならHelpを立ち上げて、例えば「方程式」とか思いつく単語を入れてみる。
すると、SolveとかNSolveとか、それっぽいものが出てくる。
例えばSolveを見ると、使い方を言葉でごちゃごちゃ説明しているが、一番わかりやすいのは一番下にある「他の例」。
「他の例」を開くと、そのいろいろな使い方の例が挙げられている。
いろいろな例を見ていくと、値の代入の仕方とか、関数の使い方もわかってくるはず。
「NSolve, Factorを用いてはいけない。」が、これもヒント。
NSolve, FactorのHelpの中には、関連する命令についても記載されているので、それらをクリックしてみればいい。
0126132人目の素数さん
2008/06/15(日) 08:07:010127132人目の素数さん
2008/06/15(日) 09:23:14それをそのまま書けばいいじゃん。
0128132人目の素数さん
2008/06/15(日) 13:36:53�dクス
ちなみに、ver6.0
0129132人目の素数さん
2008/06/16(月) 00:31:14厳密解じゃないのか・・・・・・?
0130132人目の素数さん
2008/06/16(月) 12:54:460131132人目の素数さん
2008/06/16(月) 21:48:370132132人目の素数さん
2008/06/17(火) 20:55:190133132人目の素数さん
2008/06/18(水) 18:41:51updater 上げ
0134132人目の素数さん
2008/06/20(金) 21:10:280135132人目の素数さん
2008/06/20(金) 22:16:27Mathmaticaなら板違いです。他逝ってください。
0136132人目の素数さん
2008/06/21(土) 09:35:14風船から、亜鈴の形になり、また亜鈴から風船になる。
これを繰り返すようなアニメーションは、どうしたらできるでしょうか?
<<Graphics`ParametricPlot3D`
For[a1=0,a1£2,a1++;
r0=1;
ParametricPlot[{{r1=r0+a1*LegendreP[2,Cos[q]];r1 Cos[q],r1 Sin[q]}},{q,0,2Pi},AspectRatio®Automatic];
SphericalPlot3D[r0+Re[a1]*LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint®{-3.382,0.096,0.076},ViewVertical®{0,1,0},ImageSize®300]
];
0137132人目の素数さん
2008/06/21(土) 09:44:23すいません。文字化けしました。下記の通り修正します。
こんにちは、
風船から、亜鈴の形になり、また亜鈴から風船になる。
これを繰り返すようなアニメーションは、どうしたらできるでしょうか?
<<Graphics`ParametricPlot3D`
For[a1=0,a1<2,a1++;
r0=1;
ParametricPlot[{{r1=r0+a1*LegendreP[2,Cos[q]];r1 Cos[q],r1 Sin[q]}},{q,0,2Pi},AspectRatio->Automatic];
SphericalPlot3D[r0+Re[a1]*LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint->{-3.382,0.096,0.076},ViewVertical->{0,1,0},ImageSize®300]
];
0138132人目の素数さん
2008/06/21(土) 09:48:37すいません。やっぱりまだ文字化けしてました。下記の通り修正します。
こんにちは、
風船から、亜鈴の形になり、また亜鈴から風船になる。
これを繰り返すようなアニメーションは、どうしたらできるでしょうか?
<<Graphics`ParametricPlot3D`
For[a1=0,a1<2,a1++;
r0=1;
ParametricPlot[{{r1=r0+a1*LegendreP[2,Cos[q]];r1 Cos[q],r1 Sin[q]}},{q,0,2Pi},AspectRatio->Automatic];
SphericalPlot3D[r0+Re[a1]*LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint->{-3.382,0.096,0.076},ViewVertical->{0,1,0},ImageSize->300]
];
0139132人目の素数さん
2008/06/21(土) 13:50:350140132人目の素数さん
2008/06/22(日) 19:38:23a1++の部分を、a1+=0.1に。
ParametricPlotか、SphericalPlot3Dかのどちらかにする。
実行すると20枚のグラフィックスができるので、
右端のひとつ内側の”]”部分をダブルクリックすると、ひとつに統合される。
統合された画像をダブルクリックすると、アニメーションが動き出す。
0141132人目の素数さん
2008/06/25(水) 23:48:59100%正しく素因数分解出来ていると考えて良いので
しょうか?約20桁の整数の素因数分解を行っています。
ttp://reference.wolfram.com/mathematica/ref/FactorInteger.html
PrimeQの説明の所に、10^16以上の場合は誤動作の可能性が指摘
されているので、少々気になっています。
ttp://reference.wolfram.com/mathematica/note/SomeNotesOnInternalImplementation.html#6849
0142132人目の素数さん
2008/06/26(木) 09:13:05Maple使え
0143132人目の素数さん
2008/06/27(金) 14:28:510144132人目の素数さん
2008/06/27(金) 18:44:42Mathematica 4.2 はVistaでも動きますか?
0145132人目の素数さん
2008/06/27(金) 19:06:42動くと思うかね
やってみたまえ
0146132人目の素数さん
2008/06/28(土) 00:41:360147132人目の素数さん
2008/06/28(土) 02:41:55しかし、
まんどくせ
正規ユーザなのに割れ
0148132人目の素数さん
2008/06/28(土) 15:37:47使える。
0149132人目の素数さん
2008/06/29(日) 23:55:57積を作って比較するとよい。
0150132人目の素数さん
2008/06/30(月) 01:18:05Mapleとくらべると使えないソフトだ。
0151132人目の素数さん
2008/07/05(土) 18:16:570152132人目の素数さん
2008/07/09(水) 00:06:24というかソフトって一回出来たらもうすることあんまりないんじゃないの?
マスマティカプログラマってバージョン1の頃より減ってるんでしょうか?
エクセルだって見た目は変わったけど実質はほとんど変わってないじゃん。
セルに2変数関数の数字入れてオートフィルして3Dグラフ表示させると
エクセル97のころよりは結構改善されたとは思うけど、
この程度のバージョンアップが10年かけてやる仕事なのかなあと思うし。
エクセルの本質は全く同じと言ってもいい。
3D曲面の表示だけはまあまあよくなったと言える程度。
どうなんでしょうか。
でもエクセルの3D曲面ってマウスでグリグリ回転できないよね。
マスマティカプログラマがエクセルプログラマに転向したら
すぐに実現するのにと思うのですが。
というか、エクセル作れるぐらいのスキルならマスマティカレベルの曲面の
プログラムぐらい作れるのではないんでしょうか?
どうなのでしょうか。
0153132人目の素数さん
2008/07/09(水) 03:55:57というかソフトって一回出来たらもうすることあんまりないんじゃないの?
マスマティカプログラマってバージョン1の頃より減ってるんでしょうか?
エクセルだって見た目は変わったけど実質はほとんど変わってないじゃん。
セルに2変数関数の数字入れてオートフィルして3Dグラフ表示させると
エクセル97のころよりは結構改善されたとは思うけど、
この程度のバージョンアップが10年かけてやる仕事なのかなあと思うし。
エクセルの本質は全く同じと言ってもいい。
3D曲面の表示だけはまあまあよくなったと言える程度。
どうなんでしょうか。
でもエクセルの3D曲面ってマウスでグリグリ回転できないよね。
マスマティカプログラマがエクセルプログラマに転向したら
すぐに実現するのにと思うのですが。
というか、エクセル作れるぐらいのスキルならマスマティカレベルの曲面の
プログラムぐらい作れるのではないんでしょうか?
どうなのでしょうか。
0154132人目の素数さん
2008/07/09(水) 04:22:05ヒント: π, e
0155132人目の素数さん
2008/07/11(金) 12:32:250156132人目の素数さん
2008/07/12(土) 21:12:28一行でかけるけどな・・・
NewtonRaphson[f_, x0_] := FixedPoint[# - f[#]/f'[#] &, x0]
収束判定をコントロールしたいならこうとか:
NewtonRaphson3[f_, x0_] := NestWhile[# - f[#]/f'[#] &,
x0, Abs[#1 - #2] > 0.00001 &, 2]
0157132人目の素数さん
2008/07/13(日) 11:10:020158132人目の素数さん
2008/07/13(日) 11:51:210159132人目の素数さん
2008/07/13(日) 12:16:440160132人目の素数さん
2008/07/13(日) 14:35:200161132人目の素数さん
2008/07/13(日) 19:39:420162132人目の素数さん
2008/07/15(火) 16:40:58れ。整数i は最初の行にi=4; のように与えること。
これ教えてくださいお願いします
0163132人目の素数さん
2008/07/15(火) 19:32:280164132人目の素数さん
2008/07/16(水) 07:49:52数学系の関数ライブラリが有り得ないほど充実してるプログラム言語だと思えばJavaなんぞに負ける道理はないはず
0165132人目の素数さん
2008/07/16(水) 09:44:530166132人目の素数さん
2008/07/16(水) 10:15:10整数iを1つ1つ判定していけばOK
0167132人目の素数さん
2008/07/17(木) 14:17:40MathematicaやMapleは基本的に厳密解を求めたり数式処理をするためのソフト
高速な数値解とかシミュレーションするならMatlabでいい。
C/C++を使った数値演算シミュレーションは実装(たとえばハードウェアに載せる)まで見通した場合とか
極限まで高速性を求めたいとかでないと意味がない
Javaによる数値計算なんて何の意味もない
0168132人目の素数さん
2008/07/24(木) 16:10:38現在Mathematica5.2を使っています。
A= {{1/2, 1}, {0, 1/2}}という正方行列がある時
k = 1, 2, … に対して
A^kの1、2、無限大ノルムを求めたいです。
1、無限大のノルムについては、Norm関数を用いて求めることができました。
Norm[MatrixPower[A, k], n] (n = 1, ∞)
(もしくはTotal[MatrixPower[A, k], 2])
しかし、l2ノルムに関しては
"Norm::symb: Normは記号的要素を持つ行列を受け取りました"
と表示され、求めることが出来ません。
また、このエラーに関してのノートは記述されていませんでした。
もしご存知の方がいらしたら、具体的な求め方を教えて頂けないでしょうか。
また別の質問になりますが、以前にMathematicaを用いた際に
たとえば変数xを用いる時に、"xは正の整数である"といった制約条件を
xに対して付与出来たような記憶があるのですが
本当にできたのか、どんな方法だったか、などを失念してしまいました。
併せて教えて頂ければ幸いです。
長文な上に未熟者の質問で申し訳ありませんが
どうかよろしくお願いいたします。
0169132人目の素数さん
2008/08/01(金) 10:04:20A' を Aの共役複素数求めて転置したもの として
(A'.A)の固有値の中から最大の引っ張ってくればいいんだから
↑の手順そのまま式にして使えばいいんじゃね?
線形もMathematicaもそこまで詳しくないから参考程度に。
0170132人目の素数さん
2008/08/01(金) 11:11:33後半は Assuming などでググれ
0171132人目の素数さん
2008/08/01(金) 19:46:420172168
2008/08/02(土) 02:40:45>>170
参考になりました。ありがとうございます。
このAssumingという関数で定義したからといっても
必ずしも解を確実に判定できるという訳ではないようですね。
やはり最後は人間の判断なのでしょうか・・・。
L2ノルムに関してですが
あれから色々調べてみたところ、MathWorldにて
Matrix Norm に関しては Frobenius Norm と Euclidean norm は
同値であるとの記述をみつけました。
Mathematica では Frobenius Norm に関してはサポートされているようですので
とりあえずはこちらを用いてみたいと思います。
・・・とはいえ、kに実数を代入して Forbenius Norm と Euclidean Norm を比較しても
何故か同値にならないので、一体どこまで信用できるのかは分かりませんが・・・
(私の用い方に問題がある可能性が高いですが)
0173132人目の素数さん
2008/08/02(土) 09:12:32あなたの Frobenius Norm と Euclidean Norm の定義を書いてみてくれるかな?
0174132人目の素数さん
2008/08/05(火) 01:29:300175132人目の素数さん
2008/08/05(火) 06:08:48できるんじゃね
0176132人目の素数さん
2008/08/06(水) 00:57:06Map使えばできそう。
http://reference.wolfram.com/mathematica/ref/Map.html
導関数はD。
http://reference.wolfram.com/mathematica/ref/D.html
D[ ] のスコープにそれらしいのがある。
0177132人目の素数さん
2008/08/08(金) 01:20:17Mathematicaが組み込み関数をどう計算してるのかわからないので
今日も全てMathematicaさん任せなのです
0178132人目の素数さん
2008/08/12(火) 04:18:350179132人目の素数さん
2008/08/12(火) 06:28:04コピペをどのような形式にしたらよいですか?
簡単な例も教えていただけると助かります。
0180132人目の素数さん
2008/08/12(火) 21:43:570181132人目の素数さん
2008/08/12(火) 22:39:44変換することは出来た筈。ただし、うまく変換できるかどうかは知らない。
0182132人目の素数さん
2008/08/13(水) 18:36:50math XMLとか
0183132人目の素数さん
2008/08/14(木) 05:01:50gooで質問したのですがレスが全くなく、困ってます。
http://oshiete1.goo.ne.jp/qa4245530.html
よろしくお願いいたします。
0184132人目の素数さん
2008/08/14(木) 13:55:07z = 4 とか z = -4 とかってなってる時の x と y の値は何になってんの?
0185132人目の素数さん
2008/08/15(金) 23:30:18しらね。
そんな事を調べる機能は無い。
でも10年前のマスマティカでもこんな間違いはなかったはずだが。。。
プログラムのことよく知らないんだけど、
3D表示機能なんてクラス化されて再利用できるんじゃないの?
考えられるバグの理由としてはどう思う?
プロットした点同士を結んでいけば1〜ー1からはみ出るはずがないんだが。。。
どんなアルゴリズムで書いたらこんなグラフになるんだよ。
というかブログに画像アップしといたんだけど、ブログへのリンクが消されてた(ーー;
0186132人目の素数さん
2008/08/15(金) 23:36:16http://blog.goo.ne.jp/frozenmelon
ところで、3Dのアルゴリズムってそんなに難しいの?
0187132人目の素数さん
2008/08/16(土) 01:06:440188132人目の素数さん
2008/08/16(土) 03:30:020189132人目の素数さん
2008/08/16(土) 05:29:50Mathematicaでやりゃいいじゃん。
ついでに虚数について勉強してみれば?はみ出ることもあると思うぜ。
今回のがそのせいかは知らんけど。
0190132人目の素数さん
2008/08/16(土) 09:19:49実質的にはマスマティカに関係のある問題提起をしているのだよ。
それはplot3D関数だ。
この関数を作るのは難しいのかと言う問題提起がそれである。
0191132人目の素数さん
2008/08/16(土) 09:49:11ゴクリ・・・。とっても・・・奥が・・・深いです・・・。
0192132人目の素数さん
2008/08/16(土) 23:27:35まともに動かずpremier supportに質問したらsuse 10.3はサポートしてないと
言われた。どのdistributionをサポートしてるのかと聞いたら
http://www.wolfram.com/products/mathematica/platforms/
を見ろといわれたが、「kernel 2.4 以上の主要なもの」とあるだけ。
(実はこれはすでに見ていて、だったらsuse 10.3はOKだろうと思ってた)
誰か「このdistributionでちゃんと動いてる」ての知ってたらおせーて。
0193132人目の素数さん
2008/08/17(日) 01:50:420194132人目の素数さん
2008/08/22(金) 13:00:010195132人目の素数さん
2008/09/10(水) 09:31:160196132人目の素数さん
2008/09/10(水) 09:43:01カージオイドやアステロイド
0197132人目の素数さん
2008/09/10(水) 10:48:01mathematicaのヘルプみると一杯霊が出てるけど、コピペじゃだめか?
196の言う霊もヘルプに出てそうだけど。
0198132人目の素数さん
2008/09/10(水) 10:50:590199132人目の素数さん
2008/09/10(水) 11:01:26処理に時間がかかる係数をつっこめば、遠目には絵にしか見えない。
0200132人目の素数さん
2008/09/10(水) 11:05:18「数学の情報を調べろ」ということが課題じゃないのかい。
人に聞いて解決するより、自分で調べた事や感じた事をレポートした方が主旨に沿う気がする。
0201132人目の素数さん
2008/09/10(水) 11:18:33http://demonstrations.wolfram.com/
みるとさらに山ほど例がでてるよ。ソースもみれそうだし。
0202132人目の素数さん
2008/09/10(水) 18:38:040203132人目の素数さん
2008/09/10(水) 19:09:29ここはドラえもんか、アンパンマンの顔を書いて
高評価をねらえ
0204132人目の素数さん
2008/09/10(水) 19:49:200205132人目の素数さん
2008/09/11(木) 10:05:50直線と円が描ければいいんだし。
0206202です
2008/09/11(木) 21:47:22ちゃんと自力で教科書読みながらやってるんだがなんかうまくできない。もぅなんかを丸写ししたい
0207132人目の素数さん
2008/09/12(金) 12:57:50痔力が理想だが、少しくらい霊におすがりしてもバチは当たらないよ。
0208202です
2008/09/12(金) 14:47:41誰かたすけて
0209132人目の素数さん
2008/09/12(金) 15:00:01気になる図形をクリックして、文中の「SHOW SOURCE CODE」や「DOWNLOAD SOURCE CODE NOTEBOOK」
英語が読めないことはいいわけにならない。
0210132人目の素数さん
2008/09/13(土) 12:20:220211132人目の素数さん
2008/09/13(土) 12:32:57んまりレベル高いの提出したら怪しまれるかな?あん
まりレベル高いの提出したら怪しまれるかな?あんま
りレベル高いの提出したら怪しまれるかな?あんまり
レベル高いの提出したら怪しまれるかな?あんまりレ
ベル高いの提出したら怪しまれるかな?あんまりレベ
ル高いの提出したら怪しまれるかな?あんまりレベル
高いの提出したら怪しまれるかな?あんまりレベル高
いの提出したら怪しまれるかな?あんまりレベル高い
の提出したら怪しまれるかな?あんまりレベル高いの
提出したら怪しまれるかな?あんまりレベル高いの提
出したら怪しまれるかな?あんまりレベル高いの提出
したら怪しまれるかな?あんまりレベル高いの提出し
たら怪しまれるかな?あんまりレベル高いの提出した
ら怪しまれるかな?あんまりレベル高いの提出したら
怪しまれるかな?あんまりレベル高いの提出したら怪
しまれるかな?あんまりレベル高いの提出したら怪し
まれるかな?あんまりレベル高いの提出したら怪しま
れるかな?あんまりレベル高いの提出したら怪しまれ
るかな?あんまりレベル高いの提出したら怪しまれる
かな?あんまりレベル高いの提出したら怪しまれるか
な?あんまりレベル高いの提出したら怪しまれるかな
?あんまりレベル高いの提出したら怪しまれるかな?
どこのヘボ大だよ
0212202です
2008/09/14(日) 13:35:33ここにはっていただけませんか?かなりがんばってるんですができなくて・・・
頼みますm(__)m
0213132人目の素数さん
2008/09/14(日) 14:41:08で検索しろ
0214132人目の素数さん
2008/09/14(日) 15:10:300215132人目の素数さん
2008/09/14(日) 22:51:22Importでファイルを読み込む際に、01.dat,02.datというような
連番のファイルを順々に読み込ませたいのですが
いい方法はないでしょうか。
ファイル名のパスを文字列として収納すればいいと思うのですが、
その方法がよくわかりません。
0216132人目の素数さん
2008/09/15(月) 00:56:56data = Import["http://d.hatena.ne.jp/images/diary/s/survival/2004-11-21.jpg"];
Show[data];
0217132人目の素数さん
2008/09/15(月) 03:46:03FromCharacterCode[{12489, 12521, 12360, 12418, 12435}]
0218132人目の素数さん
2008/09/15(月) 08:06:24>ファイル名のパスを文字列として収納すればいいと思うのですが、
FileNames[ DigitCharacter .. ~~ ".dat"]
はどう?
0219132人目の素数さん
2008/09/19(金) 20:29:350220132人目の素数さん
2008/09/19(金) 21:22:130221132人目の素数さん
2008/09/19(金) 23:51:04StringFormをToStringするんだ!
0222132人目の素数さん
2008/10/06(月) 23:21:52.csv から、複数列の表データを読み込んだとして、その中のいくつかの
列を入れ替えたりしたいのですが、どうするのが一番いいんでしょうか。
0223132人目の素数さん
2008/10/10(金) 22:36:220224132人目の素数さん
2008/10/11(土) 23:58:57Sort
http://reference.wolfram.com/mathematica/ref/Sort.html?q=Sort&lang=en
行列操作
http://reference.wolfram.com/mathematica/guide/PartsOfMatrices.html
http://reference.wolfram.com/mathematica/tutorial/AddingRemovingAndModifyingListElements.html
0225132人目の素数さん
2008/10/12(日) 00:06:23ToStringのアプリケーション2
http://reference.wolfram.com/mathematica/ref/ToString.html
0226132人目の素数さん
2008/10/12(日) 13:56:246のドキュメントセンターって記号の検索できなくね!?
例えば「/.」とか「:>」とかってみんな検索できる?
5のヘルプブラウザはできるんだが
0227132人目の素数さん
2008/10/13(月) 22:37:190228132人目の素数さん
2008/10/27(月) 20:59:060229132人目の素数さん
2008/10/28(火) 18:53:40はVer 4で十分だ
安く売らんかい!
0230132人目の素数さん
2008/10/29(水) 18:14:55俺は Ver3.0 だ。文句あるか。
0231132人目の素数さん
2008/10/29(水) 23:09:38御気の毒としか申しようがありません
ほんとうに御愁傷様です
0232132人目の素数さん
2008/10/30(木) 18:24:25で、君は Mathematica をどんな風に使ってるのかね。
ver3.0 では「御愁傷様」だという例を一つで良いから見せてみろ。
0233132人目の素数さん
2008/10/30(木) 19:59:47よく落ちた
0234132人目の素数さん
2008/10/30(木) 21:10:25承認しましゅって言われた
0235132人目の素数さん
2008/11/01(土) 22:12:56おめ。
何系の会社?
どこの販売代理店から買うの?
MapleやMatlabじゃなくてMathematicaの理由は?
0236132人目の素数さん
2008/11/02(日) 17:08:142週間たってからXPの復元機能つかってインストール前に戻せる?
0237132人目の素数さん
2008/11/02(日) 17:09:350238132人目の素数さん
2008/11/02(日) 17:40:130239132人目の素数さん
2008/11/02(日) 18:28:46mathematicaもmatlabも両方使ってる身としては
これらを比較する感覚が分からない
0240132人目の素数さん
2008/11/02(日) 22:57:130241132人目の素数さん
2008/11/03(月) 17:32:20Matlabで数式処理もできる。結局どっちのウェイトを重視
するかだよ。
0242132人目の素数さん
2008/11/04(火) 00:25:12数学5強ソフト:MA5
0243132人目の素数さん
2008/11/04(火) 22:48:11ノートブックでの編集時に表示形で表示させる方法を教えて頂けませんでしょう
か。ギリシャ文字をパレットより入力してしばらくは表示名で表示されていたのですが、なにかの
拍子に完全名での表示となってしまいました。形式変換や表示形式をいじくってみたのですがうま
くいきません。
0244132人目の素数さん
2008/11/06(木) 10:10:17Ctrl+Shift+tのことか?
0245132人目の素数さん
2008/11/10(月) 10:00:04ParametricPlot3Dなどで描いた曲面をイラストレーターに貼り付けたいのですが、
ただコピペをしただけでは画像として貼り付けられてしまい、
一度epsファイルで保存した後、イラストレーターで読み込まないと
線の集まり(言い方が変ですみません)として読み込めません。
Mathematica5の時は出来ました。
形式を選択してコピーの項目をいろいろ試したのですがうまくいきません。
epsファイルに保存しないでコピペをする方法を
ご存じの方がいらっしゃいましたら教えてください。
0246132人目の素数さん
2008/11/10(月) 17:36:15試さずに書いてるが
0247132人目の素数さん
2008/11/10(月) 17:37:08正規に効いてみたら
メーリングリストとか
0248132人目の素数さん
2008/11/22(土) 21:17:26, - ' ,..、、.ヾ{{フ'⌒`ヽ、
/ ,:', -‐‐` ´ '´⌒ヽ ヾ:、
. ,' ,'´ ,ィ ,ィ ,' , `ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,))
! | ! .,'-.{ ! !|; |`、.}゙!.! |. ! ヽ.
', ', |Vァ=、゙、 `゙、!-_:ト,リ', l ! | ゙',
ヽ、', l:!Kノ}. f:_.)i゙i: リ ! l ル
| l!iヾ- ' , .!__:ノ ゙ ,リ l リ'´
. ',|!!、 r‐┐ ` ノ'. /,イ / ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄
'i!゙、ヽ、 ゙ー' _, ィ,:',:''´ < Mathematica 7リリースおめ!
゙:、ィ、jヾー::: 'iヘ .ノ',リ. \_____________
,、- '´ ヽ、゙、 { `>"、
/\\ ', } //`ヽ
0249132人目の素数さん
2008/11/22(土) 22:46:390250132人目の素数さん
2008/11/23(日) 05:39:44☆ チン 〃 ∧_∧ / ̄ ̄ ̄ ̄ ̄ ̄ ̄
ヽ ___\(\・∀・)< 6.1まだー?
\_/⊂ ⊂_)_ \_______
/ ̄ ̄ ̄ ̄ ̄ ̄ ̄/|
|  ̄  ̄ ̄ ̄ ̄ ̄ ̄:| :|
| .|/
0251132人目の素数さん
2008/11/23(日) 23:44:17俺的には4の虹色のアイコンが傑作だったと思う
今回のやつは完全にウイルスだ
ところで7ってwhat's newはあるけど(http://www.wolfram.com/products/mathematica/newin7/)
まだ製品はリリースされてなくね!?
0252132人目の素数さん
2008/11/25(火) 17:17:15太くしたいんですが、どうしたらいいか知りませんか?
0253132人目の素数さん
2008/11/25(火) 21:19:03だろ
ヘルプみりゃ
太くできるんじゃね
0254132人目の素数さん
2008/11/25(火) 22:32:56太くするのは見当たらないんですが。
typoすまそ
0255132人目の素数さん
2008/11/25(火) 22:35:21描いた図をEPSでLaTeXに張り込むと
線が細すぎてきれいに表示されないので困ってます。
0256132人目の素数さん
2008/11/26(水) 18:25:590257132人目の素数さん
2008/11/26(水) 19:02:210258132人目の素数さん
2008/12/03(水) 20:38:52{1,1,-1,-1},
{I a ,-I a, -b, b},
{e^(I ka), e^(-I ka),-e^(I kd)e^(-kb),-e^(I kd)e^(kb)},
{I a e^(I ka), -I a e^(-I ka),- b e^(I kd)e^(-kb),b e^(I kd)e^(kb)}
の行列式を
A[[1,1]] Det[ A[[{2,3,4},{2,3,4}]] ]
+A[[1,2]] Det[ A[[{2,3,4},{3,4,1}]] ]
+A[[1,3]] Det[ A[[{2,3,4},{4,1,2}]] ]
+A[[1,4]] Det[ A[[{2,3,4},{1,2,3}]] ]
で計算しました。答えは出たのですが、最終的には
Det[A] = 8I a b e^{I kd} {Cos[kd] - Cos[ka] Cosh[kb] - ((b^2-a^2)/2a b) Sin[ka] Sinh[kb] }
と整理できるはずなのですが、うまくいきません。
どうすればいいですか?
0259132人目の素数さん
2008/12/04(木) 09:38:14ExpToTrigとFullSimpifyでできない?
eは\[ExponentialE]を使ってる?
0260132人目の素数さん
2008/12/04(木) 12:08:04ありがとうございます。できました。
>eは\[ExponentialE]を使ってる?
・・・・使ってませんでした。これが原因か。
4次のままで行列式が計算できなかったのも orz orz orz
0261132人目の素数さん
2008/12/16(火) 18:18:00学校の課題でPrimal-Dual法を用いて最小費用流を求めるプログラムを書け
というのが出たのですが、現在"Residual Flow"の生成が上手くいかない状況です
どなたか詳しい方いましたら教えていただけないでしょうか?
なおVerは6.0です
0262132人目の素数さん
2008/12/16(火) 21:15:47うまく行かないといわれても、どういう実装をしようとしていて、
どの部分が分からないのかを言ってくれないと、さっぱりだよ。
0263261
2008/12/17(水) 13:05:12グラフを使って最小費用流問題をPrimal-Dual法で解きたいのですが
残余ネットワークの生成の仕方が分からないということです
0264132人目の素数さん
2008/12/17(水) 23:17:02それは >>262 の質問に対する答えに全くなってないよ。
「どういう実装」と聞いているのは、例えばどのようにデータを管理しているか、
とか、そういうことを聞いている(グラフをどう持っているかとか)。
「どの部分」と聞いているのは、残余ネットワークの定義が分からないのか、
残余ネットワークを数式でどう表したらいいかわからないのか、
定義も式も分かっているけど Mathematica での書き方が分からないのか、
とかそういうのを聞いている。
0265261
2008/12/18(木) 19:24:32一応今の所は下のような感じになっています
g0 = {{1, 2}, {2, 3}, {1, 3}, {1, 4}, {4, 3}, {4, 5}, {3, 5}};
g = FromOrderedPairs[g0];
h0 = Table[Random[Integer, {1, 100}], {7}]
h1 = SetEdgeWeights[g, h0];
h2 = SetEdgeLabels[g, GetEdgeWeights[h1]];
capacity = Table[Random[Integer, {2, 15}], {7}]
h3 = SetEdgeWeights[g, capacity];
costlist = {h0, ToOrderedPairs[g]};
TableForm[costlist];
capacitylist = {capacity, ToOrderedPairs[g]};
TableForm[capacitylist];
flowlist = {Table[flow, {Length[g0]}], ToOrderedPairs[g]};
TableForm[flowlist];
0266261
2008/12/18(木) 19:32:15残余ネットワークのどの部分が分からないかということに関しては
「残余ネットワークを数式でどう表したらいいかが分からない」です。
あともう一つ分からない所があるのですが、mathematicaのグラフの辺に
2つ以上の重みをつけることは可能なのでしょうか?
質問ばかりで申し訳ないですがよろしくお願いします
0267132人目の素数さん
2008/12/18(木) 22:26:05数式で書けない限りMathematicaでも書けない。
その部分は教科書なり何なりで勉強してらっしゃい。
二種類以上の重みは、「重みのペア」を一つの重みだと考えればいいよ。
0269132人目の素数さん
2008/12/20(土) 15:49:35ひとつ前でなくてひとつ上の出力を参照する記号はありますか?
ctrl+shift+Lではなくて%みたいに短いと見やすいんですが、
%で繋いでドンドン積んでいって途中で間違えるとやり直すのが面倒だよぬ
0270132人目の素数さん
2008/12/21(日) 01:56:560271132人目の素数さん
2008/12/21(日) 11:24:15http://reference.wolfram.com/mathematica/ref/Out.html
でも、何度も使うようだったら、変数に記憶しておくほうがいいよ。
0272269
2008/12/21(日) 13:31:330273132人目の素数さん
2008/12/21(日) 22:07:490274132人目の素数さん
2008/12/22(月) 16:40:24新しいセルに手で打ち直すしかないのかな?
0275132人目の素数さん
2008/12/22(月) 20:25:360276132人目の素数さん
2008/12/23(火) 13:14:390278132人目の素数さん
2008/12/23(火) 16:24:170279132人目の素数さん
2008/12/23(火) 22:43:58もちろんバグがてんこ盛りだよん!
0280132人目の素数さん
2008/12/26(金) 16:30:01関数自体を自分で書き換えるしかないでしょうか?
f[n_] :=
Nest[If[
If[# == 0, 0,
If[
Mod[#2, #] ==
0, # + #0[# - 1, #2], #0[# - 1, #2]]] &[#, #] ==
2*#, #, #0[# + 1]] &[# + 1] &, 1, #] &@n
0281132人目の素数さん
2008/12/26(金) 22:07:34まずどっからもってきた関数なのか説明しろ
0282132人目の素数さん
2008/12/26(金) 22:37:29動かしてないが見た感じ定義をそのまんま書いたっぽい
再帰が深すぎるんじゃない?
再帰の限界を環境変数でいぢればいい
0283132人目の素数さん
2008/12/27(土) 17:28:10させようとしたのですが、
Limit[Binomial[N, r](m/N)^r (1 - m/N)^(N - r), N -> Infinity]
とすると
Gamma[1/N+1+O[N]^3]で真性特異点が見つかりましたと言われて計算できません
Binomial関数を階乗になおしても同じです。
MathematicaではGamma関数を含むと極限計算ができないんですか?
それともなにか工夫をすれば上記計算ができるようになりますか?
0284132人目の素数さん
2008/12/29(月) 14:20:42ソフトがあるのはこの掲示板でも何回か見たので使ってみた
のですが、横軸の 1 と縦軸の 1 が同じ長さになってしまいます。
Mathematica に入力したのは以下のようなものです。
Get[ "・・・/ketpic6.m"]
Get[ "・・・/ketpicapp6.m"]
XMIN = 0;
XMAX = 5;
g1 = Plot[x^2, {x, XMIN, XMAX}];
setwindow[{XMIN, XMAX}, {0, 25}];
openfile["foo.tex"]
openpicture["0.5cm"];
drwline[g1];
closepicture[1];
closefile[]
・・・のところは本当は ketpic6.m, ketpicapp6.m
のディレクトリの絶対パスです。
TeX では input で読み込んでいるだけなのですが
図は縦横比が 25 : 5 で縦長になっています。
Mathematica では AspectRatio というコマンドで変えるので
ketpic6.m, ketpicapp6.m の AspectRatio->Automatic
という部分を消してみたのですが変わりませんでした。
使っている人がいたら教えてください。
0285132人目の素数さん
2008/12/29(月) 20:45:59ver6だと求まったよ
0286132人目の素数さん
2008/12/30(火) 12:34:10わざわざやる需要あんの?
0287132人目の素数さん
2008/12/30(火) 19:58:520288132人目の素数さん
2008/12/30(火) 19:59:24○ だけでも
0289132人目の素数さん
2008/12/30(火) 21:03:43さらにバグ増加です
0290132人目の素数さん
2008/12/30(火) 21:49:390291132人目の素数さん
2008/12/30(火) 21:56:510292132人目の素数さん
2008/12/31(水) 04:16:20無料でくれるらしい
0293132人目の素数さん
2008/12/31(水) 07:05:43俺もプレミア入ってるがそんな案内まだきてないぞ?
どこの代理店?日本電子計算?HULINKS?
0294132人目の素数さん
2009/01/01(木) 18:16:350296132人目の素数さん
2009/01/02(金) 02:33:28「バージョンアップできるのは二つ前までで終わりだ
二つ前はマイナーを入れて数えるんだぞ」
って奴
0297132人目の素数さん
2009/01/07(水) 19:14:52英語版から1ヶ月経たずに出すとは翻訳チームもなかなかのもの。
ただし、誤字脱字やふさわしくない翻訳も散見されるが。
んで、肝心の内容はというと、俺に言わせれば素晴らしいの一言。
きれいなグラフがサクッと作れるし、並列計算機能が標準搭載されるようになったし。
個人的にびっくりした機能が数列を何項か与えると一般項を返す関数。
数列データベースを参照してるとのことで、例えばn次元球の体積だと
5項ではだめだけど6項与えたら一般項をちゃんと求めてきた。すごい。
0298132人目の素数さん
2009/01/08(木) 20:09:040299132人目の素数さん
2009/01/08(木) 23:54:320300132人目の素数さん
2009/01/09(金) 00:25:490301132人目の素数さん
2009/01/09(金) 11:00:090302132人目の素数さん
2009/01/09(金) 18:34:390303132人目の素数さん
2009/01/09(金) 21:02:05それ専用の関数がある
http://reference.wolfram.com/mathematica/ref/SeriesCoefficient.html
0304132人目の素数さん
2009/01/09(金) 21:58:47年度末調整で買うか。
0305132人目の素数さん
2009/01/09(金) 22:58:170306132人目の素数さん
2009/01/09(金) 23:24:470307132人目の素数さん
2009/01/10(土) 02:50:26しかしアクティベーションがうざい
0308132人目の素数さん
2009/01/11(日) 00:10:350309132人目の素数さん
2009/01/11(日) 03:33:500310132人目の素数さん
2009/01/11(日) 10:50:41シングルローカルライセンスかネットワークライセンスを
聞いてこないや。
これが原因か?
0311132人目の素数さん
2009/01/11(日) 14:24:100312132人目の素数さん
2009/01/12(月) 06:07:420313132人目の素数さん
2009/01/12(月) 17:36:1520000円とか
0314132人目の素数さん
2009/01/12(月) 17:42:59社会人?
0315132人目の素数さん
2009/01/12(月) 17:50:370316132人目の素数さん
2009/01/12(月) 18:15:280317132人目の素数さん
2009/01/12(月) 18:43:250318132人目の素数さん
2009/01/12(月) 19:52:39要組合員証
0319132人目の素数さん
2009/01/14(水) 23:51:02v4-もっさり
v5-やや改善
v5.2-サクサク!
v6-もっさり&メモリ消費
v7-高機能でサクサク!
これで、5.2&6.0両方使いだったが、v7だけに移行できる。
0320132人目の素数さん
2009/01/15(木) 00:04:59週末に一念発起してインストする。
0321132人目の素数さん
2009/01/15(木) 00:39:210322132人目の素数さん
2009/01/16(金) 08:15:550323132人目の素数さん
2009/01/16(金) 17:59:32ヘルプの日本語はよくできてる
0324132人目の素数さん
2009/01/16(金) 18:01:46そりゃそうだよなw
0325132人目の素数さん
2009/01/16(金) 21:21:56落とす奴のほうが多くなったら開発なんて続かない.
落としたやつも開発者への感謝の気持ちは忘れんなよ.
0326132人目の素数さん
2009/01/16(金) 22:21:53俺もVer3持ってるけどアカデミック版だからこそ買えた。
97年だったけど初めてインターネットのありがたさを実感したわ。
それまでどこで売っているかさえわからなかったから。
0327132人目の素数さん
2009/01/16(金) 22:23:45ドングリだっけ? LW3Dのやつとか
0328132人目の素数さん
2009/01/16(金) 23:40:36あんなものはすぐクラックされる
悪血が最狂
アドビをみろ
0329132人目の素数さん
2009/01/16(金) 23:47:55v7入らないっぽい。
OS再インストールしないとだめかも。
0330132人目の素数さん
2009/01/17(土) 00:50:59autorunのファイルでトロイを感知した
拾ってきたやつだからかなw
0331132人目の素数さん
2009/01/17(土) 13:58:20どなたかご存じですか?
生協で聞けば教えて
もらえるんでしょうか?
0332132人目の素数さん
2009/01/17(土) 16:15:410333132人目の素数さん
2009/01/17(土) 17:01:220334132人目の素数さん
2009/01/17(土) 19:17:13おれは学生だっつーの!
0335132人目の素数さん
2009/01/17(土) 19:18:23Mathematica使ってるのって
思い出してみると、>>333のような
馬鹿ばっかりだもん。
0336132人目の素数さん
2009/01/17(土) 19:40:11ver6と一緒に入れているけど、うちの環境では大丈夫。
ver6+画像処理アプリを使っており、ver7で統合されたのは良いのだが、コマンド名が結構変わっている・・・
同じコマンド名のものは2重定義されているとエラー出るし、アプリ部分は結構書き換えが必要。
それでもver6+アプリよりは快適なので、この際ver7に乗り換える!
0337132人目の素数さん
2009/01/18(日) 16:19:300338132人目の素数さん
2009/01/18(日) 17:19:15わざわざMaximaなんてものを使う意味は無くなった
金出す価値がある
0339132人目の素数さん
2009/01/18(日) 18:03:54商用のほうが若干安心できるんだよね.気持ちの問題が大きいけど.
0340132人目の素数さん
2009/01/18(日) 18:19:350341132人目の素数さん
2009/01/18(日) 18:45:40画像処理アプリって何?ウルフラムが画像処理専用のアプリを出してるのか?
リンクはってくれ
>>339
研究者だったら所属する大学か組織が買ってくれないのか?
0342132人目の素数さん
2009/01/18(日) 18:52:32もちろん研究室でも買ってるよ。
0343132人目の素数さん
2009/01/18(日) 19:01:30/ / /フ | | r-、 `ー' `ー′
L__/ r"> // | | ヽヽ
,.._ / / { { || 〉〉 ,-、
// { { い | | /ノ / / ,、
{ { ヽヽ、 ヽ`'二_二' / / / / ノ
/7 ヽヽ、 `'-ニ,/´_=;;;;=_ ヽ`ニノ / /
{ { `ー-`ニニ /r'゙_ ヽ / _`',`ー-‐'´ノ r'7
ヽヽ、 _,.ノl=・=-.∴.-=・=!'-ニ二__//
「`! `'-`ニニ二-ァ | `ー'_ !」 _ー' |ヾー―---‐'´ ,-、
l ! ,. -'´ri |l ‘~`ー'~′ l| r‐`ニニ'ー- 、_//
ヽ、ヽ--―'ニ ‐'´ {ヾ ! ,r;'ニーニヽ、!l/}ヽ、 ``ー‐'´_
__ ``'' ''´ /,トミ! ir,!-┴-!、ヾ, !'/\ヽ、_ /ノ
い、 ///ハ||′-―-、`|||' l ヽ `ヽ、二./ >>61こらまんじゃ夫馬いっちょ
\`ー―'ニ´‐',/ ! ||!/ ヽl|| lヽ \ ,-、 おまんじゃえーじゃよなw
,、 ̄ ,. -'_/ ! ||!______l|| l \`ー‐'´ ノ
── = ニ /=、。。。。。。。。。。。。。。。。r=、ヽ
─ =ニ三 (◎ ヽ-─────(◎ )
ノ◎、 |\ \ / / | /◎、
(_,rへ `ソ /> ◎) (◎く| レ' ,rへ )
─ = ニ \◎'/ / \ ヽ、◎/
ノ / \ ヽ
─ =ニ三 ( ◎( ) ◎)
─ = ー、_ら ⊂、_,r´
0344132人目の素数さん
2009/01/18(日) 22:35:14どこが代理店になるんだろ?
価格上がっていくのかな?
0345132人目の素数さん
2009/01/18(日) 22:44:050346132人目の素数さん
2009/01/19(月) 01:34:36ソースは何?
ほとんどの大学はMathematica導入してるから
生協が代理店じゃなくなったら、他の代理店はウハウハだな
日本電子計算とヒューリンクスぐらいしかないだろ
0347132人目の素数さん
2009/01/19(月) 03:06:11ネットで検索してみてください。
どこも代理店契約は今年の3月で終わります。
ただ、代理店じゃなくなっても販売は
継続すると生協のHPには書いてありますけど。
やはり今後はDL販売がメインになるみたいです。
StudentバージョンもDLで購入できるようですね。
0348132人目の素数さん
2009/01/19(月) 03:22:13学生かどうかの判断はどこでするんだろう?
0349132人目の素数さん
2009/01/19(月) 12:02:070350132人目の素数さん
2009/01/20(火) 18:44:01生協の代理店契約が切れるからって、DL販売がメインになるなんてことはないだろ
そもそもDL販売なんてものに代理店なんて必要ないし
ウルフラムとの直接のやりとりなんだから
ということはやっぱり国内の大手代理店と契約することになるんじゃないか?
0351132人目の素数さん
2009/01/20(火) 22:07:240352132人目の素数さん
2009/01/20(火) 22:13:42つうか、Mathematica使える人は、購入するぐらいの英語は使えるだろ。
0353132人目の素数さん
2009/01/21(水) 13:02:490354132人目の素数さん
2009/01/21(水) 13:09:25それでも敷居が低いに越したことは無いだろう
「どうしても欲しい」人は多少のことは何とかするだろうが
「なんとなく欲しいかな」ぐらいの人にとっては障壁は低ければ低いほどいい
0355132人目の素数さん
2009/01/21(水) 20:14:17日本電子計算とヒューリンクスだけか
線香はだめだし
0356132人目の素数さん
2009/01/22(木) 03:19:33パスワードは keygen.exe があったので、これで作ると思うんですけど。
私は WAREZ 的なものはあまり詳しくないので分かりません。
これでは安心して購入できません。
教えてください。
0357132人目の素数さん
2009/01/22(木) 08:17:19ヨドバシカメラで売ってるよ。
(うちの近くのヨドバシには置いてあった)
0358132人目の素数さん
2009/01/22(木) 08:22:23うちの近くはビックカメラだった。
0359132人目の素数さん
2009/01/22(木) 18:52:12落ちているのはニセ物だよ
今回のは買わないとダメなのさ
0360132人目の素数さん
2009/01/22(木) 20:31:59http://item.rakuten.co.jp/plusyu/4988727009970/
HULINKSは定価しか乗せてないけど、見積もり取ったら幾らぐらいに
なるんだろ。
0361132人目の素数さん
2009/01/23(金) 03:38:05解けるソフト教えてください。
0362132人目の素数さん
2009/01/23(金) 04:21:04Mathematicaはある程度の範囲のものは解いてくれるけど、
それ以上は高望みしすぎだよ。
0363132人目の素数さん
2009/01/23(金) 06:13:300364132人目の素数さん
2009/01/23(金) 09:33:17やっぱそういうのだとマスマティカ最強?
0365132人目の素数さん
2009/01/23(金) 09:56:410366132人目の素数さん
2009/01/23(金) 17:18:28数値的にです。
0367132人目の素数さん
2009/01/23(金) 21:53:220368132人目の素数さん
2009/01/23(金) 23:00:30ただし計算速度とか精度にこだわるなら別のものを使うべき。
0369132人目の素数さん
2009/01/25(日) 09:56:13フリーのScilabやMaximaでも十分だよな。
0370132人目の素数さん
2009/01/25(日) 11:23:02はい。ScilabやMaximaでも十分な用途の場合は十分です。
0371132人目の素数さん
2009/01/25(日) 16:47:290372132人目の素数さん
2009/01/25(日) 17:00:190373132人目の素数さん
2009/01/26(月) 00:16:13ヘルプの使いづらさがね。
それに最近の傾向として、
開発元は純粋な数学関連の開発をおろそかにして、
金を儲けやすい工学への応用に走ってる。
トヨタも関係してるのかもしれないが、
これはもう終わりだね。
Mapleさようなら。
0374132人目の素数さん
2009/01/26(月) 08:30:09数値計算は他にもっとすごいのいくらでもあるけど
0375132人目の素数さん
2009/01/26(月) 20:44:50こういうふうに行列の微分したいんだけどできる?
0376132人目の素数さん
2009/01/26(月) 20:56:19でもごくまれにMathematicaで解けなかった常微分方程式が
HPの数式処理電卓で解けたりするんだよなw
0377132人目の素数さん
2009/01/28(水) 03:37:31できる。例としてまず行列を用意する。
mat = Table[\[ExponentialE]^(j k x), {j, 1, 4}, {k, 1, 4}]
% // MatrixForm
次にこの行列を格納したリストmatに微分演算をする
ans = \!\(\*SubscriptBox[\(\[PartialD]\), \(x\)]\ mat\)
ans // MatrixForm
微分演算ができている
0378132人目の素数さん
2009/01/29(木) 01:25:39出力をコピペしてもとても読めません。
あと計算結果から、同じ項をなにか文字で置き換えさせたりとかもできるなら知りたいです。
例
out> {Σz^3+y^2 + Σx+y^4 - Σx^3*y^1 }/{ Σx-6+y^4 + Σx+y^4 Σ・・}
みたいなのを
{Σz^3+y^2 + A - Σx^3*y^1 }/{ Σx-6+y^4 + A Σ・・}
といった感じです。
いずれも厳しそうな感じですが、よろしくおねがいします。
http://www.mathmlconference.org/2002/presentations/smirnova/
これを使う。
なんていうのでどうでしょう?
0380132人目の素数さん
2009/01/29(木) 02:01:21TeXへの貼り付けは TeXForm[...] で吐いて微調整すれば何とかなる。
Javaに貼り付けるというのは意味がよくわからんけど、
式をソースコードの一部に使うなら CForm で吐いて調整、
コメントに貼り付けるなら InputForm かね。まあこのあたりで調べると色々あるよ。
0381132人目の素数さん
2009/01/29(木) 13:36:32複数の値を覚えさせる時、関数を使うのとリストを使うのとどちらがいいのでしょうか?
0382132人目の素数さん
2009/01/29(木) 15:21:270383132人目の素数さん
2009/01/29(木) 15:29:040384132人目の素数さん
2009/01/29(木) 15:37:420385378
2009/01/29(木) 17:41:41�dくす!研究がいっきにすすみそうです。どうしても手計算じゃすまないところを
mathematicaに頼ってみたものの項が多すぎて詰んでいました。
0386132人目の素数さん
2009/01/29(木) 17:43:01∂f/∂u+∂g/∂v=0
∂i/∂u+∂j/∂v=0
∂k/∂u+∂l/∂v=0
を境界、初期条件を入れてx,y,zについて解きたい(数値解析)したいんですが、
差分してニュートンラプソン方が良いのでしょうか?
マスマティカで可能でしょうか?
フォートランとかがいいんでしょうか?
0387132人目の素数さん
2009/01/29(木) 18:27:27ふぉとらんらんらん
ふぉーと らんらんら〜ん
0388132人目の素数さん
2009/01/29(木) 18:44:51{a,b,c,d,e,f}
があってそれを
{ {a,b}, {b,c}, {c,d}, {d,e}, {e,f} }
のようにして返す関数を作りたいのですがどうすればいいでしょうか?
組み込み関数にあるかと思ったのですが見つからなくて
0389132人目の素数さん
2009/01/29(木) 18:51:470390132人目の素数さん
2009/01/29(木) 18:52:39ありがとうございます!
できました
0391132人目の素数さん
2009/01/29(木) 23:51:43data=Import[original];
data2=Drop[data,*];
data3=オリジナル関数[data2];
ていうのを手順を何回もくりかえしたいのですが、
この一連の処理をさらに関数化したいのですが、うまくできません(ていうかできない?)
どなたか教えてください。data3をrerturnしたいです。
NewData=新しい関数[orginal]のようにしたいです。
0392132人目の素数さん
2009/01/30(金) 00:13:17data=Import[original];
data2=Drop[data,*];
data3=オリジナル関数[data2]
);
みたいな感じ?
0393132人目の素数さん
2009/01/30(金) 00:52:15391です。有難うございます。
新しい関数[orginal_]:=( data3=オリジナル関数[Drop[Import[original] ,*]; ]);
みたいに入れ子にすればできるみたいです。関数宣言のなかに変数(dataとかdata2とか)使うとむりみたいですね。
仕様?
0394132人目の素数さん
2009/01/30(金) 10:39:37Partition[{a,b,c,d,e,f},2,1]
でできるだろ?
0395132人目の素数さん
2009/01/30(金) 19:50:57Block,Module,With
あたりを使え
0396132人目の素数さん
2009/01/31(土) 21:13:46DSolveでとけますかね?
解けるやつ教えてください。
0397132人目の素数さん
2009/01/31(土) 21:24:55なんでマルチしちゃったの?
0398132人目の素数さん
2009/02/01(日) 05:51:42順列 Permutation nPr: FactorialPower[n,r]
この事実に俺は日本で一番最初に気づいたので、ここにメモしておく。
Mathematica使い始めた頃を思い出すが、
ヘルプでCombinationとして探すも出てこず、
色々調べた結果、2項係数として関数になっていることがわかった。
Permutationは探せど探せど出てこず、
結局自分で定義して使ってた。
後に分かったことは、CとかPとか使ってるのは日本だけということ。
0399132人目の素数さん
2009/02/01(日) 10:08:24nCr、nPrがそれぞれ特殊なケースでしかない。
0400132人目の素数さん
2009/02/01(日) 12:32:53例えば
Solve[2 x^2 +x-2==0, x]
とすると、
{{x -> -1.28078}, {x -> 0.780776}}
とでますが、あらかじめ、x は正の範囲に限った答えを出したい(つまり後者の答えだけ出したい)んです。
マニュアル、ヘルプを見ているのですが、よくわかりません。
初歩的な質問で申し訳ありませんがよろしくおねがいします。
0401132人目の素数さん
2009/02/01(日) 13:46:23マジで! どういう経緯でできたんだ
0402132人目の素数さん
2009/02/01(日) 18:23:33日本だけ,ということはないよ。
nCkも,一応海外でも通用する記号。
ただし,国際的に主流の記号ではないので,書物で見る頻度は,列ベクトル記法の方が圧倒的に多い。
かの Donald.E.Knuth の Concrete Mathematics(訳書:コンピュータの数学) には,
「熱心な査読者を除いて,呪わしい記法 nCk は,誰も使いたくない代物である」
と書いてある。
nCk という記号が海外で嫌われる理由は,Cというシンボルだけが目立っていて,
それより大切な n と k が”虐待されて"小さく書かれているから。
ちなみに,C_n^k という記法(Cの右下にn,右上にk)もある。
中国では,この記法が主流だそうだ。
0403132人目の素数さん
2009/02/01(日) 20:23:18(n)
(r)
と表されることはわかったんだけど、
じゃあnPrとかnHrとかは海外ではどうやって表現されてるわけ?
nCrは「n個のものからr個とった組合せ」を表現しているわけだし、
nPrは「n個のものからr個とった順列」を、
nHrは「n個のものからr個とった重複組合せ」を表現してる。
もしたとえばnPrが
(n)
(r)r!
とかで表現されてるんだとしたら、これは本来の意図が伝わりにくくて、nPrの方が優れた記法だと思うんだが。
nHrもしかりで、
(n+r-1)
(r)
としか表現できないんだとしたら(単一の表現がなくて)、nHrの存在ってありがたくね!?
0404132人目の素数さん
2009/02/01(日) 20:25:53「n個のものからr個とった重複順列」
ってシンボルとして存在してたっけ?
なんかあったような気がするんだがど忘れしてしもた
0405132人目の素数さん
2009/02/01(日) 21:41:590406132人目の素数さん
2009/02/01(日) 23:50:59< > や [ ] や { } を使うと重複組合せなり何なりを表すルールだったはず
0407132人目の素数さん
2009/02/02(月) 04:51:41Homogeneous Product だから nHr という記法が使われる。
高校の教科書には掲載されてないけど,「大学への数学」なんかでよく使われる。
>>403
Knuth の本だと,
(n)_r で nPk
(( n r )) で nHk
[ n r ] で第一種スターリング数(n個のものを r 個の円形テーブルに並べる方法の数)
{ n r } で第二種スターリング数(n個のものを r 個の組に分ける方法の数)
を表している(本当は行でなく列記法で)。
>>404
ただの n^r では?
0408132人目の素数さん
2009/02/02(月) 06:20:27とっくに絶滅したと思ってた…
0409132人目の素数さん
2009/02/02(月) 15:30:37マスマティカで多重総和を計算させるとき, 外側の総和の添字の記号って内側の総和の添字、
中の関数にも影響しますかね?
うまく文字で説明できてないと思うんで式を書いておくと、例えば
Σ_{i=0}^t Σ_{k=0}^i F(k)
っていう式があったときに、外側の総和のiが2だったとすると、
内側の総和のiも2になって順次計算されていくか?ということで。
C言語でいうfor文みたいになっているのかなーと。
0410132人目の素数さん
2009/02/02(月) 15:56:220411132人目の素数さん
2009/02/02(月) 18:30:47あれっ、高校の教科書に載ってたと思うけどなぁ、重複組合せ。ちなみに俺25。
>[ n r ] で第一種スターリング数(n個のものを r 個の円形テーブルに並べる方法の数)
多項式の係数としての意味は知ってたんだが、こんな解釈ができたのか。知らなかった。
>ただの n^r では?
いや、そうなんだけど、シンボルがたしかあったなぁと思って。
調べたらΠだった。nΠr
0412132人目の素数さん
2009/02/02(月) 23:43:4825歳だと,複素数平面があった旧課程ですね。
当時の学習指導要領
http://www.nicer.go.jp/guideline/old/h01h/chap2-4.htm
によれば,
>(イ) 組合せ
>[用語・記号]nPr,nCr,階乗,n!
なので,少なくとも高校範囲内ではなかったはず。
ただし,出版社によっては教科書の「発展事項」コーナーなどに掲載されていた可能性も否定できない。
>調べたらΠだった。nΠr
そういえば,そんな記号があったね。思い出しました。
0413132人目の素数さん
2009/02/03(火) 00:27:271年生の時に高校との橋渡しとして習うけど、
すぐに2項係数のベクトル表記を教えられ、「なるべくこっちを使ってください」と言われた。
授業では2項係数をほとんど使わなかった。
0414132人目の素数さん
2009/02/05(木) 09:37:16z[x_, y_] := x^(1/2) y^(1/2);
Plot3D[z[x, y], {x, 0, 10}, {y, 0, 10}, PlotPoints -> 50,
ViewPoint -> {0.70, -3.0, 0.75}]
z=3 で xy平面に平行な平面を描きたいのですが
どうやって描けばよいですか?
どなたかお願いします。
0415132人目の素数さん
2009/02/05(木) 10:25:08とするのです。
0416132人目の素数さん
2009/02/05(木) 10:54:53完璧です。 綺麗にお絵描きできました。
ありがとう。
もうひとつお願いします。
z=3 で xy平面に平行な平面のみをプロットする方法も
お願いします.
0417414
2009/02/05(木) 11:03:02{z[0, y], 3} とするのですね.
0418132人目の素数さん
2009/02/05(木) 11:04:17>>416
でした.
0419132人目の素数さん
2009/02/05(木) 12:17:08z[x, y] = 3 というのは x, y の関数なので
z = 3 は
Plot3D[3, {x, -1, 1}, {y, -1, 1}]
でいいのです。
0420132人目の素数さん
2009/02/05(木) 15:19:20完璧です。
ありがとう.
0421132人目の素数さん
2009/02/05(木) 19:16:23として保持させてからShowする必要があった。
最近のMathematicaはかなり改善されてるが、まだまだ改良の余地はある。
がんばれウルフラム!
0422132人目の素数さん
2009/02/05(木) 20:14:04>>414 or >>416 に関連して
y=5 で xz平面に平行な平面のプロット方法
をどなたかお願いします.
0423132人目の素数さん
2009/02/05(木) 20:17:48Mathematica 持ってないけど。
0424132人目の素数さん
2009/02/05(木) 20:30:03感激です。
素早いレスありがとう。
調べてみたいと思います。
0425132人目の素数さん
2009/02/07(土) 12:54:15なんだっけ?
0426132人目の素数さん
2009/02/07(土) 15:07:380427132人目の素数さん
2009/02/07(土) 19:48:59さっさとバージョンアップしろ
7はすごいぞ
3D回転もデフォルトだ
0428132人目の素数さん
2009/02/07(土) 20:36:030430132人目の素数さん
2009/02/07(土) 21:49:16細部まで綺麗だし
そのうえに
あんなことも
こんなこともできるんだ!
0431132人目の素数さん
2009/02/07(土) 21:54:140432132人目の素数さん
2009/02/07(土) 22:20:30全部わかる。確かに何のソフトかわからんw
0433132人目の素数さん
2009/02/07(土) 23:16:05v7はv5.2よりはメモリを消費するので、v7で計算できないことがv5.2ではできることもある。
グラフィックの質が格段に良くなっているので、計算点数減らしてでもv7を使いたくなるが。
0434132人目の素数さん
2009/02/08(日) 00:37:12マジか
下位ヴァージョンのカーネルも選べるようにしてくれたらいいのに
0435132人目の素数さん
2009/02/08(日) 00:59:49確かに欲しいけど。
0436132人目の素数さん
2009/02/08(日) 17:47:05おまえらどんだけ少ないメモリ積んでんだよ
0437132人目の素数さん
2009/02/08(日) 18:07:110438132人目の素数さん
2009/02/08(日) 20:46:05Manipulateしても
ちゃんと着いてくる
なかなかやるな 7
0439132人目の素数さん
2009/02/09(月) 01:22:183Dのグラフィックスが使いやすくなった。
マウスで3Dの視点をお好みの位置に決めた後、
その位置がどういう Viewpoint なのかを教えてくれる
機能があれば良いのだが..
0440132人目の素数さん
2009/02/09(月) 06:28:190441132人目の素数さん
2009/02/09(月) 09:27:470442132人目の素数さん
2009/02/09(月) 19:32:08windows 7のOSで
Mathematica7これ最強
0443132人目の素数さん
2009/02/09(月) 21:00:57↓はWeatherDataで取得した風速ベクトルのデータだけど、たかが9点を
元にした流線を描かせただけでメモリを1.8GBも消費して(俺の環境では)
落ちる。なんとかならんもんかね。
ListStreamPlot[{{{139.767, 35.683}, {0.740407, -4.74255}},
{{139.833, 35.633}, {12.5036, -7.3423}},
{{139.78, 35.552}, {8.58425, 7.34851}},
{{139.662, 35.837}, {0.349356, 1.35571}},
{{139.528, 35.672}, {4.63184, -1.2594}},
{{139.401, 35.737}, {0.896234, -1.07553}},
{{139.348, 35.748}, {4.86188, 4.16199}},
{{139.384, 35.428}, {-0.593812, -1.26783}},
{{140.386, 35.765}, {-9.48151, -6.14744}}}]
0444132人目の素数さん
2009/02/09(月) 21:18:44>ListStreamPlot[{{{139.767, 35.683}, {0.740407, -4.74255}},
>{{139.833, 35.633}, {12.5036, -7.3423}},
>{{139.78, 35.552}, {8.58425, 7.34851}},
>{{139.662, 35.837}, {0.349356, 1.35571}},
>{{139.528, 35.672}, {4.63184, -1.2594}},
>{{139.401, 35.737}, {0.896234, -1.07553}},
>{{139.348, 35.748}, {4.86188, 4.16199}},
>{{139.384, 35.428}, {-0.593812, -1.26783}},
>{{140.386, 35.765}, {-9.48151, -6.14744}}}]
確かに落ちる
2群に別れてるからかな
0445132人目の素数さん
2009/02/10(火) 13:36:28f0に具体的な数字を入れずに描く方法はありませんか?
f0はそのままグラフ上に反映されて欲しいです.
0446132人目の素数さん
2009/02/10(火) 20:33:300447132人目の素数さん
2009/02/10(火) 21:01:36みたいな定数を表示させたいんじゃない?
0448132人目の素数さん
2009/02/10(火) 21:04:31http://www.wolfram.com/products/mathematicahomeedition/
・価格は295ドル
・機能はプロフェッショナル版と全く同じ
・商用利用・授業での利用等は不可(宿題に使うのは可)
・ダウンロード販売のみ
・現時点では英語版のみ。アメリカとカナダからのみダウンロード可
個人でライセンス持ってる俺の立場は…
0449132人目の素数さん
2009/02/10(火) 22:01:09>(宿題に使うのは可)
宿題に使ったら
バカになるじゃないか
0450132人目の素数さん
2009/02/10(火) 22:27:39アメでは大学院のクラスでもhomeworkが出るよ。日本で言うレポートのこと。
高価なMatlabでも安価な学生版があるわけだし、Mathematicaにあってもおかしくない。
0451132人目の素数さん
2009/02/10(火) 22:43:10学生でもなく、企業で使うのでもなく、家で趣味として数学を楽しむ人
がターゲット。ってどれぐらい需要があるんだろ。
0452132人目の素数さん
2009/02/10(火) 23:50:040453132人目の素数さん
2009/02/11(水) 00:39:54を使用し、その結果を職場のProfessionalで走らせる。
これライセンス違反になるの?
0454132人目の素数さん
2009/02/11(水) 01:00:19これって日本語版もあるの?
0455132人目の素数さん
2009/02/11(水) 01:17:57まだ今月の5日に発表されたばかりで、日本語版は無い。
そのうち出るかもしれんが。
0456132人目の素数さん
2009/02/11(水) 08:58:240457132人目の素数さん
2009/02/11(水) 09:46:54ライセンスがProfessionalとHome分あるなら何も問題無いでしょ。
0458132人目の素数さん
2009/02/11(水) 10:23:59円高メリットが、、、
0459443
2009/02/11(水) 11:31:25偏在してると駄目なんだろうか。
data = Cases[{{WeatherData[#, "Longitude"],
WeatherData[#, "Latitude"]}, WeatherData[#, "WindSpeed"]*
({Cos[#], Sin[#]} &[WeatherData[#, "WindDirection"]])} & /@
WeatherData[{"Tokyo", 50}], {{_Real, _Real}, {_Real, _Real}}];
g1 = ListStreamPlot[data];
g2 = Graphics[{EdgeForm[Black], LightGray,
CountryData["Japan", "FullPolygon"]}];
Show[g2, g1, ImageSize -> 400, AbsoluteOptions[g1, PlotRange],
AspectRatio -> Automatic]
0460445
2009/02/11(水) 16:30:58たとえばy=axのグラフを書こうとするとaに何か値を入れなければなりませんよね.
このとき,aが未定定数のまま,プロットができないかという質問です.
0461132人目の素数さん
2009/02/11(水) 18:37:310462132人目の素数さん
2009/02/11(水) 20:19:35違うと思います。
>>460
a が分からないと右上がりかどうかも分からないので、グラフの形は分かりません。
a に具体的な値をいろいろ当てはめて同じグラフにプロットすることはできます。
0463132人目の素数さん
2009/02/12(木) 00:58:35Manipulate[Plot[Sin[a x], {x, 0, 2 \[Pi]}], {a, 1, 10}]
0464132人目の素数さん
2009/02/12(木) 01:07:53このページのコードを実行してみたんだけど、
切れ目が入ってしまって、画像のような出力は得られない。
どうやらu=πのところっぽいんだけど、なんでだろ?
0465132人目の素数さん
2009/02/12(木) 03:37:460466132人目の素数さん
2009/02/12(木) 06:14:41UnitStepをIf文に直してuの範囲を0〜4Piにしたらそれらしいものに
なったけど、まだ何か違うな
0467132人目の素数さん
2009/02/12(木) 06:20:520468132人目の素数さん
2009/02/12(木) 12:18:13定員割れしてフトコロ寂しいんだろ。察してやれ
0469132人目の素数さん
2009/02/12(木) 16:02:44そうです.
>>462
言われてみれば,そりゃそうですね.
>>463
こりゃいいですね.ありがとう.
0471132人目の素数さん
2009/02/12(木) 19:55:27最初問題なく出力されたけど、
起動しなおしてもう1回やってみたら切れ目が・・・
なんじゃこりゃ
0472132人目の素数さん
2009/02/14(土) 20:55:51でもさ、通常より格段に安いとは言ってもそこそこの値段がするものなんだからパッケージが欲しいなあ
カスタマーサービスに連絡しろって書いてあるけど実際もらえるのかな
0473132人目の素数さん
2009/02/16(月) 11:29:33行列式とか、複数の数式をまとめてカウントする方法って無いでしょうか?
\!\(\*
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{"8", "5", "0", "6", "0"},
{"3", "0", "0", "2", "0"},
{"0", "0", "4", "0", "0"},
{"2", "0", "2", "2", "2"},
{"7", "0", "0", "0", "0"}
}], "\[NoBreak]", ")"}],
Function[ BoxForm`e$,
MatrixForm[ BoxForm`e$]]]\)
こんな感じの行列式の1つの数をカウントする方法をご教授願いたいです
CountだとList要素で、しかも1つのlistしか出来ないみたいなので
0474132人目の素数さん
2009/02/16(月) 11:41:42Count[\!\(\*
TagBox[
RowBox[{"(", "", GridBox[{
{"8", "5", "0", "6", "0"},
{"3", "0", "0", "2", "0"},
{"0", "0", "4", "0", "0"},
{"2", "0", "2", "2", "2"},
{"7", "0", "0", "0", "0"}
}], "", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]\), a_ /; a > 0, 2]
0475132人目の素数さん
2009/02/16(月) 12:02:05例えば
{5,0,0,0,0}という数字の羅列があったとして
0の数字の個数が知りたい場合は
Count[{5,0,0,0,0},0]にすればいいですよね?
それでは、{5,0,0,0,0}、{6,2,0,0,3}、{5,2,0,0,0}などの乱数リストが
100個あった場合、0の数の合計を知りたい時どうすればいいのか判らないのです
0476132人目の素数さん
2009/02/16(月) 13:33:32リストのリストなら、例示したように、Countの第3引数のレベルを2にすればいいんじゃないの?
Count[リストのリスト, a_/;a==0, 2]
で0の数がでるでしょ?
0478132人目の素数さん
2009/02/16(月) 14:54:21解決しました!ありがとうございました!
0479132人目の素数さん
2009/02/18(水) 23:52:580480132人目の素数さん
2009/02/19(木) 07:57:460481132人目の素数さん
2009/02/19(木) 18:33:270482132人目の素数さん
2009/02/19(木) 18:45:480483132人目の素数さん
2009/02/19(木) 19:20:13orz
0484132人目の素数さん
2009/02/22(日) 01:19:15Mathematica 4 をVistaでも使ってる人いる?
0485132人目の素数さん
2009/02/22(日) 04:35:170486476
2009/03/06(金) 16:44:540487132人目の素数さん
2009/03/07(土) 12:23:260488132人目の素数さん
2009/03/07(土) 15:30:110489132人目の素数さん
2009/03/07(土) 22:58:21Mathematicaはスゴイよ。最適化の問題を解こうと思ってMatlabでmpsフォーマットのファイルを
読み込むライブラリか何かネットにないかな?と思って探したけどなかった。しかし、
Mathematicaはデフォルトでサポートw
0490132人目の素数さん
2009/03/08(日) 15:59:24ver6 です
Mathematica のサイトにある
Needs["PlotLegends`"]
ListPlot[Table[{x, f[x]}, {f, {Sin, Cos, Log}}, {x, 0, 10, 0.5}],
PlotLegend -> {"Sine", "Cosine", "Log"},
LegendPosition -> {1.1, -0.4}, Joined -> {True, True, False},
PlotMarkers -> Automatic]
↑すらできません
PlotLegend は Plot, ListPlot で使えるらしいです
http://reference.wolfram.com/mathematica/PlotLegends/ref/PlotLegend.html
0491132人目の素数さん
2009/03/08(日) 23:27:10使う気しないんだけど、なんであんなデザインなんだろう。
0492490
2009/03/08(日) 23:46:15では1つのグラフに複数のデータをプロットする場合は
どのように区別させているのですか?
0493132人目の素数さん
2009/03/09(月) 03:21:43確かに。
俺は Text とか Line とか使って手動で凡例作ってる。
0494132人目の素数さん
2009/03/09(月) 23:32:39ブン・ウルフラム氏が、新たな検索技術「Wolfram|Alpha」を限定公開した。
Wolfram|Alphaは5月に公開予定で、現在は少数のユーザーにのみ公開している。
このプロジェクトは、アルゴリズムや自然言語解析などを活用して、「円周率の300
けた目は?」といった事実に関する質問に回答する「計算型ナレッジエンジン」を目
指している。ウルフラム氏は具体的な仕組みを明らかにしていない。
Wolfram|Alphaのデモに参加した検索エンジンの専門家ノバ・スピバック氏は、この
検索エンジンは「Googleのように答えを含む文書を返すのでもなく、Wikipediaのよう
な巨大なデータベースでもなく、Powersetのように単に自然言語を解析して文書を
検索するのでもない」としている。
Wolfram|Alphaは、「ティンブクトゥはどこの国にあるか」「水素原子の陽子は幾つか」
「今月のシアトルの平均降水量は?」「ISSはどこか」「GOOG(Google)の株価が300ド
ル以上だったのはいつか」といった事実に関する質問に回答してくれるという。
スピバック氏は、Wolfram|Alphaは幅広い質問への答えを提供し、データベースから
答えを探すのではなく回答を計算するという点でGoogle よりもずっとスマートだとし
ている。「WebにとってGoogleと同じくらい重要かもしれない」「Googleの人がこれを見
たら驚いて、欲しくなるだろう」と賞賛している。
http://headlines.yahoo.co.jp/hl?a=20090309-00000039-zdn_n-sci
0495132人目の素数さん
2009/03/10(火) 00:16:160496132人目の素数さん
2009/03/10(火) 18:38:370497132人目の素数さん
2009/03/10(火) 19:44:500498132人目の素数さん
2009/03/10(火) 21:20:53物理に寄与した?
0499132人目の素数さん
2009/03/10(火) 21:45:400500132人目の素数さん
2009/03/11(水) 02:11:00俺は入教してしまうだろう
0501132人目の素数さん
2009/03/11(水) 02:30:010502KingGold ◆3waIkAJWrg
2009/03/11(水) 04:30:170503132人目の素数さん
2009/03/12(木) 00:12:310504132人目の素数さん
2009/03/21(土) 19:32:250505132人目の素数さん
2009/03/22(日) 08:56:35以上
0506132人目の素数さん
2009/03/22(日) 10:19:110507132人目の素数さん
2009/03/22(日) 13:01:40弟も役員でウハウハ
一家安泰
0508132人目の素数さん
2009/03/25(水) 20:24:350509132人目の素数さん
2009/03/29(日) 23:39:440510132人目の素数さん
2009/04/03(金) 11:26:47どっかにチュートリアルとか使い方の解説はないかな。
ヘルプ見ても個別のコマンドの短い説明が載ってるだけで、ノートブックでの統合的な使い方が
わからない。Webで探したけど、不思議なほど見つからないんだよね。
0511132人目の素数さん
2009/04/03(金) 21:43:530512132人目の素数さん
2009/04/05(日) 22:35:40Mathematica Navigatorという本に入門程度は書いてあった。
日本のMathematica本はみんな古くて、ver6以上扱ってるのないな。
0513132人目の素数さん
2009/04/06(月) 18:40:34もう本は出ないか
0514132人目の素数さん
2009/04/06(月) 23:15:07一から勉強したいとか、面白い機能がないかパラパラ探すにはまだ
本のほうがいいな。その点」Wolframの本はよかったんだが。
0515132人目の素数さん
2009/04/07(火) 20:16:07ソフトに入ってるんじゃなかったっけ
0516132人目の素数さん
2009/04/09(木) 10:08:30トピックの寄せ集めって感じ。実質的な中身は同じなんだろうけど、本に比べると読みづらいんだな。
0517132人目の素数さん
2009/04/09(木) 11:47:400518132人目の素数さん
2009/04/11(土) 14:19:290519132人目の素数さん
2009/04/11(土) 15:13:20論文一本でも多く出しておくことのほうが有利に評価される。
0520132人目の素数さん
2009/04/11(土) 15:27:23全然有利にならない。RDBMSのSQL使える方がよっぽどウリになる。
0521132人目の素数さん
2009/04/11(土) 15:50:250522132人目の素数さん
2009/04/11(土) 19:27:37あるから少しは有利になるかもしれんけど、操作法なんて
すぐ覚えられるからなあ。むしろ使いこなすだけの数学的
知識が必要。
0523132人目の素数さん
2009/04/11(土) 21:53:55調べて見ても同様のケースについては全然報告がないので当惑している。
今度の学会に、今から他言語で一から組み直すにはもう時間がない・・・・
かといって本学のVer6.0の契約は今月いっぱいまでだろう、・・・・・正直困った。
具体的には
・フロントエンド、カーネルが起動するが1+1すら演算中のまま答えを返さない、ずっとそのまま。
・ドキュメントセンターを開こうとする、あるいはネットワークのテストを行おうとするとうまくいかない。
・関数自体の動作は起動のたびにjavaw.exeというプロセスを止めると回復する、
これでもドキュメントセンター、ネットワーク系ははおかしいまま。
・Administratorでインストール、起動した場合はすべてうまく動作するようである。
・カーネルだけの動作ならどの権限でもうまく動作する。
何か指針がありましたら、ご教示ください・・・・。よろすく・・・・。
0524KingGold ◆3waIkAJWrg
2009/04/12(日) 07:29:340525132人目の素数さん
2009/04/12(日) 15:22:38インストールして一回目の起動ではそういう現象も有った。
>1+1でずっと演算中…
pc再起動したらなくなったよ
ちなみにWindowsXP 64bit
確かver6は平行して使えると思ったけど??
0526132人目の素数さん
2009/04/12(日) 21:27:50>pc再起動したらなくなったよ
インストール以来再起動することは何度もあるんですが
やはりうまくいかないままです。
0527132人目の素数さん
2009/04/24(金) 00:26:34向こうも理由が分からないならVer6.0をもう少し使わせてもらえないか掛け合ってみるとか?
0528523
2009/04/27(月) 03:58:05システム課の方と戦って、なんとか動くようになりました
プロキシの設定をうまく参照できていないようですねぇ・・・。
0529132人目の素数さん
2009/04/29(水) 19:02:03起動のたびにいちいち確認するのな
シングルアロンならそんな問題はない
0530132人目の素数さん
2009/04/29(水) 20:48:510531132人目の素数さん
2009/05/06(水) 17:18:55Ordering2D[mat_] :=
Transpose@{Quotient@##, 1 + Mod@##} &[
Ordering@Flatten@mat + # - 1, #] &[Dimensions[mat][[2]]]
(* テスト *)
MatrixForm /@ (test =
Flatten[Table[RandomInteger[10, {row, col}], {row, 10}, {col, 10}],
1])
Table[
{# == Sort@#} &[test[[n, Sequence @@ #]] & /@ Ordering2D[test[[n]]]],
{n, Length@test}]
変則的なやつ{{1,2},{3}}とか、3次元以上{{1,2},{{3,{4}}},5}
にも使えるOrdering関数を寛大でエロい人下さい!
0532132人目の素数さん
2009/05/06(水) 23:59:30Flatten[Position[list, #] & /@ Sort@Flatten@list, 1]
0534132人目の素数さん
2009/05/08(金) 22:31:210535132人目の素数さん
2009/05/09(土) 15:39:07自分のホームページで宣伝するしかないのかな?
0536132人目の素数さん
2009/05/09(土) 19:11:21* 中核となる画像処理関数の性能アップ
* 基本数学,数学授業,文章作成のアシスタントパレットの日本語化
* 素早い画像処理を可能にする右クリックメニュー
* 新しいチュートリアル,「How to」ガイド,スクリーンキャスト
* ドキュメントに含まれる数千に及ぶ新規例題
* 向上したドキュメント検索
* Windows 7の手書き数式認識機能との統合
* 新しいgridMathematica Serverとの統合
0537132人目の素数さん
2009/05/26(火) 23:25:43下記HPに記載の計算をwindows PCで評価するのですが、
メモリーが足りなくて途中で終了してしまいます。環境は
次の通りですが、最後まで計算させる方法を教えてください。
また計算結果(グラフ)がどのような結果なのか?教えて下さい。
なお、計算時間は5時間以上必要です。
計算プログラム
http://www.geocities.jp/etctransformation/mathematica1-7.html
PC環境
Dell
Inspiron 530
Windows Vista Home Basic
メモリ 100GB
プロセッサ Intel Core 2 Duo CPU 2.2GHz
0538132人目の素数さん
2009/05/27(水) 08:21:50>>メモリ 100GB
。 。
/ / ポーン!
( Д )
0539132人目の素数さん
2009/05/27(水) 18:58:48Windows Vista などというクソは
100GBのメモリを乗せても
3GBとしか使えないのだ
出直してこい!
0540猫でつ ◆ghclfYsc82
2009/05/27(水) 19:41:36クソばかりだが、リナックスに入れ替えようとするとデバイスドライバが無い
から困りますよね。この「糞みたいな状況」は何とかなりませんかね。
猫の希望はソラリスですが、もしダメならデビアンでもエエのですが。
誰かどうしたらエエか教えて下さい。
(恐らくクノーピクスで解決なのではと思っていますが。)
0541132人目の素数さん
2009/05/27(水) 21:33:480542猫でつ ◆ghclfYsc82
2009/05/27(水) 21:43:540543132人目の素数さん
2009/05/28(木) 00:17:570544132人目の素数さん
2009/05/28(木) 00:21:060545猫でつ ◆ghclfYsc82
2009/05/28(木) 00:47:30今だったらプローブとかもあるんだろうけど。
とにかくノートって安く小さくする為に変な部品を使ってますよ。
0546132人目の素数さん
2009/05/28(木) 00:51:29あんまり苦労せんがのお・・・
0547猫でつ ◆ghclfYsc82
2009/05/28(木) 01:27:31そやけど本当はソラリスが一番好きなんですが。
0548132人目の素数さん
2009/05/28(木) 21:50:56http://www.amazon.co.jp/dp/486190465X
Linux Perfect Vol.3
http://www.amazon.co.jp/dp/4048679244/
週刊アスキー別冊 さくさくUbuntu! 新バージョン9.04 CD-ROM付き
0549猫でつ ◆ghclfYsc82
2009/05/28(木) 22:10:32そんでソラリスだけは今でもかなり「相性」があるみたいで、
ちょっとな〜 って感じなんですが。
でもドライバが揃い過ぎて複雑怪奇になったヤツよりかは
マシなんですかねw (一番酷いのがマイクロソフトでしょ。)
0550やんやん ◆yanyan72E.
2009/05/29(金) 07:40:14このレスを見てから、2日間ぐらいMacProをぶん回して
計算しているんだけどいまだに計算おわんないよ?
まぁ計算は進んでいるみたいだから
2週間ぐらいまでなら最後まで計算させてみるけれど。
ちなみにMacProのスペックは
Core2Duo Xeon 2.66MHz x2
メモリ6GB
0551132人目の素数さん
2009/05/29(金) 17:10:09>Core2Duo Xeon 2.66MHz x2
ばーか
そんな石ねーよ
0552やんやん ◆yanyan72E.
2009/05/29(金) 17:27:20Core2アーキテクチャのデュアルコアXeonって意味ね。
具体的にはWoodcrest。
0553猫でつ ◆ghclfYsc82
2009/05/30(土) 11:18:03数値計算ですか?
0554132人目の素数さん
2009/05/30(土) 20:04:28Cっぽく書くと速くなる位しか分かりません
0555132人目の素数さん
2009/05/30(土) 23:23:06HPはサービス内容の宣伝はいいだけど、肝心の価格が見当たりませんでした。
0556132人目の素数さん
2009/05/30(土) 23:35:48カビのはえたような時代遅れのCっぽく書くいたって遅くなるだけだろwww
4GLライクに記述できる部分は4GLライクに記述汁。
0557132人目の素数さん
2009/06/01(月) 14:39:56コード晒せば
0558132人目の素数さん
2009/06/02(火) 21:27:19長時間の計算、有難うございます。申し訳ございません。一部コードが間違っていました。修正しましたので再度計算頂きましたら幸いです。
http://www.geocities.jp/etctransformation/mathematica1-7.html
但し、5日間くらい必要だと思われます。表示した図を教えてください。
0559132人目の素数さん
2009/06/04(木) 10:06:07どうしたらよいでしょうか?
フラッシュというのでできますでしょうか?
それともjavascriptでできますでしょうか?
0560132人目の素数さん
2009/06/04(木) 18:46:35なんでもできるぞ
かんがれ
0561猫でつ ◆ghclfYsc82
2009/06/04(木) 18:56:18そんでジャバやったら関数でもライブラリでもいっぱいあるでしょうが
0562132人目の素数さん
2009/06/04(木) 21:50:14サーバーサイドでしょうか?
クライアントサイドでしょうか?
マスマティカをhtmlに埋め込むか、サーバーに置くか
という選択肢もありでしょうか?
0563変態ロリーマン ◆ghclfYsc82
2009/06/04(木) 22:05:22した方がベターじゃないんですか
だから使用するポートは80番その他だけにして、
後はセッションを張った時だけに
ローカルのCPUパワーを喰う様にした方が
良いのではないでしょうか
0564132人目の素数さん
2009/06/13(土) 15:27:25ある陰関数をContourplotでPlotして得られる曲線があるのですが
曲線上のある点(x,y)を境に2つに分けて、片方を点線、片方を実線で
プロットしてやりたいのですが、どうすればいいでしょう。
0565132人目の素数さん
2009/06/15(月) 20:05:110566132人目の素数さん
2009/06/15(月) 22:14:030567132人目の素数さん
2009/06/15(月) 22:51:05微分方程式の解き方ですが
yl''[x]==-k x (0<=x<=L/2)
yr''[x]==-k(L-x) (L/2<=x<=L)
B.C. y[0]=0,y[L]=0,yl[L/2]=yr[L/2],yl'[L/2]=yr[L/2]
これをDSolveを使って解く方法を教えて下さい。
あるいは、DSolveを使わない方法でもかまいません。
どうか、お願いします。
0568132人目の素数さん
2009/06/16(火) 06:50:26どうなるでつか?
0569132人目の素数さん
2009/06/16(火) 07:08:340570132人目の素数さん
2009/06/16(火) 08:33:19どっかに日本語ヘルプが読めるv7体験版ありませんか?
ベクトル場グラフィックをがどの程度なのかきになるんです。
0571132人目の素数さん
2009/06/16(火) 20:54:500572132人目の素数さん
2009/06/17(水) 09:09:09v6と違ってちゃんと日本語だったんた。ありがとう。
使用期間もあとちょっとだし、てきとーにいじって遊んでみる。
0573猫ぱんだ ◆ghclfYsc82
2009/06/19(金) 08:18:38ウォルフラムのヲッサンはお商売が巧いなぁ
0574132人目の素数さん
2009/06/19(金) 11:33:040575132人目の素数さん
2009/06/19(金) 20:41:160576132人目の素数さん
2009/06/22(月) 10:10:35実験データを片対数グラフにプロットし,
ListLogLinearPlot[・・・・・,Joined->True] でデータ同士を線で結んだんだけど,データとデータの間が直線で結ばれてしまう('・ω・`).
Excelの散布図の様に,なめらかな曲線でデータ同士を結ぶにはどうしたらイイ?
0577132人目の素数さん
2009/06/26(金) 14:59:16{r、0、10}、{φ、-π、π}の範囲でグラフィックにしたいんですけど、
VectorPlot[]
だと、デカルト座標(x、y)として読みこまれてしまって、
ほしい図がでてこないんですが、どうしたらいいですかね。
0578132人目の素数さん
2009/06/26(金) 15:21:220579132人目の素数さん
2009/06/27(土) 12:46:370580132人目の素数さん
2009/06/28(日) 13:07:18これでどうでしょう。
rotate[x_, y_, t_] := {x Cos[t] - y Sin[t], x Sin[t] + y Cos[t]}
Ar[r_, t_] := 0
At[r_, t_] := 1;
With[{r = Sqrt[x^2 + y^2], t = ArcTan[x, y]},
VectorPlot[rotate[Ar[r, t], At[r, t], t], {x, -10, 10}, {y, -10, 10}]]
0581132人目の素数さん
2009/06/29(月) 15:25:04PlorPlotだとベクトルが描けないような気がします。
>>580
できました!ありがとう!
当初、極座標成分ベクトルをデカルト座標成分ベクトルに直してプロットするところまでは
思いついたんですが、逆三角関数の処理がうまくできなくて困ってました。
自分の書いたプログラムの方も
ArcTan[y/x]→ArcTan[x,y]
に直したら正しい図が出てきましたw
0582132人目の素数さん
2009/06/30(火) 23:12:03もう解決してしまったかもしれないが、一応。
DSolve[{yl''[x] == -k x, yr''[x] == -k (L - x), yl[0] == 0,
yr[L] == 0, yl[L/2] == yr[L/2], yl'[L/2] == yr'[L/2]}, {yl[x],
yr[x]}, x]
Block[{L = 1, k = 1}, Plot[If[x <= L/2, yl[x], yr[x]] /. %, {x, 0, L}]]
0583132人目の素数さん
2009/07/01(水) 10:07:03やりたいことはこういうことでしょうか?
data = {{1, 1.0}, {10, 2.2}, {100, 3}, {1000, 4.4}};
Fit[data, {1, Log[x]}, x]
g1 = LogLinearPlot[%, {x, MapThread[Min, data][[1]]/10,
MapThread[Max, data][[1]]*10}];
g2 = ListLogLinearPlot[data];
Show[g1, g2]
0584132人目の素数さん
2009/07/05(日) 03:14:500585132人目の素数さん
2009/07/16(木) 23:12:55y^2/b^2+z^2/c^2-x/a=0
この関数をxyz平面にプロットしたいんだけどどうすればいいんでしょう
0586132人目の素数さん
2009/07/16(木) 23:13:49ごめんなさい
0587132人目の素数さん
2009/07/17(金) 18:05:46xyz平面とは何か
0588132人目の素数さん
2009/07/17(金) 18:25:18まあまあ、いいじゃないですか。
>>585
バージョン6以降が使えるのだったら、
Manipulate[
ContourPlot3D[x^2/a^2 + y^2/b^2 - z/c == 0, {x, -2, 2}, {y, -2, 2}, {z, -10, 10}]
, {a, 1, 4}, {b, 1, 4}, {c, -1, 1}]
などとすればよいです。
0589132人目の素数さん
2009/07/21(火) 21:47:08こんにちは、
y=t^4 Exp[-t^2];をフーリエ変換すると、
y2=FourierTransform[t^4 Exp[-t^2],t,w];
Print[y2];
になります。
この計算を、まともに下記で確認したいのですが、答えが一致しません。
下記のどこを修正すれば良いのでしょうか?
y=.;
y=t^4 Exp[-t^2];
y1=Integrate[y*Exp[-I*w*t],{t,-Infinity,Infinity}];
Print[y1];
下記HPのNo2の計算を、MATHEMATICAで確認したいのですが、プログラムを教えてください。
http://www.geocities.jp/dirac_equation/img010.pdf
0590132人目の素数さん
2009/07/21(火) 23:22:31うーん、マジレスしないといけないんだろうか?
y1 = FourierTransform[t^4 Exp[-t^2], t, w];
y2 = 1/Sqrt[2 Pi] Integrate[t^4 Exp[-t^2] Exp[-I w t], {t, -Infinity, Infinity}];
y1 - y2
とりあえずここで止めとく。反応があれば考えるかも。
0591132人目の素数さん
2009/07/22(水) 07:20:500592132人目の素数さん
2009/07/22(水) 08:25:27>>590
お返事有難う御座います。
y2 = 1/Sqrt[2 Pi] Integrate[t^4 Exp[-t^2] Exp[-I w t], {t, -Infinity, Infinity}];
の積分計算がうまくできません。if==何とかが出ます。設定が悪いのでしょうか?
0593132人目の素数さん
2009/07/22(水) 18:44:06if==何とかじゃわからん。エラーメッセージを全部書きなさい。Mathematica のバージョンも書くこと。
もし初心者だったら、掛け算にはスペースを使わず、* を使ったほうがいい。
0594132人目の素数さん
2009/07/23(木) 22:52:23下記の通りです。Ver4.2です。
y2=1/Sqrt[2 Pi] Integrate[t^4 Exp[-t^2] Exp[-I w t],{t,-Infinity,Infinity}]
General::stop: 計算中,Series::esssのこれ以上の出力は表示されません.
SeriesData::argbu: SeriesDataは1つの引数で呼ばれました.3個と6個の間の引数が想定されています.
If[Im[w] == 0, ((1/16)*Sqrt[Pi]*(12 - 12*w^2 + w^4))/E^(w^2/4), Integrate[E^(-t^2 - I*t*w)*t^4, {t, -Infinity, Infinity}]]/Sqrt[2*Pi]
0595132人目の素数さん
2009/07/24(金) 06:17:40問題なく答えが出てるじゃないか。w の虚部が零か非零かで場合わけされてるんでしょ。
0596132人目の素数さん
2009/07/29(水) 11:16:54Plot[x, {x, -1, 1}]
といった具合にPlotでグラフ描画すると全部シャープ記号(#)で描画されて
めちゃめちゃ荒いアスキーアートになってしまいます。
普通に細い線でグラフを描くにはどうすればいいのでしょうか
0597132人目の素数さん
2009/07/29(水) 12:43:22実行環境を詳しく
0598132人目の素数さん
2009/07/29(水) 22:03:41kernel上でやってたのが原因でした・・
お騒がせしてすみませんでした;;
実行環境はWindowsXP,Mathematica5.0 kernelでした
0599132人目の素数さん
2009/08/01(土) 16:35:09教えて下さい。
普通にコピーしてテキストファイルに貼り付けて、いらない部分をデリートした後、
Excelでスペース区切りで開くと確かにデータは読み込めて、グラフにもすることが出来るのですが、
縦軸と横軸に何か縮尺がかかっており、また縦軸にオフセットが乗っています。
どうすれば正常に読み込ませることが出来るのでしょうか?
0600「猫」∈社会の屑 ◆ghclfYsc82
2009/08/01(土) 17:55:16数学は出来しませんねん!
それよりも自分の頭を使って何か考えたらどうですねん!
0601132人目の素数さん
2009/08/01(土) 18:10:09てめー、関係ない所まで首突っ込んで
何も悪くない人間まで罵倒するようになったか。
Mathematicaが仕事で必要な人間もいるんだよ。
自分と他人、他人Aと他人Bが、それぞれ違う立場の
人間であることを理解できないなら人間やめろ。
あ、もうやめてるか。
0602「猫」∈社会の屑 ◆ghclfYsc82
2009/08/01(土) 18:16:36意見したらアカンのか!
0603「猫」∈社会の屑 ◆ghclfYsc82
2009/08/01(土) 18:29:43自分と他人のそれぞれ違う立場が理解出来ひんんは
どっちやねん!
そういう事を何時もしとんはオマエとちゃうんか!
しかも此処は誰でもカキコ出来る自由掲示板やないけ!
オマエ、口に気を付けへんとエラい事になるゾ
判ったか!
0604そら
2009/08/02(日) 00:42:53>そんな「ソフトの使い方」なんか覚えても数学は出来しませんねん!
それは使う人に依ると思う。あなたが必要ないなら覚えなければ良い話。
0605132人目の素数さん
2009/08/02(日) 12:05:13こんなんじゃだめですか?
Table[{x, Sin[x]}, {x, 0, 2 Pi, 0.1}]
Export["graph.xls", %]
Ver 5.0 が xls 形式をサポートしていないようだったら、
Table[{x, Sin[x]}, {x, 0, 2 Pi, 0.1}]
Export["graph.csv", %]
とすればいいのかも。(Ver 5.0 が手元にないのでわからない。)
0606132人目の素数さん
2009/08/03(月) 20:47:42mathematicaでのカイ二乗統計量
x二乗検定のプログラムをmathematicaで作成せよ。
(1) データ2項目のカテゴリカルデータ
1: 3 2
2: 2 1
3: 1 1
(2) クロス集計表を作成せよ。
(3) 独立性検定に用いるx二乗統計量を計算せよ。
(4) p値を求めよ。(x二乗分布の分布関数から)
(3)は、d1 = RandomChoice[{0.8, 0.2} -> {1, 2}, 100];(*A大 1:男、2:女*)
d2 = RandomChoice[{0.4, 0.6} -> {1, 2}, 100];(*1:大学院進学 2:その他*)
d0 = {d1, d2};
d = Transpose[d0]と
入れた後に、
MeanTest[d, 1, FullReport -> True]
としました。
この結果、
{FullReport ->
TableForm[{{{1.15, 1.67}, {4.1797832761154154`, 14.177446878757825`}, StudentTDistribution[99]}}, TableHeadings -> {None, {"Mean", "TestStat", "Distribution"}}], OneSidedPValue -> CDF[StudentTDistribution[99], {-4.17978, -14.1774}]}
と出ました。
検定統計量は、mean1.15に対し、teststatが4.17978、mean1.67に対し、teststatが14.1774となりました。
検定統計量は2つでることは間違いではないでしょうか。
早急にお返事ください。
また、(4)はNeeds["HypothesisTesting`"]
ChiSquarePValue[0.3, 1]でいいんでしょうか?
0607132人目の素数さん
2009/08/04(火) 23:26:53mathematica4.1,windows2000を使用しております。
BaseForm[200,2]を評価すると、そのままBaseForm[200,2]と表示されて答えが出ません。どうしたら計算するでしょうか?
0608132人目の素数さん
2009/08/05(水) 10:12:481+1 などの他の計算は実行されるでしょうか?
0609132人目の素数さん
2009/08/05(水) 19:48:06>1+1 などの他の計算は実行されるでしょうか?
大丈夫です。
0610132人目の素数さん
2009/08/06(木) 17:06:33ってことはカーネルは動いているってことなのかな?
情報が少なすぎてよくわからん。スペルミスとかではないですよね?
0611132人目の素数さん
2009/08/20(木) 17:57:09Clear[dt]
Hold[dt[7]] /. dt[hikisuu_] :> hikisuu
がHold[7]になりますが、なぜ
Hold[Dt[7]] /. Dt[hikisuu_] :> hikisuu
でHold[7]とならず、Hold[Dt[7]]となってしまうのでしょうか?
Hold[Dt[7]] /. _?(Dt == # &)[hikisuu_] :> hikisuu
で無事Hold[7]になりましたが、上の原因が知りたいです。
(Windows VistaでMathematica 7.0.0です)
0612132人目の素数さん
2009/08/20(木) 21:42:23全微分ってのは、中で使われてる変数の名前に応じてその値が異なるから、
Dt[x]:>dx みたいにして使うことはあるけど、
例:http://reference.wolfram.com/mathematica/tutorial/TotalDerivatives.ja.html の最後
Dt[x_]で定義するのがそもそもおかしいな。
0613611
2009/08/20(木) 22:42:24レスありがとうございます。例が不適切&説明不足でした。
Dtは定義済みの関数の一つという意味で使いました。
パターンがマッチしてるのになんで置き換えてくれないんだ〜と
0614132人目の素数さん
2009/08/20(木) 23:46:23定義済み関数の中でも、Dtは特殊なんだ。変数の名前に依存するから。
Sin とかでやってみれば違いがわかる。
0615611
2009/08/23(日) 22:32:10すいません確認が遅れてしまいました。レスありがとうございます。
Sinで試したらおっしゃる通りうまくいきました。定義されてるかどうかで決まるわけではなかったんですね。
Hold[Dt[7]] /. Hold[Dt[hikisuu_]] :> hikisuu
にするとどの関数でも変換するようなので、とりあえずHold付けておきます。
重ねて、レス下さった方々有り難うございました。
0616132人目の素数さん
2009/09/02(水) 13:18:32xが0のときはどういう解釈なんでしょう?
0617132人目の素数さん
2009/09/04(金) 17:23:311/0はComplexInfinityと一応定義されてるからじゃ?とおもったけど
単純に1/x == 1/xがTrueだからみたいね
ForAll[x, x == 0, 1/x == 1/x] // Trace
0618132人目の素数さん
2009/09/04(金) 18:51:09そんで方向が分からない時はイコールは未定義と。
全然知らんかった。
0619132人目の素数さん
2009/09/07(月) 21:33:20卒業後もそれを使えるの?
0620132人目の素数さん
2009/09/14(月) 01:35:38どう書けばいいのでしょうか。
0621132人目の素数さん
2009/09/17(木) 12:00:360622132人目の素数さん
2009/10/04(日) 03:20:15ってわかる人いない?
0623132人目の素数さん
2009/10/06(火) 22:54:440624132人目の素数さん
2009/10/08(木) 16:52:50f[z_] := 1/Log[c] Log[Abs[(z - c)/(c z - 1)]] /. {c -> (3 - Sqrt[5])/2}
ContourPlot[Abs[f[x + y I]], {x, -1.2, 1.2}, {y, -1.2, 1.2},
PlotRange -> {0, 1}, PlotPoints -> 50, Mesh -> False,
Epilog -> {{PointSize[0.05], Hue[0.5],
Circle[{0, 0}, 1]}, {PointSize[0.05], Hue[0],
Circle[{1/3, 0}, (3 - Sqrt[5])/2]}
}]
0625132人目の素数さん
2009/10/13(火) 18:04:500626625
2009/10/13(火) 18:07:230627132人目の素数さん
2009/10/17(土) 08:55:32z=g(x,y)
の2曲面の交わる線がどんな形してるかを調べたいのですが
やり方が思いつきません。
z=f(x,y)
y=2x+3
のように曲面と直線の交わる線ならやり方はわかります。
x=t
y=2t+3
z=f(t,2t+3)
とやってParametricPlot3Dすればいいと思うのですが、
曲面と曲面の交わる線の表示がわかりません。
よろしくお願いします。
0628132人目の素数さん
2009/10/17(土) 22:33:13ParametricPlot3D[{x,y,z}/.Solve[{z==f(x,y),z==g(x,y)},{x,y}],{z,-1,1}]
とか
0629132人目の素数さん
2009/10/17(土) 22:52:06ぐぐると去年の11月ごろの記事で、「Mathematicaは来年初頭に対応する」というのが
ぞろっと出るけど、それ以降のがでてこない。wolframのサイトでも見つからないし。
うまくいかなかったって感じなのかな。
0630132人目の素数さん
2009/10/18(日) 05:03:01とやってみましたができませんでした。。
0631132人目の素数さん
2009/10/18(日) 09:53:36???
0632132人目の素数さん
2009/10/18(日) 13:49:54お願いします
ver4です。
0633132人目の素数さん
2009/10/21(水) 12:35:29LogPlot[Pe[-Sqrt[R]], {R, 5, 15}]
がほとんど直線になります
QPSKという変調方式の誤り率なんですが
教科書では曲がっておりました
なぜでしょう???
0634132人目の素数さん
2009/10/21(水) 15:43:420635132人目の素数さん
2009/10/21(水) 17:37:14何をしたいのかさっぱりわからん
ソルウ゛にz要らんやん
0636132人目の素数さん
2009/10/24(土) 00:08:49どういう形をしているのかを知りたいというわけです。
よろしくお願いいたします
0637132人目の素数さん
2009/10/25(日) 19:32:02ちょっとズルですが、以下のようなやり方はどうでしょう?
f[x_, y_] := Sin[x + y]
g[x_, y_] := x y
h[{x_, y_}] := {x, y, g[x, y]};
points = ContourPlot[f[x, y] == g[x, y], {x, -10, 10}, {y, -10, 10}][[1, 1]];
ListPointPlot3D[Map[h, points]]
0638132人目の素数さん
2009/10/26(月) 07:00:59しかしできませんでした。。
0639132人目の素数さん
2009/10/26(月) 13:07:48ContourPlot[Sin[x + y] == x y, {x, -3, 3}, {y, -3, 3}]
では?
0640132人目の素数さん
2009/10/26(月) 18:53:54mathematicaで出力をベクトルで表示することって可能ですか?
a={ax,ay,az}
b={bx,by,bz}
c={bx,by,bz}
d={bx,by,bz}
って入力したとして、
Dot[a,Cross[b,c]]+Dot[a,Cross[c,d]]
とかの出力は
az(-by cx +bx cy)+・・・
って感じでベクトルの成分で出力されてしまうんですが
a・(b×c+c×d)
みたいにベクトルのままで出力されてほしいんです。
これってできますか?
0641132人目の素数さん
2009/10/28(水) 19:15:05最後に、//TraditionalForm
0642132人目の素数さん
2009/11/07(土) 00:18:31ベクトルとマトリックスをリストによって表現するのはwolframのミスではなかったか。
ベクトルの集合を作って、この集合に何か作用させようとするときに、リストを使うと簡単に行えない。
0643132人目の素数さん
2009/11/07(土) 00:35:33どのようにしたらフィッティングできるでしょうか?
FindFitとNDSolveを組み合わせるのでしょうか?
ご教授よろしくお願いいたします。
ちなみにver7を使っています。
0644132人目の素数さん
2009/11/07(土) 22:48:34明示的に列行列や行行列にするのはだめ?
{{1}, {0}} とか{{1,0}}とか。
0645132人目の素数さん
2009/11/08(日) 15:00:02ドキュメントセンターのFindFitの項目に例題がいっぱい載っている。
例えば
data = {{1.0, 12.}, {1.9, 10.}, {2.6, 8.2}, {3.4, 6.9}, {5.0, 5.9}};
model[a_?NumberQ, k_?NumberQ] :=
model[a, k] =
First[x /. NDSolve[{x'[t] + k x[t] == 0, x[0] == a}, x, {t, 0, 6}]]
fit = FindFit[data, model[a, k][t], {a, k}, t]
Show[Plot[model[a, k][t] /. fit, {t, 0, 6}], ListPlot[data]]
のようにすればよい。
0646132人目の素数さん
2009/11/08(日) 15:01:06間違えた。643への返信でした。
0647132人目の素数さん
2009/11/08(日) 16:15:30あまりの優しさに感動しました。
本当にありがとうございました。
またなにかございましたらよろしくお願いいたします。
0648132人目の素数さん
2009/11/09(月) 04:04:330649132人目の素数さん
2009/11/09(月) 04:10:450650132人目の素数さん
2009/11/09(月) 04:19:460651132人目の素数さん
2009/11/09(月) 04:36:26D[x a^y,x]
0652132人目の素数さん
2009/11/09(月) 04:52:130653132人目の素数さん
2009/11/09(月) 05:09:21Series[(y^2)(a^x)+(x)(a^(y^2)),{x,0,3},{y,0,3}]
0654132人目の素数さん
2009/11/09(月) 10:15:270655132人目の素数さん
2009/11/09(月) 14:54:140656132人目の素数さん
2009/11/10(火) 20:02:350657132人目の素数さん
2009/11/10(火) 20:21:25自分でもっと探してみな。
0658132人目の素数さん
2009/11/19(木) 14:15:210659132人目の素数さん
2009/11/19(木) 17:07:59たとえば、
(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o)
とならんでいるものを、
(a,b,c)(d,e,f)(g,h,i)(j,k,l)(m,n,o)
みたいに3つずつのくくりにしたいです。
お願いします。
0660132人目の素数さん
2009/11/19(木) 17:14:530661132人目の素数さん
2009/11/22(日) 13:00:12せっかくの円高なのに、、、。
0662132人目の素数さん
2009/11/24(火) 11:17:36カルキングはクソ。買って損した
0663132人目の素数さん
2009/11/24(火) 12:55:160664132人目の素数さん
2009/11/24(火) 13:09:56機能はプロフェッショナルと同じらしい。
しかし値段は無限と有限の違いだよ。
アメリカでは2月から売られてるようだが、日本では気配の気の字も感じられないな。
http://www.wolfram.com/products/mathematicahomeedition/qa.html
Q: How is Mathematica Home Edition different from the professional version of Mathematica?
Mathematica Home Edition is a fully functional version of Mathematica Professional with the same features.
0665132人目の素数さん
2009/11/27(金) 18:17:59図を見ている場所を知る方法ってありますか?
ViewPointで設定できることはわかるのですが、逆に見ている状態がどこにいるのかが知りたいです。
比較検討のために必要なので、わかる人がいましたらお願いします。
0666132人目の素数さん
2009/11/27(金) 20:38:38マジで!?
これって知人に海外旅行のついでに買ってきてもらうとかいうわけにはいかないの?
0667132人目の素数さん
2009/11/28(土) 09:39:14どうなんだろうね。
しかしダウンロード販売してるのに、わざわざそんなことしなきゃいけないのかね。
日本の代理店とかに聞くと、ウルフラム社が日本とかに売らないようにしてる、とか言うが、
日本の代理店の方が高い方しか売らないでくれとお願いしてるんじゃないかと疑っている、、。
これは公正取引法違反じゃないのか?
0668132人目の素数さん
2009/11/28(土) 15:45:51Home Edition すらいらないじゃんってなるかも。
http://www.hirax.net/keywords/log/IronRuby/latest
これと、そこから貼られているリンク先の記事参照。
0669132人目の素数さん
2009/11/28(土) 17:35:21そのうち手が打たれそうな気もする。
0670132人目の素数さん
2009/12/06(日) 00:07:09何かを計算させている最中ではなく数式を打ち込んでいる最中に止まります
0671132人目の素数さん
2009/12/06(日) 00:19:150672132人目の素数さん
2009/12/06(日) 00:57:090673132人目の素数さん
2009/12/06(日) 07:30:08確実に不良セクタの確認をするには専用のツールを使うしかない。
スタンダードなのはHD Tune。無料版と有料版があるが無料版でOK。
http://www.forest.impress.co.jp/article/2007/02/22/hdtune.html
もしこれで不良セクタが見つかった場合はHDDをフォーマットしてから
(クイックフォーマットはダメ)再インストールするのが確実。
一応OS標準の機能としてHDDのエラーチェックがありこれで修復もできるが
これは不良セクタで読めなかったファイルを勝手にダミーで置き換える上に
どのファイルを置き換えたかレポートしないのでほとんどの場合解決にならない。
0674132人目の素数さん
2009/12/06(日) 09:03:260675132人目の素数さん
2009/12/06(日) 21:08:52下記は強制振動のHPです。In[4]の変数を変えて、振幅をどんどん増加させるには、変数をどのようにすれば良いでしょうか?
http://www.eng.toyo.ac.jp/~nwada/AppliedPhysicsMathVer6/AP_ForcedOscillation_No7/5ForcedOscillationVer6.html
0676132人目の素数さん
2009/12/07(月) 21:50:03http://www.wolfram.com/products/mathematicahomeedition/
の一番下に,
Mathematica Home Edition is not licensed … …
… in the U.S., Canada, U.K., and Europe.
For other territories, please contact Wolfram Research.
とあるけど,「日本のユーザだけどHome Edition を買いたい!」
と尋ねてみた人いる?
0677132人目の素数さん
2009/12/07(月) 22:08:40半年くらい前だけど、日本でいつ売るのか聞いたことがある。
いつからという質問には答えてくれず、...not available in Japan と
そっけなかった、、。
みんなでもっと騒げば売り出すかも、、。
0678132人目の素数さん
2009/12/07(月) 22:45:34俺は、先生に買ってもらったMathematica がノートPCに
入っているから Home Edition は買わないけど、日本から
買った人いるのかどうかには関心がある。
もし、俺が本当に買うつもりなら、ダメもとで他の連中と同じ
値段で、俺にも売るようにプッシュしてはみるけど。
相手も商売だから、こちらからプッシュすれば、売るような
気がするけど。 そのために、please contact Wolfram Research
っていう一文があるんじゃないのかなぁ.
0679132人目の素数さん
2009/12/08(火) 00:56:02ただ拡大するだけなら線形だから入力信号のAを大きくするのとちゃうの?
共鳴に着目するなら Ω だけどね
0680132人目の素数さん
2009/12/08(火) 01:22:53もしかしたら、変なものをつけない?
マウスとキーボードだけにしても同じなの?
それからマウスもとりかえてみて
0681132人目の素数さん
2009/12/11(金) 18:19:37厳密解ではなく、近似解である可能性はあるのでしょうか?
RootをさらにToRadicalで計算した結果、無理数などで解が表されたら、それは厳密解ですよね??
0682132人目の素数さん
2009/12/11(金) 21:26:42Reduceは厳密解を求めるものであって近似解を求めるものじゃないから解が表示されたら、それは厳密な普通の解じゃない?
0683132人目の素数さん
2009/12/11(金) 21:35:520684132人目の素数さん
2009/12/12(土) 19:51:290685132人目の素数さん
2009/12/22(火) 22:10:47縦軸の桁が増えたりすると大きさが変わるじゃないですか
Plot[x, {x, 0, 1}]
Plot[100000 x, {x, 0, 1}]
0686132人目の素数さん
2009/12/23(水) 11:04:57↓にはプロットから座標を抽出するためにはCtrl+マウス左と書いてあるのですが実行しても何も起こりません
ttp://library.wolfram.com/howtos/coords/
ver5.2ではうまくいくのですが、何か操作法が変わったのでしょうか?
そうであれば正しい操作を教えて頂きたいです
よろしくお願いします
0687132人目の素数さん
2009/12/24(木) 12:08:40Plot[x, {x, 0, 1}]
Plot[10 x, {x, 0, 1}, PlotRange -> {{0, 1}, {0, 1}}]
0688132人目の素数さん
2009/12/25(金) 12:01:340689132人目の素数さん
2009/12/25(金) 18:29:360690132人目の素数さん
2009/12/26(土) 18:20:180691132人目の素数さん
2009/12/28(月) 00:11:27をつくってそれを呼んだときにアンダーフローがおきてしまうのですが
なにか解決方法知りませんか? アンダーフローがおきたら小数点5位まで丸め込むとかできないですかね?
0692132人目の素数さん
2009/12/29(火) 02:51:470693132人目の素数さん
2009/12/29(火) 13:36:290694132人目の素数さん
2009/12/31(木) 18:40:130695132人目の素数さん
2009/12/31(木) 18:52:260696132人目の素数さん
2009/12/31(木) 18:54:48ありがとうございます。
古いバージョンのせいか
Simplify[x (1 + y^2/x^2)^(1/2), x > 0] でうまくいきました。
0697132人目の素数さん
2010/01/01(金) 01:42:362次元上ではImageRotateでできるらしいんですが...
0698132人目の素数さん
2010/01/02(土) 01:53:30座標視点をかえる
アニメーション機能をつかえば
0699132人目の素数さん
2010/01/02(土) 03:07:110700132人目の素数さん
2010/01/02(土) 04:18:30表示その他ですこしかわるがあまり差は出ないと信じる
0701132人目の素数さん
2010/01/04(月) 20:42:510702132人目の素数さん
2010/01/11(月) 17:51:590703132人目の素数さん
2010/01/12(火) 01:31:44Plot[x^2, {x, 0, 1}, Ticks -> None]
0704132人目の素数さん
2010/01/12(火) 02:45:150705132人目の素数さん
2010/01/17(日) 07:00:010706猫は淫獣 ◆ghclfYsc82
2010/01/18(月) 01:00:56★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
小沢先生、頑張って下さい。私は最後まで味方になります。
猫
0707132人目の素数さん
2010/01/20(水) 12:55:494次元データのcontourを書きたいのです.
ControurPlot3Dを使えば良さそうなのは分かるのですが,使い方がわかりません.
陰関数f(x,y,z)=0の球面上にg(x,y,z)の大きさでcontourを描きたいんです.
http://tech.ckme.co.jp/img/sphere8cdata_matlab.png
こんな感じです
よろしくお願いします!
0708132人目の素数さん
2010/01/21(木) 22:17:47x J0(x) = 2 J1(x) [J1(x), J0(x)はベッセル関数BesselJ[n,x])
を満たすxを求めるにはどうすればいいですか?
NSolveは使えませんでした。
0709132人目の素数さん
2010/01/21(木) 23:33:05FindRoot[x*BesselJ[0, x] == 2*BesselJ[1, x], {x, 5}]
FindRoot[x*BesselJ[0, x] == 2*BesselJ[1, x], {x, 8}]
etc.
0710132人目の素数さん
2010/01/22(金) 10:27:51出来ました。ありがとうございます。助かりました。
0711132人目の素数さん
2010/01/25(月) 21:04:23たとえば Solve[{3*x^2+y==0,x+3*y^2==0},{x,y}] だと虚数解も出ます
実数解だけ表示したいのです
0712132人目の素数さん
2010/01/25(月) 21:06:00とやってみたけど Solve は Assumption なんてオプションは知らない
と言われました
よろしくお願いします
0713132人目の素数さん
2010/01/25(月) 22:22:220714132人目の素数さん
2010/01/25(月) 22:53:28Select[x/.%,Element[#,Reals]&]
0715132人目の素数さん
2010/01/26(火) 00:03:120716714 correction
2010/01/26(火) 00:10:25Select[{x,y}/.%,Element[#,Reals]&]
0717132人目の素数さん
2010/02/09(火) 21:29:05Import["Capture.avi", "Animation"]
Import::fmterr: データをvideo形式でインポートすることができません. >>
と出てくるんだけど原因に心当たりある方いませんか?
使ってるver6のドキュメントでは当然のように読み込めてるみたいだし
もしくは読み込めるvideo形式に何か制限があるのかもしれないけど
0718132人目の素数さん
2010/02/10(水) 07:46:25Maxima
歴史的には、こっちの方がふるい。
0719132人目の素数さん
2010/02/11(木) 07:51:11のはず。あんだけ高いんだから。
0720132人目の素数さん
2010/02/11(木) 08:15:12>>あんだけ高いんだから。
Home Editionは300ドル弱。そんなに高くない。
日本では買えないが、、、。
0721132人目の素数さん
2010/02/11(木) 14:01:59RotationMatrix3Dだとz-x-zのオイラー角らしく、y軸回りにならないのですが。
z軸90度回転行列で前後をはさめばいいことですが、見難くなるので。
0722132人目の素数さん
2010/02/11(木) 19:39:24RotationMatrix3Dは無い。
RotationMatrix[x,{0,1,0}]
ではどうか?
だめなら1の位置を変えて試すべし。
0723132人目の素数さん
2010/02/11(木) 20:58:08Ver.3.0なので、RotationMatrixはないみたいです。
RotationMatrix3Dで何とかするしかないですね。
0724723
2010/02/11(木) 22:48:07平面上の回転行列
{{Cos x,-Sin x},{Sin x,Cos x}}
を(u,v,w)空間の(u,w)だけに作用させる。
あとは自分で考えよう!
0725132人目の素数さん
2010/02/12(金) 12:36:20それって自分で関数を書けってことですよね?
引数をとって行列を帰す関数なんて、どうやって書くんでしょう?
0726132人目の素数さん
2010/02/12(金) 13:00:05行列の積だけなんだけどな、
これで分からなければMathematicaやめた方がいい。
0727132人目の素数さん
2010/02/12(金) 13:57:55行列たって、ListのListだろ?それを返す関数なんて別に難しくはないだろ?
0728132人目の素数さん
2010/02/12(金) 17:25:19>>727
どうも。単純にこれでよかったんですね。
Ry[s_]:={ {Cos[s], 0, -Sin[s]},{0, 1,0}, {Sin[s], 0,Cos[s]}};
Cだと行列は帰せないから配列を変数で渡さないと、
とか考えてしまって。
mathematicaは思い出したようにたまに使うだけで、
で行列を扱うのがほとんど初めてだったので、
だいぶまごつきましたけど、おかげさまで解決しました。
0729132人目の素数さん
2010/03/05(金) 09:36:35デフォルトの順序を変更するなどして、b-aのままにしておく方法があったら教えて下さい。
やりたいのは、Cのプログラムを出力することで、
-a+bと書くよりb-aと書いた方が速いと思うからです。
イメージ的には
Plus[b,-a, Order->f]
みたいにして項の順序を決定する関数を指定できたら、と思うのですが。
0730132人目の素数さん
2010/03/05(金) 09:39:26試してみれば?
0731132人目の素数さん
2010/03/06(土) 00:53:16その二つでどちらが使いやすいとおもいますか?
0732132人目の素数さん
2010/03/06(土) 11:29:14使いやすさとは関係ないけど、日本語の情報量の多さで言えば断然Mathematicaということに
なる。
0733132人目の素数さん
2010/03/06(土) 17:40:54自分がもっているメイプルは入力がマセマティカみたいにwiyswigでないから使いにくい
新しいバージョンはwysiwigみたいだが
0734132人目の素数さん
2010/03/06(土) 17:51:23メイプルは
だからカナダなんだね
0735132人目の素数さん
2010/03/06(土) 17:59:15役にたつとおもいますか?
0736132人目の素数さん
2010/03/06(土) 18:11:480737132人目の素数さん
2010/03/06(土) 18:18:16うれしいでしょう?
0738132人目の素数さん
2010/03/06(土) 18:23:260739132人目の素数さん
2010/03/06(土) 18:39:180740132人目の素数さん
2010/03/06(土) 18:57:390741132人目の素数さん
2010/03/06(土) 18:59:290742132人目の素数さん
2010/03/06(土) 19:08:23むしろ科学への応用がメインなのでは。
0743132人目の素数さん
2010/03/06(土) 19:10:54英語が読めない旧来型日本人の勉強法だとそうでしょうね。
0744132人目の素数さん
2010/03/07(日) 08:38:24そのほか思い浮かんだアイデアをほとんどその瞬間で検証できるんだから
勉強が進みますよ
0745132人目の素数さん
2010/03/07(日) 10:35:26何をどう計算すべきかを考えるのが人間の仕事
0746132人目の素数さん
2010/03/07(日) 10:41:36使っている人間が賢ければ計算機も賢くなる。
単純計算しかできないなんてことは無いと思うよ。
0747132人目の素数さん
2010/03/07(日) 13:05:160748132人目の素数さん
2010/03/07(日) 17:54:57すべてのn次(対称)行列Aの対して、性質expr(A)が成立つ
ことを、ForAllを使って調べたいのですが、
「全ての行列」ってどう表現すればよいか
どなたか教えてもらえますか。
0749132人目の素数さん
2010/03/10(水) 00:49:512DではLogLogPlotやLogLinearPlotは存在するようですが、3Dで見つけることができませんでした。
ver7を使っています。
よろしくお願いいたします。
0750132人目の素数さん
2010/03/10(水) 08:36:060751132人目の素数さん
2010/03/10(水) 21:25:22こんな感じ?「実対称行列の固有値は実数」
n = 2;
matrix = Table[a[i, j], {i, 1, n}, {j, 1, n}];
alist = Flatten[matrix];
expr[x_] := ( Element[Eigenvalues[x], Reals]);
ForAll[Evaluate[alist],
Element[alist, Reals] && matrix == Transpose[matrix],
expr[matrix]] // FullSimplify
n=2までしか解けなかったorz
0752132人目の素数さん
2010/03/11(木) 03:21:24具体的には、
x = t - sin(t)
y = 1 - cos(t)
で書かれたサイクロイドの微分をしたいと思っています
よろしくお願い致します。
0753132人目の素数さん
2010/03/11(木) 10:39:53yの値をLogで取ってしまうと軸に表示される値はLogを取った値が表示されますよね?
純粋なyの値を軸に表示させることはできるのでしょうか?
やはりyの値だけLogをとらないとだめなのでしょうか・・・
0754132人目の素数さん
2010/03/11(木) 12:45:550755132人目の素数さん
2010/03/11(木) 12:50:54軸の表示だけなら場当たり的だけどTicksを手で書いてしまう手はある
0756132人目の素数さん
2010/03/11(木) 13:16:32Plot3D[,,,,Ticks->{Automatic,LogScale[2,5],Automatic}]
とかはどう?
ここに霊がでてるけど:
http://mmi.tudelft.nl/~wouter/publications/model/texture.html
0757756
2010/03/11(木) 14:00:26本当に霊かと思ったら、
Needs["Graphics`Graphics`"]
がいるようだ。
でもLogScaleの使い方はいまいちよくわからないので
753が調べてまとめて教えておくれ。
0758132人目の素数さん
2010/03/11(木) 18:04:20実質的には
ForAll[{a,b,c},A={{a,c},{c,b}};Eigenvalues[A]//FullSimplify]
だけになる。
実際のM演算は実質的にはなにもしていないのでは?
0759132人目の素数さん
2010/03/11(木) 19:00:360760132人目の素数さん
2010/03/11(木) 19:23:13微積分の本を読んだ方が良い、、、。
0761132人目の素数さん
2010/03/11(木) 19:56:26dy/dx、d^2y/dx^2、d^3y/dx^3、・・・
と求めたいんですよ
1回微分なら自分でもできますが、高階微分となるとちょっと手計算では・・・
0762やんやん ◆yanyan72E.
2010/03/11(木) 20:21:46高階微分でもヤコビアンの微分の項が加わるだけでしょ.
式だけ手計算で出して面倒なところだけ
Mathematicaにやらせればいいんじゃない?
0763132人目の素数さん
2010/03/11(木) 21:07:39d/dx=1/x' d/dt
だから、例えば
x[t]=t-Sin[t]
y[t]=1-Cos[t]
d[f_]:=D[f,t]/x'[t]
とでも定義すれば
d^n y/dx^n は Nest[d,y[t],n]
で計算できるけど、これじゃダメなの?
0764748
2010/03/12(金) 02:29:18イメージはそんな感じです。もっと簡単に整数の行列で、
1: ∃matrix={a[i,j]}.∀i.a[i,i]=1…対角成分が1の行列は存在する
2: ∀matrix={a[i,j]}.∃i.a[i,i]=1…どんな行列も対角成分のどれかが1
を考えるとき、751のalistを使って
1.Resolve[Exists[alist, Element[alist, Integers], ForAll[i, Element[i, Integers], a[i, i] == 1]]]
2.Resolve[ForAll[alist, Element[alist, Integers], Exists[i, Element[i, Integers], a[i, i] == 1]]]
のようにすると、1はTrue、2はFalseとなってほしいのですが、
論理式が表示されるだけでTrue/Falseを計算してくれないのが
どうしてなのかわからないのです(ちなみに5版です)。
0765752
2010/03/12(金) 02:52:58ありがとうございます
私の手計算では、2回微分のとき、
d^2y/dx^2 = dt/dx d/dt(dy/dx)
= 1/(1-cost) ・ (cost (1- cost) - sint sint)/(1-cost)^2
= (cost - cos^2 t - sin^2 t)/(1-cost)^3
= (cost - 1)/(1-cost)^3
= -1/(1-cost)^2
という感じになるんですが、
mathematicaにsimplifyさせると
-1/4 (Csc(t/2))^4
になっちゃうんですよね・・・
この両者をTrueQで比較してもうまく行かないんですが、
どうしたらいいでしょうか?
0766132人目の素数さん
2010/03/12(金) 06:49:28それくらいなら手でも検算できるが(cost=1-2sin(t/2)^2だから)、
mathematicaでaとbがおなじかどうかは、
a-b//Simplify
として運良く0になればわかる。
0767132人目の素数さん
2010/03/12(金) 06:57:01なるほどwww
昔からmathematicaの三角関数の変形は融通がきかなくて困りますが
こういうチェック方法も面白いですね
0768132人目の素数さん
2010/03/12(金) 10:15:570769Fランク受験生
2010/03/12(金) 12:21:46具体的に計算機がどう動くかイメージできないいじょう数式処理ソフトもだめだとおもう。
計算機がどうやればできるかだね
ただ宣言すればできるなら、数学はかわるとおもうけど 裏打ちを支えるのが数学だから相互に発展するんだろうとおもいます。
0770132人目の素数さん
2010/03/12(金) 14:17:43Nest「d、y、2」だとy’がでてきてそれいじょうすすまないんですけど
どうしたらいいでしょうか?
0771132人目の素数さん
2010/03/12(金) 14:27:43x[t]=もx[t_]=だが、、、。
0772132人目の素数さん
2010/03/12(金) 14:30:360773132人目の素数さん
2010/03/12(金) 16:40:41D[f[t], {t, 4}]
で4回微分出来たけど
0774やんやん ◆yanyan72E.
2010/03/12(金) 16:45:030776132人目の素数さん
2010/03/12(金) 22:13:090777132人目の素数さん
2010/03/12(金) 22:49:470778132人目の素数さん
2010/03/12(金) 23:07:21-(1/4) Csc[t/2]^4
Nest[d, y[t], 3] // Simplify:
1/8 Csc[t/2]^8 Sin[t]
Nest[d, y[t], 4] // Simplify:
-(1/16) (4 + 3 Cos[t]) Csc[t/2]^10
...
Nest[d, y[t], 10] // Simplify
-(1/4096)(3193801 + 4517832 Cos[t] + 1553157 Cos[2 t] +
233370 Cos[3 t] + 11340 Cos[4 t]) Csc[t/2]^28
とできるけどな。あってるかどうか知らんが、、、。
0779761
2010/03/12(金) 23:28:18d[f_] := D[f, t]/D[x[t], t]
って書けばちゃんと動きましたよ
Simplify[Nest[d, y[t], 2]]とします
0780773
2010/03/12(金) 23:53:42こういう計算をMathematicaで済ますのは感心しないな。
計算するためにはアルゴリズムを考えなきゃいけないし、
アルゴリズムが分かれば自力で計算できる。
この場合Mathematicaは答え合わせとして使うべき。
0781132人目の素数さん
2010/03/12(金) 23:55:445回微分ぐらいしてみてください
関数形がひどいことになるんでw
0782Fランク受験生
2010/03/13(土) 00:10:44Nest[d,y[t],2] 779
d[y]=0
d[y[t]]=sin(t)/(1-cos(t))
d[y[z]]=0
tを意識しないといけない定義なんだね
一般性にかけるのかな?
0783132人目の素数さん
2010/03/13(土) 00:15:45陽にy=f(x)と書けないからね
0784132人目の素数さん
2010/03/13(土) 00:29:49にするとd[y[z],z]=sin(z)/x'(z) になっちゃうなあ
0786773
2010/03/13(土) 00:52:42Solve[x == 1 - Cos[t], t]
out=
{{t -> -ArcCos[1 - x]}, {t -> ArcCos[1 - x]}}
後者を採用すると、
y[x_] := t - Sin[t] /. t -> ArcCos[1 - x]
y[x]
out=
Sqrt[1 - (1 - x)^2] + ArcCos[1 - x]
D[y[x], {x, 5}] // Simplify
out=
(3 (-5 + 12 x - 12 x^2 + 8 x^3))/((-2 + x)^4 x^3 Sqrt[-(-2 + x) x])
でどう?
0787Fランク受験生
2010/03/13(土) 01:31:35わたしも形はちがっても差は0になりました。 眠いですね
0788132人目の素数さん
2010/03/16(火) 21:24:130789132人目の素数さん
2010/03/19(金) 17:37:21非暗号化ならうまくいくけど。誰かgmailのsmtpでSendMail送信できた人いますか。
0790132人目の素数さん
2010/04/13(火) 10:47:430791132人目の素数さん
2010/04/15(木) 05:51:260792132人目の素数さん
2010/04/15(木) 07:00:19だれか消費者庁に通報しろ!
0793132人目の素数さん
2010/04/28(水) 19:05:45高校数学なんてすっかり忘れてる(´・ω・`)
Mathematica様(と皆様)だけが頼りです
Sin[x]+Cos[x] = 1/3 のとき、Sin[x]^3 + Cos[x]^3 を求めよ
なんて、Mathematicaで解ける?
0794132人目の素数さん
2010/04/28(水) 19:38:270795132人目の素数さん
2010/04/28(水) 23:31:38三角関数の展開で出来るはずだよ。
0796132人目の素数さん
2010/04/29(木) 05:50:53すいません、出先なのでお礼だけ
GW明けに試してみます
0797132人目の素数さん
2010/04/29(木) 07:48:400798132人目の素数さん
2010/04/29(木) 09:40:26(sin(x) * cos(x) = 〜
の式があると簡単になるよね。
0799132人目の素数さん
2010/05/04(火) 20:35:350800132人目の素数さん
2010/05/04(火) 23:19:35sx^3+sy^3=1/27-sy/3+sy^2=1/27-sx/3+sx^2=1/2{(2/27)-(sx+sy)/3+1)}=77/153
0801132人目の素数さん
2010/05/04(火) 23:20:383/27
0802132人目の素数さん
2010/05/14(金) 16:39:23gcd[x_, y_] := ({x0 = x, y0 = y};
While[y0 > 0, gcd[x0, y0] = gcd[y0, Mod[x0, y0]]]; Return[x0])
こういう雰囲気でいけないかなと試行錯誤してます。
0803132人目の素数さん
2010/05/15(土) 23:55:31グローバル変数x0、y0を使い、関数自身をダイナミックに変更してるようで
ちょっとよくわからないのですが、簡単に、
gcd[x_, 0] := x
gcd[x_, y_] := gcd[y, Mod[x, y]]
ではいけないのですか?
0804132人目の素数さん
2010/05/16(日) 07:04:45そう書きたいのですが、再定義禁止(言葉が悪くかったです、オーバーロード禁止?)という条件で
やりたいです。
0805132人目の素数さん
2010/05/16(日) 07:45:43Fortranの様に手続き的に,しかも再帰呼び出しなしということですか?
だったらこんな感じでどうですか?
gcd[x_, y_] :=
Module[{x0 = x, y0 = y, r},
While[y0 > 0, r = Mod[x0, y0]; x0 = y0; y0 = r]; x0]
0806132人目の素数さん
2010/05/16(日) 07:55:27gcd[x_, y_] :=
Module[{x0 = x, y0 = y}, While[y0 > 0, {x0, y0} = {y0, Mod[x0, y0]}];
x0]
0807132人目の素数さん
2010/05/16(日) 09:21:12>>806
ありがとうございます!
0808132人目の素数さん
2010/05/19(水) 14:32:50を作りたいのですが教えてくれませんか?
(ヒント:nがpqの形になるための必要十分条件は、nの約数の個数が4、
かつnが立方数ではない(すなわち、n^(1/3)が整数でない)ことである)
ヒントがあってもわかりません。お願いします。
0809やんやん ◆yanyan72E.
2010/05/19(水) 18:40:02ヒントを使うよりも素直に2〜47の相異なる2つの素数の積のうち
100以下のものをリストアップした方がはやい気がする。
一個一個の数字の約数の個数を数えるのは骨が折れる。
0810132人目の素数さん
2010/05/19(水) 19:28:31FactorIntegerで素因素分解して、2つの素数の積でなりたち、指数が1、のを選べばできると思うけど、
そんな解き方でいいんだろうか?
a=Table[FactorInteger[i], {i, 100}];
b=Cases[a,x_/;Length[x]==2];
c=DeleteCases[b, {___, {_, x_ /; x > 1}, ___}]
0811132人目の素数さん
2010/05/19(水) 19:39:180812132人目の素数さん
2010/05/23(日) 21:13:32(2階線形、2変数)偏微分方程式を長方形の境界を持つ領域(0<x<a,0<y<b)で解きたいとき、
境界条件として、向かい合う辺での値を与えて(u(0,y)=f(y),u(a,y)=g(y))解くようなことはできないでしょうか?
0813132人目の素数さん
2010/06/02(水) 14:11:05戦闘の数字のリストを作れ。またそのうち、1が何個あるか調べよ。
わからないのですみませんが誰か教えてくれませんか?お願いします。
0814132人目の素数さん
2010/06/02(水) 19:15:490815132人目の素数さん
2010/06/03(木) 18:32:04頭を使ってもわからないから質問してるんじゃないのかな?
僕は初心者だからもちろんわからん。
0816132人目の素数さん
2010/06/03(木) 23:14:21力づくでダサくてもいいなら、
ToStringで数を文字列とし、
StringTakeで文字列の一部をとりだす。
文字列としての数を、数に変換したければ、ToExpressionでできる。
それで所望の関数f5[n]ができたら、1から100までのリストxに
Map[f5,x]とすれば数字リストが得られる。
1が何個あるかは、Countをつかえばよい。
>>815にもだが、
このような力づくの方法ならそれほど難しくはない。
mathematicaの解説本があるからそれを一度読んだ方が良い。
mathematicaのhelpも自分で読むべきだし。
0817132人目の素数さん
2010/06/04(金) 02:04:310818132人目の素数さん
2010/06/04(金) 22:58:050819132人目の素数さん
2010/06/04(金) 23:08:59しかしヘルプが充実しているので、ヘルプを見ればだいたい分かる。
また、デモンストレーション・プロジェクトを眺めていると、自分のやりたいことのヒントが見つかる。
0820132人目の素数さん
2010/06/06(日) 03:52:55本の方が画面より見やすいから
古い版のやつでも3000円位で売っていたら買いだろうね。
0821132人目の素数さん
2010/06/08(火) 19:20:09で、答えは?本読んでもよくわからん。
0822まんこ
2010/06/09(水) 04:37:01f5[n_] := StringTake[ToString[2^n], 1];
Count[Table[f5[n], {n, 1, 100}], "1"]
0823まんこ
2010/06/09(水) 04:38:04f5[n_] := StringTake[ToString[2^n], 1];
Count[Table[f5[n], {n, 1, 100}], "1"]
0824132人目の素数さん
2010/06/09(水) 04:39:330825132人目の素数さん
2010/06/09(水) 13:44:23ただし、東西南北をそれぞれ、{1、0}{ー1、0}{0、ー1}{0、1}で
表すことにする。
(ヒント:りすとdirectotions={{1,0},{-1,0},{0,-1},{0,1}}を作り、1、2、3、4の範囲の
乱数kを発生させ、directotionsのk番目の要素をとるとランダムな方角が得られる)
すいませんが練習だと思って解いてくれませんか?課題が出せないんでお願いしますm(_ _)m
0826132人目の素数さん
2010/06/09(水) 18:12:51自分でやれ
0827132人目の素数さん
2010/06/10(木) 00:11:39directotions = {{1, 0}, {-1, 0}, {0, -1}, {0, 1}};
NestList[# + Part[directotions, IntegerPart[Random[ ] 4 + 1]] &, 0, 100]
0828132人目の素数さん
2010/06/10(木) 07:12:36RandomChoiceというのもあるよ。
NestList[# + RandomChoice[directotions] &, 0, 100]
0829132人目の素数さん
2010/06/10(木) 07:18:14にきれいな答えが出てるね。
0831132人目の素数さん
2010/06/14(月) 17:32:24例えば、f4[1234]=10.
わかる人います?教えてくれませんか?
0832132人目の素数さん
2010/06/14(月) 17:52:16中学生かよ
自分でやれ
0833132人目の素数さん
2010/06/14(月) 18:48:420834132人目の素数さん
2010/06/14(月) 19:07:26別にいいじゃん。教えてやれよ。
わからないならともかく。
0835132人目の素数さん
2010/06/14(月) 19:14:490836132人目の素数さん
2010/06/14(月) 19:22:270837132人目の素数さん
2010/06/14(月) 19:25:190838132人目の素数さん
2010/06/16(水) 00:58:448の噂とかを尋ねて、たまたま通りかかったばぃ
>>836
上手いな…おまん若いじゃろ!!そんな関数知らん
かった、ありがとうなぁ
>>831
問題投げかけといて、音沙汰なしかい!
ぼてくりこすぞ〜
0839132人目の素数さん
2010/06/16(水) 01:40:080840>>831
2010/06/16(水) 13:49:460841132人目の素数さん
2010/06/16(水) 13:53:11てか教えて
正整数nのケタ数を与える関数f2を書け。たとえば、f2[40172]=5
0842132人目の素数さん
2010/06/16(水) 14:57:39問題間違えました。
0より大きく1より小さい既約分数で分母がn以下のものを
小さい方から並べてできるリストを出力する関数f5[n]を書け。たとえば、
f5[n]={1/4,1/3,1/2,1/3,1/4}
(このような数列をnに対するFarey数列といい、いろいろおもしろい性質が知られている。
たとえば、隣り合う項をa/b,c/dとすると、必ずab-cd=-1となっている。)
すいませんがお願いします。
0843132人目の素数さん
2010/06/16(水) 15:01:23かね。(未検証)
0844132人目の素数さん
2010/06/16(水) 15:02:260845132人目の素数さん
2010/06/16(水) 15:05:24f5[4]をしたら1になりました。違うかもしれません。
すいません。
0847132人目の素数さん
2010/06/16(水) 16:06:23たとえば、f[1234]=4321
r(n)=n/4となるような最小の正の整数を求めよ。
お願いします。意味わかりません。
0848132人目の素数さん
2010/06/16(水) 18:03:57死ね
0849132人目の素数さん
2010/06/17(木) 01:51:10意味わからんぞぃ!
例えばのfはrでええんかぁ
んでn/4はリバースする前のnが4の倍数なんか?
4の倍数を利用しないブルートフォース版じゃが
n = 1;
While[n/4 != FromDigits[Reverse[IntegerDigits[n]]], n++];
Print[n]
0850132人目の素数さん
2010/06/17(木) 04:11:150851132人目の素数さん
2010/06/17(木) 12:39:38ありがとうございます!!!!!!
0852132人目の素数さん
2010/06/17(木) 15:31:46演習問題の書込み→回答は
このスレッドの価値を上げますね。だれも損しない。
0853132人目の素数さん
2010/06/17(木) 18:48:190854132人目の素数さん
2010/06/18(金) 02:26:310855132人目の素数さん
2010/06/18(金) 17:34:030856132人目の素数さん
2010/06/18(金) 22:50:190857132人目の素数さん
2010/06/19(土) 06:33:42知らないバカ
0858132人目の素数さん
2010/06/19(土) 16:07:28刃傷沙汰起こす前に病院行けよ
0859132人目の素数さん
2010/06/19(土) 23:44:32上の問題はたしかに入門編レベルじゃけん、Mathematica使いには簡単じゃろ!
しかしなぁMathematica使いなら、すて台詞ばかりはいてどうするんかぃ?
そや上のf5[n]について、プログラム板でのマルチ投稿(マルチは気に入らんが…)
があって、そのレスが以下ばぃ!
-----------------------------------------------------
292 :デフォルトの名無しさん:2010/06/16(水) 15:22:37
>>291
あまり美しくないが。
f5[n_] := DeleteDuplicates[Sort[Flatten[Table[Table[Which[0 <= a/b <= 1, a/b, True, {}], {a, 0, n}], {b, 1, n}]]]]
-----------------------------------------------------
これでもええんじゃけど、思考的にうかぶのは上の形を関数3個で表わした形じゃろ?
おぅ!>>857!!
こんなん入門書0章じゃけんど応用がきかんとだめばぃ!
みんなにわかりやすぅ〜く説明したってくれんかぁ!!
いっとくけんど、くだらんとか簡単すぎて書けんとかなしじゃど!!
おりさまか?おりはココの住人でもないし、Mathematica使いでもないばぃ!
0860132人目の素数さん
2010/06/20(日) 17:47:38本屋に行っても「例題で解く〜」くらいしかいい本がないんだよね〜
0861132人目の素数さん
2010/06/20(日) 18:04:170862132人目の素数さん
2010/06/21(月) 23:55:01つまらんなぁ…まぁがんがれやぁ
つか、おりもあきたけん
0863Fラン受験者
2010/06/22(火) 02:53:21>>849
ANS=879^n12 n=0....infinity
8712 87912 879912 8799912 .......
0864132人目の素数さん
2010/06/23(水) 01:42:18>>849
と問われたけん、仁義のれすばぃ!
>ANS=879^n12 n=0....infinity
が一般論なんか?んなら突然の879
を説明せんかい!!
つか、ぶぅぶぅ〜
0865Fラン受験者
2010/06/23(水) 02:06:13879^(0)12 は 8712と読んでください。 n=1..infinity です。
8712
87912
879912
8799912
87999912
87999912
879999912
.......
証明もお願いします。
0866132人目の素数さん
2010/06/23(水) 02:36:17>証明もお願いします。
おりか?むっつりスケベの
>861
にまかす!!
こういう板はむずい…
0867132人目の素数さん
2010/06/23(水) 13:58:29{3,1,1,5,2,1,2}み適応し、{3,1,5,2}になることを確かめよ。
(Unionを使うと{1,2,3,5}になるので注意)
なんか方法ある?
0868132人目の素数さん
2010/06/24(木) 02:15:460869132人目の素数さん
2010/06/24(木) 12:20:430870132人目の素数さん
2010/06/24(木) 17:51:55文句しか付けられないのかお前は。
能力か人望のどちらか位は備えていないと人生キツいだろう。
0871132人目の素数さん
2010/06/24(木) 18:52:060872132人目の素数さん
2010/06/24(木) 19:38:08「その数は2通りの仕方で2つの正の立方数(立方数とはn^3)の和で表される
自然数のうち最小のものです」さてその数は何か?まずf[n]を作ること。
答えはたしか91になればいいと思います。
僕にとってムズイので教えてください。
0873132人目の素数さん
2010/06/24(木) 19:39:4991 = 6^3 + (-5)^3 = 4^3 + 3^3
0874132人目の素数さん
2010/06/24(木) 19:42:04「まずf[n]を作ること。」は消して。いらんこと書いた。
0875132人目の素数さん
2010/06/25(金) 01:09:57>答えはたしか91…
おぅ!おまんラマヌジャンしっとるんか〜
つか、タクシーの番号が二桁ってなんかい!
逸話も間違っとるし、-5は自然数じゃなか!
0876132人目の素数さん
2010/06/25(金) 01:23:45タクシーの番号を見てラマヌジャンが
「その数は2通りの仕方で2つの正の立方数(立方数とはn^3)の和で表される
自然数のうち最小のものです」さてその数は何か?
とね。
正の立方数で考えてくれませんか?
0877132人目の素数さん
2010/06/25(金) 01:24:440878132人目の素数さん
2010/06/25(金) 02:32:051729 = 12^3 + 1^3 = 10^3 + 9^3
んで、あとはかんがえれ
0879132人目の素数さん
2010/06/25(金) 13:16:440880132人目の素数さん
2010/06/25(金) 14:41:090881132人目の素数さん
2010/06/25(金) 16:14:120882132人目の素数さん
2010/06/25(金) 20:38:50c:\work\TEST.nbで保存したファイルを開いたときに,
c:\workを取得する簡単な方法がないでしょうか?
Export[]を使うときにフルパスを一々指定しないで済ませたいのですが…
0883132人目の素数さん
2010/06/26(土) 01:08:25..\
0884132人目の素数さん
2010/06/26(土) 04:16:20タクシーの番号を見てラマヌジャンが
「その数は2通りの仕方で2つの正の立方数(立方数とはn^3)の和で表される
自然数のうち最小のものです」さてその数は何か?
0885132人目の素数さん
2010/06/26(土) 07:59:39ラマヌジャンはもっと高いところにいた
そんなエピソードは忘れろ
0886132人目の素数さん
2010/06/26(土) 08:37:590887132人目の素数さん
2010/06/26(土) 13:18:170888132人目の素数さん
2010/06/26(土) 16:09:24NotebookDirectory[]
は?
0889882
2010/06/26(土) 20:23:05>>888
レスありがとうございます.
申し訳ありません.
Mathematicaのバージョンを書き忘れてました.
NotebookDirectory[]はまさに欲していた機能でしたが,
ver6以上でした.
私が使っているのはver5.1なのでまだ実装されていなかったようです.
とりあえずこのまま使い続けることにします.
0890132人目の素数さん
2010/06/26(土) 22:10:020891132人目の素数さん
2010/06/28(月) 02:39:060892882
2010/06/28(月) 11:07:23ver6以前でも使える方法がいくらか紹介されていました.
ttp://www.mofeel.net/1164-comp-soft-sys-math-mathematica/5574.aspx
ttp://blog.hulinks.co.jp/2010/02/mathematica-gis-2.html
これによると,例えば
SetDirectory["FileName" /. NotebookInformation[EvaluationNotebook[]] // First // ToFileName];
list = Table[{x, Sin[x]}, {x, 0, Pi, 0.1}];
Export["sine.txt", list, "Table"]
で sine.txt を現行ノートブックと同じフォルダに生成することができました.
因みにversionは
Version Number: 5.1.1.0
Platform: Windows
です.ありがとうございました.
0893132人目の素数さん
2010/06/29(火) 00:32:100894132人目の素数さん
2010/06/29(火) 01:36:33t=p^3+q^3=r^3+s^3
なら、p^3+q^3-(r^3+s^3)=0
を満たす四組が必要条件じゃろ!
んでブルートフォース ばぃ!!
つか、この手の問題ってそんなに
綺麗に解ける問題なんか?
0895132人目の素数さん
2010/06/29(火) 16:11:05Select[PowersRepresentations[#, 2, 3] & /@ Range[20000],
Length[#] == 2 &]
{{{1, 12}, {9, 10}}, {{2, 16}, {9, 15}}, {{2, 24}, {18, 20}}}
Map[Total, r^3, {2}]
{{1729, 1729}, {4104, 4104}, {13832, 13832}}
で、
1^3+12^3=9^3+10^3=1729
なんか?
0896132人目の素数さん
2010/06/30(水) 12:57:526より大きな偶数nに対し、n=p+q(p<=q)となるような素数の{p,q}を
与える関数goldbach[n]を書け。
goldbach[1000]を実行せよ。
goldbach[n_] := Table[{i, n - i}, {i, 3, n/2}]
までやったんだけど、goldbach[10]をやると、
{{3, 7}, {4, 6}, {5, 5}}になるんだ。どうすりゃ{4,6}消える?
PrimeQはどこで入れるの?最初から違う?
0897132人目の素数さん
2010/06/30(水) 13:36:04goldbach[n_] :=
Select[Table[{i, n - i}, {i, 3, n/2}], (PrimeQ[#] == {True, True}) &]
0898896
2010/06/30(水) 13:46:29すごいね。
ありがとう。マジで感謝!!
0899896
2010/06/30(水) 14:51:17例えば12321は対称数である。
正整数nに対し、それが対称数ならnを出力し、
そうでなければ、nとその逆転数を加え、それが対称数かどうか調べる。
この操作を対称数が得られるまで繰り返し、
得られた対称数を出力するような関数f[n]を書け。
たとえば、
78→78+87=165→165+561=726→726+627=1353→1353+3531=4884
なのでf[78]=4884である。
わかんないです。助けてください。
0900132人目の素数さん
2010/06/30(水) 15:02:57f[n_] := Module[{r = FromDigits[Reverse[IntegerDigits[n]]]},
If[r == n, Print[n], f[n + r]]]
0901896
2010/06/30(水) 15:21:050902896
2010/06/30(水) 16:26:10この予想が正しいか、実験して確かめよ。
ごめんなさい。お願いします。
0903132人目の素数さん
2010/06/30(水) 17:30:01While[Apply[Or, PrimeQ[n - 2^Range[Floor[Log[2, n]]]]] == True,
n += 2];
Print[n]
127
なめとんのかぁ?ぜんぜんだめじゃん!
0904132人目の素数さん
2010/06/30(水) 17:54:07De Polignacは「5以上の任意の奇数は(2^m)+p
(m>=1,pは素数)と書ける」と予想した。
この予想が正しいか、実験して確かめよ。
0905132人目の素数さん
2010/06/30(水) 18:08:51だからn=127で既に予想が崩れてるって...
予想がただしけりゃ、
n - 2^Range[Floor[Log[2, n]]] /. n -> 127
のリストには素数が少なくとも一つはあるはず
じゃろ?しかしこれがFalseでWhileが終了...
0906132人目の素数さん
2010/06/30(水) 18:46:370907132人目の素数さん
2010/06/30(水) 21:50:01プログラムbirthday[n]を書け。
birthday[50]を1000回試行して平均確率を求めよ。
確率苦手です。
0908132人目の素数さん
2010/06/30(水) 23:03:50birthday[n_] :=
Table[Random[Integer, {1, 365}], {n}] // (Length[Union[#]] != Length[#]) &
Table[birthday[50], {1000}] // (Count[#, True]/Length[#] &) // N
0909132人目の素数さん
2010/06/30(水) 23:22:010910132人目の素数さん
2010/07/01(木) 07:35:570911132人目の素数さん
2010/07/01(木) 11:58:03確率の授業で必ず習うってことはない
0912132人目の素数さん
2010/07/01(木) 19:55:26と表される。ここにa1,a2,....,akは-1,0,1のいずれかである。
たとえば25=3~3+0*3^2+(-1)*3+1
nに大志、上のようなリスト{a1,a2,....,ak}を与えるような関数を書け。
25に対してはリスト{0,-1,1}を与えるはずである。
またn=1000に対して実行してみよ。
お願いします。
0913132人目の素数さん
2010/07/01(木) 21:51:04n = 25; lis = {};
While[n > 0, m = Mod[n, 3];
lis = Prepend[lis, If[m != 2, m, -1]];
n = If[m != 2, Quotient[n, 3], Quotient[n, 3] + 1]];
Drop[lis, 1]
でどうですか?n=1000でも試してみて下さい
0914132人目の素数さん
2010/07/01(木) 23:13:27nにm=3^(k+1)+3^k...+1を加算すると、akは{-1,0,1}+1->{0,1,2}で、
一般的な3進数のためIntegerDigits[#,3]で係数を取得する。つぎにmを
引いて(各係数から1を引く結果akの範囲は{-1,0,1})nに戻したあと、
//〜で整形処理ばぃ
ex.
f[n_] := IntegerDigits[n + (3^(Floor[Log[3, n]] + 2) - 1)/2, 3] - 1 //
If[First[#] == 0, Rest[#], #] & // Rest
0915Fらん受験生
2010/07/02(金) 00:31:33n = 10000
Clear[aa, bb, cc];
aa = IntegerDigits[n, 3];
k = Length[aa] + 1; bb = Table[0, {k}]; cc = bb; Do[
bb[[k - i]] = aa[[i]], {i, Length[aa]}]
Do[Which[bb[[i]] == 2, {bb[[i + 1]] += 1; bb[[i]] = -1},
bb[[i]] == 3, {bb[[i + 1]] += 1; bb[[i]] = 0}], {i, k}]
Do[cc[[k - i + 1]] = bb[[i]], {i, k}]
n=10000 のとき
{1, 1, 1, 2, 0, 1, 1, 0, 1}ーー>{1, -1, -1, -1, -1, 0, 1, 1, 0, 1}
になりますね。
>>914 さんの方法の初等的なやりかたですが、(頭のいい人はちがうなあああ)
これはAndre Weil の整数論初歩の練習問題(証明せよという形式)に出てきますね。
0916Fらん受験生
2010/07/02(金) 01:00:03天秤と{1,3、。。。、3^n-1}gの錘群をつかえばすべてのNgが計量できることを証明せよ
のようなものでした。
シロウト向けの講義いい加減な記憶ですから、責任はとりません。
0917132人目の素数さん
2010/07/02(金) 02:35:19n=2345 {0, 1, -1, 0, 0, -1, -1}
>>915
n=2345 {0, 1, 0, 1, -1, 0, 0, -1, -1}
アラ違う?!
0918913
2010/07/02(金) 08:56:10ちょっと一本取られました。綺麗な解答ですね。
>>917
914,915では正解でした。検算
In:=
3^7 + 3^5 - 3^4 - 3 - 1
Out:=
2345
0919132人目の素数さん
2010/07/02(金) 08:57:100920912
2010/07/02(金) 13:22:320921132人目の素数さん
2010/07/02(金) 13:35:36f[n_, l_: {}] :=
If[n == 0, Rest[l],
f[Quotient[n + 1, 3], Prepend[l, Mod[n + 1, 3] - 1]]]
f[25]
{0, -1, 1}
どや?
0922132人目の素数さん
2010/07/02(金) 13:37:590923132人目の素数さん
2010/07/02(金) 13:58:07そして、その数値と求積の公式や積分を利用して得られる理論値との違いを調べてください。
すなわち答えることは以下の4点です。
(1)どうゆう図形について考えましたか?
(2)その図形の面積をモンテカルロ法で求めるとどうゆう結果がでましたか?
(3)その図形の面積を数学で求めるとどうなりますか?
(4)あなたがモンテカルロ法で計算した結果は理論値に対して何%のずれがありますか?
0924↑
2010/07/02(金) 13:59:340925Fらん受験生
2010/07/02(金) 15:02:04921のやりかた f「2345」={0, 1, -1, 0, 0, -1, -1}
917のやりかた cc= f「2345」= {0, 1, 0, 1, -1, 0, 0, -1, -1}
で921はあやまり?(けたくだりについて考え不足でした)
917は再帰ではなく、Cでもかけます。 めもり、計算速度の効率がいい。 でもつまんない><;;
>>923 はべつにCでもBasic でもいいんじゃない。 とりあえず円でも計算したら。
乱数の性質と精度と計算時間とビット長の関係でも報告したら、いいんじゃない。
0926132人目の素数さん
2010/07/02(金) 15:28:42ん?まちがっとるかぁ?
2345 = {1x3^7+0x3^6+1x3^5-1x3^4+0x3^3+0x3^2-1x3^1-1}
-> {1, 0, 1, -1, 0, 0, -1, -1} == {1, a1, a2, ... ak}
で、{a1, a2, ... ak}の{0, 1, -1, 0, 0, -1, -1}を返す
のとちがうんかぁ
0927132人目の素数さん
2010/07/02(金) 17:54:28>{a1, a2, ... ak}の{0, 1, -1, 0, 0, -1, -1}を返す
のとちがうんかぁ
それはずるい
はなしならん
0928132人目の素数さん
2010/07/02(金) 18:08:16マセじゃないといけないらしいです。
モンテカルロ法が今ひとつわからないです。お願いします。
0929Fらん受験生
2010/07/02(金) 21:25:40ss[x_, y_] = x^2 + y^2 - 1;
pp := (mx = Random[]; my = Random[]; ss[mx, my]);
tt[n_] := Table[pp, {n}]
Area[n_] := Count[Negative[tt[n]], True]/n // N
Truearea = Pi/4
ErrorDiff[n_] := Area[n] - Truearea
Table[ErrorDiff[n],{m,1000,10000,1000}]
で検討してみてください
0930132人目の素数さん
2010/07/02(金) 21:33:34はは 畳の上の水練ですな クリティカルパスの設計で修羅場をくぐっていませんね
0931918
2010/07/02(金) 22:24:04クリティカルパスの設計そのものが全く分かりません。
プログラムは昔BASICを自己流でやっただけなので、ヘルプを見ながら何とか作っている状態です。
でも、このスレはとても勉強になるので感謝しています。
0932132人目の素数さん
2010/07/02(金) 23:55:15すげー。ありがとう。Sランに改名しろよ。
0933132人目の素数さん
2010/07/03(土) 01:17:05あんなぁ、おりが>>914で『手抜き版』と謳ったけん、
921で>>913さんの純正の正面突破ロジックを、
おりなりに圧縮&再帰化した『再帰版』じゃけど!
>>925の疑義がでたけん、簡単な説明をしたばぃ!
んなら…
-------------------------------------------------
>{a1, a2, ... ak}の{0, 1, -1, 0, 0, -1, -1}を返す
のとちがうんかぁ
それはずるい
はなしならん
-------------------------------------------------
って、どこがずるくて、はなしならんのか説明せんかい!
なんかソロバン塾で、計算もしないでご高説を垂れるお方
が沸いてますが、徒な無意味エントロピー増大お疲れさま。
0934132人目の素数さん
2010/07/05(月) 08:08:12n段目が辺の長さが1の正方形であるようなピラミッドがある。
各段は長さ1の正方形のレンガが敷き詰めてある。
レンガの総計が平方数であるようなピラミッドは何段か?
お願いします。
0935132人目の素数さん
2010/07/05(月) 12:19:000936132人目の素数さん
2010/07/05(月) 13:04:280937132人目の素数さん
2010/07/05(月) 13:26:330938132人目の素数さん
2010/07/05(月) 15:22:55Mathematicaと関係ない、やぽー知恵袋で聞けば?
0939132人目の素数さん
2010/07/05(月) 16:00:42でもこの問題をMでやれと言われました。
0940132人目の素数さん
2010/07/05(月) 16:19:19っMatlab
0941Fらんく受験生
2010/07/05(月) 21:26:52n=24 のとき 70x70になる
0942132人目の素数さん
2010/07/06(火) 00:35:270943132人目の素数さん
2010/07/06(火) 00:48:31トランプのカードを
{H,1},{H,2}.....,{H,13}(ハート)
{S,1},{S,2}.....,{S,13}(スペード)
{D,1},{D,2}.....,{D,13}(ダイヤ)
{C,1},{C,2}.....,{C,13}(クラブ)
で表す。
(1)52枚のカードからランダムに5枚選び出す関数を書け。
(2)5枚の手札の(ポーカーにおける)点を下の表に従って与える関数を書け。
(3)(1)を1000回繰り返したときの平均点数を求めよ。
(表)
ワンペア 5点 ツーペア 10点 スリーカード 15点
ストレート 20点 フラッシュ 25点 フルハウス 30点
フォーカード 40点 ストレート・フラッシュ 50点
ロイヤル・ストレート・フラッシュ 100点
0944132人目の素数さん
2010/07/06(火) 00:49:210945132人目の素数さん
2010/07/06(火) 22:24:490946132人目の素数さん
2010/07/07(水) 03:53:060947943
2010/07/07(水) 10:21:000948132人目の素数さん
2010/07/07(水) 19:28:170949132人目の素数さん
2010/07/07(水) 20:13:02まずは「ポーカー 確率」で検索して、自分で考えて下さい。
0950132人目の素数さん
2010/07/07(水) 22:21:170952132人目の素数さん
2010/07/08(木) 11:50:25教えてくれないのですか?わかるなら教えてください。
0953132人目の素数さん
2010/07/08(木) 16:33:420954Fランク受験落第生
2010/07/09(金) 16:32:05すみません。
この問題は非常によい問題です。(教育の観点から)
問題の記述と解決の努力のなかでいろいろ学ぶところの多いとおもいます。 出題者の先生もそれを期待しているのでしょう。
時間がかかっても自分で試行錯誤してとかれることを進めます。
ポーカをあまりよくしませんが、ルールの表現と判定などいろんなやり方があると思います。
0955132人目の素数さん
2010/07/09(金) 17:29:59とっかかりの参考例…某ソースの流用版!
carddeck =
Flatten[Map[({{H, #}, {S, #}, {D, #}, {C, #}}) &, Range[13]], 1];
f := RandomSample[carddeck, 5]
EvaluatePokerHand[a_List] := With[{
dist = Sort /@ Map[Length, Split /@ Sort /@ Transpose[a], {2}],
straight =
MemberQ[Sort /@ Take[Partition[Range[13], 5, 1, 1], 10],
Sort[Last[Transpose[a]]]]
},
If[straight,
If[ dist[[1]] == {5},
If[Sort[Last[Transpose[a]]] == {1, 10, 11, 12, 13},
100(*"royal flush"*), 50(*"straight flush"*)], 20 (*"straight"*)],
If[ dist[[1]] == {5}, 25(*"flush"*),
Switch[dist[[2]],
{1, 4}, 40(*"four of a kind"*), {2, 3},
30(*"full house"*), {1, 1, 3},
15(*"three of a kind"*), {1, 2, 2},
10(*"two pair"*), {1, 1, 1, 2}, 5(*"pair"*), {1, 1, 1, 1, 1},
0(*"nothing"*)
]]]]
(EvaluatePokerHand /@ Table[f, {1000}]) // Mean // N
0956132人目の素数さん
2010/07/09(金) 17:37:200957132人目の素数さん
2010/07/09(金) 17:44:410958132人目の素数さん
2010/07/14(水) 13:48:10Nの倍数ばかりでてくる数列を作成してください。
そして、その数列に登場する数がすべてのNの倍数に
なっていることを証明してください。
ただし、作成する数列としては、ak=NkのようにNの倍数を含むものや、
連続するN数の積をとるといったような、
明らかにそれとわかるものは除外します。
すいませんがわかる方教えてくださいますようお願いします。
0959132人目の素数さん
2010/07/15(木) 13:05:340960132人目の素数さん
2010/07/15(木) 19:01:180961132人目の素数さん
2010/07/16(金) 13:03:41Nの倍数ばかりでてくる数列を作成してください。
そして、その数列に登場する数がすべてのNの倍数に
なっていることを証明してください。
ただし、作成する数列としては、ak=NkのようにNの倍数を含むものや、
連続するN数の積をとるといったような、
明らかにそれとわかるものは除外します。
すいませんがわかる方教えてくださいますようお願いします。
0962132人目の素数さん
2010/07/16(金) 13:13:29当たり前すぎてワケ分かりません。
0963132人目の素数さん
2010/07/16(金) 13:37:58途中まででいいんで。
0964132人目の素数さん
2010/07/16(金) 13:48:14倍数である=割り切れる=余りが出ない だから、Mod関数を使えば?
0965132人目の素数さん
2010/07/16(金) 13:56:130966132人目の素数さん
2010/07/16(金) 14:09:130967132人目の素数さん
2010/07/17(土) 02:24:26ぼけ
0968132人目の素数さん
2010/07/18(日) 04:39:310969132人目の素数さん
2010/07/19(月) 03:13:200970132人目の素数さん
2010/07/19(月) 10:02:420971132人目の素数さん
2010/07/19(月) 17:34:050972132人目の素数さん
2010/07/20(火) 18:07:46ばーか
0973132人目の素数さん
2010/07/20(火) 18:36:02Nの倍数ばかりでてくる数列を作成してください。
そして、その数列に登場する数がすべてのNの倍数に
なっていることを証明してください。
ただし、作成する数列としては、ak=NkのようにNの倍数を含むものや、
連続するN数の積をとるといったような、
明らかにそれとわかるものは除外します。
すいませんがわかる方教えてくださいますようお願いします。
レス数が950を超えています。1000を超えると書き込みができなくなります。