〓 Mathematica 伍 〓
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/11/11(金) 09:02:28.82, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,）)
! | ! .,'-.{ ! !|; |`､.}ﾞ!.! |. ! ヽ.
', ', |Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l ! | ﾞ',
ヽ､', l:!Kノ}. f:_.)iﾞi: ﾘ ! l ル
| l!iヽﾞ- ' , .!__:ﾉ ﾞ ,ﾘ l ﾘ'´
. ',|!!､ r‐┐ ｀ ノ'. /,ｲ ／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣
'i!ﾞ､ヽ､ ﾞｰ' _, ィ,:',:''´ ＜ Mathematica に関する話題はここに書いてね！
ﾞ:､ｨ､jヾｰ::: 'iﾍ .ノ',ﾘ. ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
,､- '´ ヽ､ﾞ､ { ｀＞"､
/＼＼ ', } //｀ヽ
過去スレ
〓Mathematica〓
http://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/
〓 Mathematica 四 〓
http://kamome.2ch.net/test/read.cgi/math/1197270001/ (鯖飛びでログ消滅)
http://mimizun.com/log/2ch/math/1197270001/
〓 Mathematica 5 〓
http://kamome.2ch.net/test/read.cgi/math/1285859504/
0795１３２人目の素数さん
2014/06/07(土) 18:10:56.76セミナーというかMathematica紹介のプレゼン的な内容だったから、>>794に載るような
ことをする人じゃないと思う
講師というか、いわゆるエヴァンジェリストｗなのかな
0796１３２人目の素数さん
2014/06/07(土) 19:34:29.660797１３２人目の素数さん
2014/06/07(土) 19:35:11.85>>794 の誰かだろ
0798１３２人目の素数さん
2014/06/10(火) 01:06:33.78Timing関数で時間を図るとしてどういう風にすればいいんでしょうか？
入力値と出力値は決まってます
0799１３２人目の素数さん
2014/06/10(火) 18:04:17.920800１３２人目の素数さん
2014/06/12(木) 18:09:41.01今日発表かな。
0801１３２人目の素数さん
2014/06/12(木) 18:38:33.030802１３２人目の素数さん
2014/06/12(木) 19:48:02.060803１３２人目の素数さん
2014/06/12(木) 20:37:29.260804１３２人目の素数さん
2014/06/13(金) 06:02:11.28https://pbs.twimg.com/media/Bp6wm6OIcAAoDvm.png
https://pbs.twimg.com/media/Bp1Z6R1CAAA59Wl.png
https://pbs.twimg.com/media/Bpwg88YIMAAxbf6.png
https://pbs.twimg.com/media/BprLc16CQAANcZ4.png
https://pbs.twimg.com/media/BpW4iGEIUAAEG7m.png
https://pbs.twimg.com/media/BooRdTsIIAAqtNO.jpg
0805１３２人目の素数さん
2014/06/25(水) 01:15:13.390806１３２人目の素数さん
2014/06/27(金) 22:39:56.520807１３２人目の素数さん
2014/06/28(土) 02:36:08.71自己レスだが、単位系換算サポートの上での議論でゲージ原理やら接続、テンソルやらが議論にのぼったんだろうけど
ゲージ原理が時計表示のゲージに曲解されて実装されたのがゲージのような気がする。
0808１３２人目の素数さん
2014/07/07(月) 17:06:30.13Mathematica では違うのな。
逆とか気持ち悪いぜ……
0809１３２人目の素数さん
2014/07/07(月) 22:20:19.41下記はどこが悪いでしょうか？
x = 10;
If[x == Integer, Print["Ok"];
0810１３２人目の素数さん
2014/07/07(月) 23:00:22.58IntegerQでは如何
0811１３２人目の素数さん
2014/07/08(火) 09:16:31.32If[IntegerQ[x] == True, Print["Ok"]];
で出来ました。助かりました。有難う御座いました。
0812１３２人目の素数さん
2014/07/09(水) 01:13:15.38あほ
あっホウ
0813１３２人目の素数さん
2014/07/12(土) 01:03:23.37learning, computational geometry, geographic computation, and
device connectivity--as well as deepening capabilities and
coverage across the algorithmic spectrum.
For a complete list of new features and enhancements, including
multiple undo, visit:
multiple undoｷﾀｺﾚ
0814１３２人目の素数さん
2014/07/13(日) 13:45:21.82Mathematica の日本語リファレンス、ページが死んでないか？
0815１３２人目の素数さん
2014/07/15(火) 19:04:10.290816１３２人目の素数さん
2014/07/17(木) 06:59:13.84Mathematicaをご利用いただき，誠にありがとうございます．
Mathematica 10日本語版が今週リリースされる予定です
0817１３２人目の素数さん
2014/07/17(木) 17:30:30.330818１３２人目の素数さん
2014/07/17(木) 18:40:35.11A new version of Mathematica Japanese Edition for Students for Sites has been released and is available
through your organization's site license. Your existing license has been upgraded, as follows:
0819１３２人目の素数さん
2014/07/18(金) 01:22:57.39出たねバージョン１０
なんかおもしろいこと発見したらみんな報告してや
0820狸 ◆2VB8wsVUoo
2014/07/18(金) 14:46:20.03>46 名前：KingMathematician ◆LoZDre77j4i1 ：2014/07/17(木) 07:43:40.73
> 面が積もる, つまり 3 次元測度.
>
>4 名前：KingMathematician ◆LoZDre77j4i1 ：2014/07/17(木) 22:08:46.17
> y=800x-20x^2/円.
>
>6 名前：KingMathematician ◆LoZDre77j4i1 ：2014/07/17(木) 22:12:59.46
> Beside, I was strict.
>
>42 名前：KingMathematician ◆LoZDre77j4i1 ：2014/07/17(木) 22:17:01.62
> とりあえず, 応用数学を修得すれば良かろう.
>
0821１３２人目の素数さん
2014/07/18(金) 17:26:42.32[ESC]+lambda とか打つときに候補が出て補完されて。
0822１３２人目の素数さん
2014/07/18(金) 20:13:11.290823１３２人目の素数さん
2014/07/18(金) 20:57:14.340824１３２人目の素数さん
2014/07/18(金) 21:21:30.750825１３２人目の素数さん
2014/07/18(金) 23:46:53.45復活してるな。バージョンアップのためのメンテナンスだったのか。
0826１３２人目の素数さん
2014/07/19(土) 00:09:14.360827１３２人目の素数さん
2014/07/19(土) 08:34:03.13https://pbs.twimg.com/media/Bs1SW5JCIAAwGWB.png:large
0828１３２人目の素数さん
2014/07/19(土) 11:21:52.79より時間がかかってしまうのは、
積分だから今までより多くの変数変換の適用を試みてると考えればまだわかる
ただ、問題は最後の結果。
解ける問題が解けなくなってしまっているのは。。。
ウルフラムに報告しろ！
0829１３２人目の素数さん
2014/07/19(土) 11:48:05.18V10でこの問題解けるようになったか誰か試してみてくれる？
http://www15.atwiki.jp/volume_calculation/
俺が作った計算用のコマンドは以下（他に計算方法があれば何でもOK）
n = 2;
Integrate @@
Join[{Boole[
Sum[c[i]^2, {i, n}]^(-1/2) >=
Sum[c[i]^2/a[i]^2, {i, n}]^(-1/2) + 1]},
Table[{a[i], 0, 1}, {i, n}], Table[{c[i], -1, 1}, {i, n}]]
これをV9で計算すると数値の結果が返ってこなくて、
>>827のV10の最後の結果みたいな感じになる。
n=1にするとちゃんと正しい答え1が返ってくる。
また、NIntegrateにすると、いろいろ警告はでるけど、
正しい答え(5/4)π-8/3の近似値1.26071と出てくる。
V10で解けるようになったかが知りたい。
あと、Mathematica関係なくこの問題解ける人いる？
0830１３２人目の素数さん
2014/07/19(土) 18:24:00.070831１３２人目の素数さん
2014/07/20(日) 14:11:42.16同じ計算させてるのに1回目と2回目で結果が違う。
ちなみにV9なんだけど、V10もってる人ふるまい同じか試してみてくれない？
n = 2;
ss = Array[s, n];
u[n_] := s[n]/Sqrt[1 - Sum[s[i]^2, {i, n - 1}]]
v[1][ss_] := ss[[1]]
v[n_][ss_] :=
u[n] + Integrate[
v[n - 1][ss[[;; n - 1]]/Sqrt[1 - ss[[n]]^2/a[n]^2]], {a[n], u[n],
1}, Assumptions ->
And @@ Thread[0 < ss] && Sum[s[i]^2, {i, n}] < 1]
v[n][ss](*1回目*)
v[n][ss](*2回目*)
0832１３２人目の素数さん
2014/07/20(日) 14:47:45.85V9だけど結果同じ
0833832
2014/07/20(日) 14:48:26.59○ V8
0834１３２人目の素数さん
2014/07/20(日) 15:09:22.91http://up3.viploader.net/omoshiro/src/vlomoshiro037682.png
0835１３２人目の素数さん
2014/07/20(日) 15:31:30.47確認ありがとう。
「結果同じ」っていうのは
●俺と結果が同じ、つまり、計算結果は両者で違う
のか
●計算結果が両者で同じ
のどちらの意味かな？
>>834
ワオ！正しい方に結果が落ち着いてるね
バグ報告しようかと思ったけどV10でなおってるならいいや
確認ありがとう。
V8とV10は正しく計算されるのにV9だけおかしいのかな。。。
http://up3.viploader.net/omoshiro/src/vlomoshiro037683.png
0836１３２人目の素数さん
2014/07/20(日) 15:36:08.020837１３２人目の素数さん
2014/07/20(日) 15:38:48.86一回目と二回目で等しい結果が出たという意味ね。
0839１３２人目の素数さん
2014/07/20(日) 15:49:10.80確認ありがとう。
俺もV9って書いてたけど正確にはV9.0.1で、
パソコンがおかしくなってんのかと思って再起動してみたけど、
結果は同じ（答えが変わる）だった。
ちなみにWindows版ね。
同じバージョンなのにふるまいが変わるとか一体どうなってんだｗ
0840832
2014/07/20(日) 15:51:38.14いや，よくみたら>>834とは符号逆だな
0841832
2014/07/20(日) 15:56:15.479.0.1.0
全然検証してないけど見た目の結果は確かにV8と違うな
0842１３２人目の素数さん
2014/07/20(日) 16:24:47.23V8は間違ってる方に落ち着いてるねｗこれはたちが悪いｗ
V9.0.1.0の結果は俺のと違って、ちゃんと正しい答えに落ち着いてるのね
$Version俺もやってみたけど832さんのとまったく一緒のが出てきた
ちなみにただしいv[2]は
v[2][s1,s2]=s1√(1-s2^2)+s2√(1-s1^2)
符号はマイナスのやつが正しくて、
第一項と第二項を共通項くくりだして約分すれば上式になることが確認できる。
ちなみにn=2に限らず結果は対称式になる
0843１３２人目の素数さん
2014/07/20(日) 17:07:52.5910でQuit[]してスタート->初回が異なる
キャッシュをリセット->初回の値が返る
http://imgur.com/Gxup1AK
0844１３２人目の素数さん
2014/07/20(日) 17:23:14.18Mathematicaを起動してない状態から起動して
>>831をコピペして実行するとどうなる？
答えが変わるんならV10もV9と同じ振る舞いということになるな
>>834のV10はMacっぽいから、
Macだと両者正しい値がでるのか、
もしくは>>834が他のノートブック開いててそこで計算させてた状態で
コピペ→実行だとキャッシュがらみでうまくいったのかも
0845843
2014/07/20(日) 17:58:18.2910起動 -> 831コピペ&実行 -> 答えは異なる(初回プラス)
なお、
10では仮定Thread[0 < ss]に等号を入れても、8&9のように
ConditionalExpressionは返らず、共にマイナス側が返って
きます。よくわからん。。。
0846１３２人目の素数さん
2014/07/20(日) 22:30:51.74だけどせっかく一般的に作ったのにn=3に関してはだめっぽい（時間かかりまくる）
なんかうまい方法ないかな
n = 2;
ss = Array[s, n];
u[n_] := s[n]/Sqrt[1 - Sum[s[i]^2, {i, n - 1}]]
v[1][ss_] := ss[[1]]
v[n_][ss_] :=
v[n][ss] =
u[n] + Integrate[
v[n - 1][ss[[;; n - 1]]/Sqrt[1 - ss[[n]]^2/a[n]^2]], {a[n], u[n],
1}, Assumptions -> And @@ Thread[0 < ss]]
v[n][ss]
v[n][ss] = Simplify[Apart[v[n][ss]], Sum[s[i]^2, {i, n}] < 1]
V[n] = Evaluate[v[n][ss] /. Table[s[i] -> Slot[i], {i, n}]] &
jacobian[n_] :=
r^(n - 1) Product[Sin[\[Theta][i]]^(n - i - 1), {i, n - 1}]
ss = Table[(1 - r) Product[Sin[\[Theta][i]], {i, j - 1}] If[j == n, 1,
Cos[\[Theta][j]]], {j, n}]
w[n] = 2^n Integrate[V[n] @@ ss jacobian[n], {r, 0, 1},
Sequence @@ Table[{\[Theta][i], 0, \[Pi]/2}, {i, n - 1}]]
0847１３２人目の素数さん
2014/07/20(日) 22:34:50.2710でもなおってないんだ
ま、等号を含めるかどうかで体積は変わらないはずだから
そもそもそれもバグだよねｗ
V8が自信満々に間違った答えを出力してることからも
Integrateについてはあんまり信用しない方がいいね
0848１３２人目の素数さん
2014/07/21(月) 02:22:57.43> >>834のV10はMacっぽいから、
俺834だけど、起動した直後にコピペして計算結果出しただけだからキャッシュとか関係ないだろう。
0849１３２人目の素数さん
2014/07/21(月) 21:28:48.19ttp://www15.atwiki.jp/volume_calculation/
はずかしながらノートとnbファイルをアップしました
漸化式を展開することによってv[n],w[n]を
それぞれ部分に分解するというアプローチをとった。
厳密解に関しては
v[3]は全部求まったけどw[3]は一部分だけ求まった
1/12 \[Pi] (-5 + 3 \[Pi] - 6 Log[2])
近似値に関しては時間をかければいくらでもいけるっぽい
nbファイルの最後にグラフがある
0850１３２人目の素数さん
2014/08/03(日) 07:03:19.15変なコーディングしてんのかね
0851１３２人目の素数さん
2014/08/15(金) 11:37:42.30下記ＨＰのプログラムは、コンプトン散乱を計算したものです。
通常は、
sl[q] = (?u[0]*q0 + ?u[1]*(-q1) + ?u[2]*(-q2) + ?u[3]*(-q3) + ms);
で計算して
{p1 -> 0, p2 -> 0, k0 -> p3, k1 -> 0, k2 -> 0, k3 -> -p3, q0 -> p0, q1 -> p3*Sqrt[1 - z^2],で、纏めます。
すると、正常の解が得られます。
ここで、
sl[q] = (?u[0]*q0^0.5 + ?u[1]*(-q1^0.5) + ?u[2]*(-q2^0.5) + ?u[3]*(-q3^0.5) + ms);
{p1^0.5 -> 0, p2^0.5 -> 0, k0^0.5 -> p3, k1^0.5 -> 0, k2^0.5 -> 0, k3^0.5 -> -
のように変更します。
すると、全く違う答えになってしまいます。どこが間違っているのでしょうか？
計算時間は５分程度です。
http://www.geocities.jp/dirac_equation/index.htm
0852１３２人目の素数さん
2014/08/15(金) 11:51:50.46なんで 1/2 じゃなくて 0.5 とか使うのかね
0853１３２人目の素数さん
2014/08/15(金) 16:47:49.98>どこが間違っているのでしょうか？
自分が何を計算してるのかを全く理解せずに式だけいじろうとしてるところ
0854１３２人目の素数さん
2014/08/15(金) 17:30:49.00＞自分が何を計算してるのかを全く理解せずに式だけいじろうとしてるところ
ご指摘の通りです。
マンデルスタム変数の導出方法を教えて下さい。
p0→(s+m^2)/(2Sqrt[s]),
p1→0,
p2→0,
p3→(s-m^2)/(2Sqrt[s]),
k0→p3,
k1→0,
k2→0,
k3→-p3,
q0→p0,
q1→p3*Sqrt[1-z^2],
q2→0,
q3→p3*z,
j0→p3,
j1→-p3*Sqrt[1-z^2],
j2→0,
j3→-p3*z,
z→1+t/(2p3^2),
t→2m^2-s-u
0855１３２人目の素数さん
2014/08/15(金) 17:58:24.78>>745の
>マンデルスタム変数：http://www.geocities.jp/dirac_equation/desnhishoumetu.pdf の(8.55)式
の通り
0856１３２人目の素数さん
2014/08/15(金) 18:40:46.13>マンデルスタム変数：http://www.geocities.jp/dirac_equation/desnhishoumetu.pdf の
＞通り
後で、再度勉強します。
思い返しますと、マンデルスタム変数のmathematicaコードを教えて頂く前に、自分なりに、不完全ですがコードを作って計算しておりました。それも、見直します。
コンプトン散乱で、
p0,p1,k0,k1, k2,k3などが、すべて２乗倍になったときの計算結果は
（１）q0をq0^2などに変更したコンプトン散乱その１
（２）q0をq0^2などに変更したコンプトン散乱その２
（３）どちらも間違い
の内のどれでしょうか？
http://www.geocities.jp/dirac_equation/index.htm
のcompton r1に計算結果を追加しました。
0857１３２人目の素数さん
2014/08/15(金) 21:14:47.19デタラメやる前にちっとは勉強したら
0858１３２人目の素数さん
2014/08/15(金) 21:43:31.49数学として考えた場合、答えは、どうなると考えられますか？
私は
（２）q0をq0^2などに変更したコンプトン散乱その２
が正しいと考えます。
p0,p1,k0,k1, k2,k3などが、すべて２乗倍になるだけで、マンデルスタム変数を使って変換する式は、そのままのはずです。
すると、ルートの場合は、どうしたら計算できるのか？
がわからないです。
0859１３２人目の素数さん
2014/08/15(金) 22:55:53.13ベクトルやベクトルの内積や合同変換群や不変量などの数学的な性質を
一切無視した置き換えだからデタラメだと言われてるんだが。
「ローレンツ不変量」とか知らんのだろう？
0860１３２人目の素数さん
2014/08/16(土) 07:43:49.82＞「ローレンツ不変量」とか知らんのだろう？
問題を整理します。
ローレンツ変換に対して、不変にならないのは、どの量でしょうか？
１．p0,p1,k0,k1, k2,k3などが、すべて２乗倍になったときのこの式でいうｓl(q) 等
２．マンデルスタム変数等の変換の仕方
まず、１ですが、不変であるような気もします。クラインゴルドン方程式でローレンツ変換に対して、不変になるからです。
２ですが、デタラメデス。導出方法が理解できず、適当に変換方法を、２乗しています。
そこで、いままでの計算は忘れて、（２乗倍とか単なるルート倍とかはやめて）、物理的に意味のある計算を提案します。
①p0を、p01＾0.5とp02＾0.5の積として分解して計算を進める。
②p01＾0.5、p02＾0.5のマンデルスタム変数等の変換を行い
③最後にp01＾0.5＝p02＾0.5＝p0＾0.5として
④通常のコンプトン散乱の結果が得られることを確認する
これなら、「ローレンツ不変量」その他、すべて、満たしているはずです。
「なんで、わざわざそんな面倒なことをするのか？」と思われるかもしれませんが、自動車の分解が面白いのと同じで、特に意味はなく、面白そうだからです。
0861１３２人目の素数さん
2014/08/16(土) 11:35:15.48下記を計算してみました。
ＨＰのcompton 途中
です。
コードのどこが悪いでしょうか？
①p0を、p01＾0.5とp02＾0.5の積として分解して計算を進める。
②p01＾0.5、p02＾0.5のマンデルスタム変数等の変換を行い
③最後にp01＾0.5＝p02＾0.5＝p0＾0.5として
④通常のコンプトン散乱の結果が得られることを確認する
0862１３２人目の素数さん
2014/08/16(土) 14:03:29.920863１３２人目の素数さん
2014/08/16(土) 16:11:01.71＞Mathematicaにはデバッガはないの？
バグをとるのではなく、
①理論を本当に理解して正しい計算式をコードとして入力しているか
②Mathematicaでちゃんと答えを導出できるテクニックを持っているか
になると思います。
更に下記を計算してみました。
ＨＰのcompton 更に途中
もう一歩のような気がします。消えてほしい変数が残ってしまいます。
しかし、消えてほしい変数が消えても答えが間違っているかもしれませんが、
0864１３２人目の素数さん
2014/08/16(土) 22:40:02.38おまえのスペックは？
高校生が夏休みの自由研究でもやってんの？
Mathematicaはある答えを導出するために計算式入力して
はいMathematicaさん計算してください！って使い方だけじゃないのよ
開発元のウルフラムもV10からは「ウルフラム言語」って読んでるくらいで、
普通のプログラミング言語のようにプログラムすることも可能。
ただ、個人的に思ってることはMATLABとかはデバッガがあるのに、
Mathematicaは基本的に机上デバッグするしかなくて効率わるいなぁと。
0865１３２人目の素数さん
2014/08/16(土) 22:46:06.70信じられんmapleでもデバガーあるというのに
やっぱりmaple>>>>>>>>mathematicaだろうな
0866１３２人目の素数さん
2014/08/16(土) 22:59:55.74Mapleはデバッガあるんだ！
Mapleの使い勝手ってMathematicaと比べてどう？
俺も昔Maple使ってたんだがV6とかぐらいで今はだいぶ進化してるんだろうな
その当時のMapleの印象としては
●ヘルプがくそすぎる
●しょうもない関数が多すぎ（複利計算の公式通りの関数とかｗ）
って感じだったからMathematica使うようになったんだよね
その当時デバッガがあったかは覚えてないけど、
最近のMaple vs. Mathematicaって公平にジャッジしてどうなんだろ。
ヘルプもちゃんと翻訳されてマシになってきてるかもしれんな。
0867１３２人目の素数さん
2014/08/17(日) 00:29:21.56>①理論を
全然理解してないって昔から指摘されてるのに、全く進歩しないねぇ
0868１３２人目の素数さん
2014/08/17(日) 01:29:58.20mapleは最近、ドキュメントモードとワークシートモードの２種類が
あるようになったよ。ドキュメントモードはワープロみたいな文章が
作れる。ワークシートモードは従来のコマンド入力するモードです。
クリッカブルマスと言って、コマンド入力しなくても計算できる機能
もあるんだよ。とりわけドキュメントモードの時にやくだつ。
ヘルプはmaple18で完全日本語化したよ
ヘルプについては使いにくいけど、mathematicaも変わらないと思う。
mapleはコマンドについての奥が深いと感じる。mathematicaの方が
コマンドを書きやすい。
公平にジャッジするとどうなるかはわからない。
mapleはプログラミング言語に近いけど、mathematicaはそうではない。
0869１３２人目の素数さん
2014/08/17(日) 03:10:24.32mapleもうV18までいってんだｗ
昔はmathematicaと同程度のバージョンアップだったのに
最近は開発サイクル変えたのかな
クリッカブルマスはmathematicaで言えばパレット入力みたいなもんかな？
ドキュメントモードに相当する機能（見栄えを美しくする）も
mathematicaで出来た気がする
ヘルプやっと日本語化されたんだねｗ
V18って最新版だから最近ようやくって感じか
使ってた当時、使いづらいなと思ったのが、
ヘルプ上で実行できなかったこと。
これも改善されてるのかな。
昔使ってた時のなごりで今でもメルマガが届くんだが
Mapleは応用に力入れてるのかなって感じがする
MapleSimとかをしきりにアピールしてて、
確かにこれに相当するMathematicaのアドオンは知らない
（あるのかもしれないけど知名度低い）
モデルベースデザインをうたって産業界を巻き込もうとしてる感じ
トヨタとかではけっこう実績もあったりするらしい
mathematicaは産業界への応用に力いれてる気配はなく、
純粋な学問としての数学を突き進んでる感じ
あとは環境整備とかにも力いれてるね
だんだん方向性に違いが出てきてるから、
自分にあったものを、って感じかな
0870１３２人目の素数さん
2014/08/17(日) 08:55:31.22y = Sqrt[(m^4 + 2*m^2*s - s*(s + 2*u))/p3a]*Sqrt[(m^4 + 2*m^2*s - s*(s + 2*u))/p3b];
y1 = Simplify[y //. {p3a^(1/2)*p3b^(1/2) -> p3}];
を実行させると
Sqrt[(m^4 + 2*m^2*s - s*(s + 2*u))/p3a]*Sqrt[(m^4 + 2*m^2*s - s*(s + 2*u))/p3b]
となります。
{p3a^(1/2)*p3b^(1/2) -> p3}
が効かないです。
これを効かせて
m4+2 m2 s-s (s+2 u)/p3
とさせるには、
{p3a^(1/2)*p3b^(1/2) -> p3}
をどのように修正すれば良いでしょうか？
0871１３２人目の素数さん
2014/08/17(日) 11:50:47.63いや、パレット入力とはちがう。
mathematicaだとパレットのギリシャ文字のパイ
をクリックすると入力行にそのギリシャ文字が現れる。
他方mapleだとドキュメントモードの場合は
ギリシャ文字のパイが同様のことをして現れる。
他方、ワークシートモードの場合は
パレットのギリシャ文字のパイをクリックしても
入力行に現れるのはコマンドのPiです。
クリッカブルマスというのはそうじゃなくて、
たとえば多項式の一部をまとめたい場合に、
そこをドラッグする。するとどう変形したいかの候補が
現れてそれを選択すると、もとの式全体の、しかし
一部分がまとめられた、式が現れる。
maplesimは持ってないから知らない。
mathematicaはmathematicaalpha？というのとか、
raspberrypiへの無償提供とかが話題になってるよね。
もともとmapleはAlgol70?だかを参考にして作っているから
プログラミング指向が強いんだよ。mathematicaのプログラミング
は普通のプログラミング言語とかなり違うでしょう？
そういうのが嫌だったな。それにプログラミングの本は多くないし。
0872１３２人目の素数さん
2014/08/17(日) 11:59:55.79ヘルプ上での実行はできなかったんじゃなかったかな？
だけど、それについては慣れだよ。mathematicaに慣れると
ヘルプ上での実行はあったほうがいいとおもうし、mapleに
慣れるとそんなものいらないよ、別にワークシート開いてコピー
ペーストで実行すればいいから、と思うようになるし。
バージョンアップのペースが昔とくらべてどうとかは知らない。
ドキュメントモードだけど、数式計算もできるワープロレベル
だよ。自分はワークシートモードで使ってるけど。
mathematicaは最初の頃のバージョンについてはかなりおそまつ
だったというのはよく聞く。mapleはmapleVあたりで開発者の
一人がスピンオフさせられてしまったんだよね。それでmapleの
魅力が減ってしまった。mathematicaの商業的な成功は宣伝に力を
入れたことだ、悪貨が良貨を駆逐するようにmathematicaが
mapleよりも売り上げを大きく上回るようになったともよく言われるよね。
自分はmapleが好きなんで、これからもmapleの開発が続いてほしいと
おもってる。mathematicaを昔使っているときと比べてもmapleの勉強を
するときは遥かに楽しいよ。奥が深いとおもうから。
0873１３２人目の素数さん
2014/08/17(日) 12:04:31.31maple6だとデバッガはあったはずだよ。４についての本をみたけど
それにもデバッガが出て来てたから。
0874１３２人目の素数さん
2014/08/17(日) 12:45:28.82答えが一致しません。
原因不明。
0875１３２人目の素数さん
2014/08/17(日) 13:29:15.68http://www.geocities.jp/dirac_equation/index.htm
の
答えが一致しない
の
Compton r3
を計算して頂けないでしょうか？
8*(m^4 - s*u + m^2*(3*s + u))
が正解なのですが
4*m^2*(14*m^2 - m^4/s + s + 2*u)
となります。
Mathematica９を使用しているのですが、他のバージョンでしたら、ちゃんと計算するかも？その可能性は、ほとんどないのですが、何回見直しても、原因がわからんのです。
0876１３２人目の素数さん
2014/08/17(日) 16:46:21.48多分、結果はver9と同じとは思いますが、今からver10trialをインストールして試してみます。
0877１３２人目の素数さん
2014/08/17(日) 21:38:58.98下記を教えて下さい。
①ある式ｙ１があり、その変数を別の変数に変換させました。
②別の変数を、元の変数に戻そうとしました。その式をｙ２とします。
③ｙ１とｙ２の表示は同じになるはずですが、異なります。
④しかし、見かけは異なるのに、ｙ１-ｙ２＝０になります。
ｙ１とｙ２の表示を同じにするには、元に戻す変換を、どのようにすれば良いでしょうか？
http://www.geocities.jp/dirac_equation/index.htm
の
変換
に、ファイルが入っています。
0878１３２人目の素数さん
2014/08/18(月) 12:53:52.630879１３２人目の素数さん
2014/08/18(月) 20:04:48.44有難う御座います。
おっしゃる通り、やってみましたら、出来ました。
感謝、感謝、感謝で御座います。
ご教示頂いた法則は、全く知らなかったです。
今まで、mathematicaの変換ができなくて、よく悩んでいたのですが、これが原因だったかもしれません。
昨日からずーと悩んでいたのですが、非常に助かりました。同時に、「もし教えて頂いてなかったら、一生解決しなかったのかも？」と「ぞ～」としました。
http://www.geocities.jp/dirac_equation/index.htm
に出来た計算結果を入れております。
本当に有難う御座いました。
0880１３２人目の素数さん
2014/08/18(月) 20:41:02.25すません。出来た思ったのですが、私の早とちりでした。
再度、やってみますと、結果は変わりませんでした。
どこが間違っているでしょうか？再度、ご教示頂きましたら幸甚です。
http://www.geocities.jp/dirac_equation/index.htm
0881１３２人目の素数さん
2014/08/18(月) 21:05:29.110882１３２人目の素数さん
2014/08/18(月) 21:21:56.88http://www.geocities.jp/dirac_equation/index.htm
やはり、変数を、クリアまたは初期化すれば、出来るんでしょうね？
0883１３２人目の素数さん
2014/08/18(月) 22:11:46.44ルールがp0a -> p0/p0b -> (p0a p0b)/p0b -> p0aとか…
0884１３２人目の素数さん
2014/08/18(月) 23:26:03.61お返事有難う御座います。
そういうことで御座いますか、、、やっと解りました。
何度も何度も、すいません。
今度こそ、感謝、感謝、感謝で御座います。
計算結果は、ＨＰに入れました。
更に、すいませんが、新しい質問です。
誠に勝手ながら861に戻らせて頂きます。
①p0を、p01＾0.5とp02＾0.5の積として分解して計算を進める。
→終わりました。
②p01＾0.5、p02＾0.5のマンデルスタム変数等の変換を行い
③最後にp01＾0.5＝p02＾0.5＝p0＾0.5として
④通常のコンプトン散乱の結果が得られることを確認する
②以降の計算結果をＨＰ「新しい質問です」に入れたのですが、通常のコンプトン散乱の結果が得られません。どこが間違っているのでしょうか？
夏休みの自由研究と思って付き合って頂きましたら幸甚です。
0885１３２人目の素数さん
2014/08/19(火) 22:23:37.23答えは、y1^0.5=(8*(m^4-su+m^2(3s+u))^0.5
ですが、導出方法が不明ですね。
0886１３２人目の素数さん
2014/08/19(火) 22:34:12.11f=8*(m^4 - s*u + m^2*(3*s + u))
となる(p=p,q,k,j,i=0,1,2,3)。
piaとpibのruleでは、piを纏めると
pi -> pia pib -> ci(m,s,u)^2
となるため、fとはならない。
piaとpibのruleの対称性から、仮に
一方をpib -> 1とするとfがえられる。
>p0を、p0a＾0.5とp0b＾0.5の積として分解...
ならばp0 -> p0a p0bはおかしくないですか?
0887１３２人目の素数さん
2014/08/20(水) 20:28:59.47＞ならばp0 -> p0a p0bはおかしくないですか?
その通りで御座います。
また、
答えは、y1^0.5=(8*(m^4-su+m^2(3s+u))^0.5
にならないことにも気がつきました。
私が求めていたことは、下記のように、修正すれば得られることも気づきました。（勝手にmathematicaで次の計算は出来ないだろうと思い込んでおりました。やってみますとちゃんと計算してくれました。）
sl[q] = (?u[0]*q0 + ?u[1]*(-q1)^(1/2) + ?u[2]*(-q2)^(1/2) + ?u[3]*(-q3)^(1/2)^(1/2) + ms);
sl[p] = (?u[0]*p0 + ?u[1]*(-p1)^(1/2) + ?u[2]*(-p2)^(1/2) + ?u[3]*(-p3)^(1/2)^(1/2) + ms);
sl[k] = (?u[0]*k0 + ?u[1]*(-k1)^(1/2) + ?u[2]*(-k2)^(1/2) + ?u[3]*(-k3)^(1/2)^(1/2));
sl[j] = (?u[0]*j0 + ?u[1]*(-j1)^(1/2) + ?u[2]*(-j2)^(1/2) + ?u[3]*(-j3)^(1/2)^(1/2));
貴重なご指摘に感謝致します。
0888１３２人目の素数さん
2014/08/22(金) 19:06:39.75ファイルフォーマットを変えるのが
うざい
新しくすると古いバージョンで読めない
余計なことをするな
0889１３２人目の素数さん
2014/08/22(金) 19:08:44.030890狸 ◆2VB8wsVUoo
2014/08/22(金) 19:09:21.50>4 名前：KingMathematician ◆LoZDre77j4i1 ：2014/08/21(木) 19:27:24.63
> とりあえず, vector analysis から修得しようか.
>
狸
痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕
逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢
痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕
逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢
痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕
逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢
痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕
逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢
痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕
逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢
痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕
逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢
痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕
逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢・逮捕・痴漢
0891１３２人目の素数さん
2014/08/24(日) 17:17:42.70jacobian[n_]:=r^(n-1) Product[Sin[\[Theta][i]]^(n-i-1),{i,n-1}]
v[n_]:=2^n Integrate[jacobian[n],{r,0,1},Sequence@@Table[{\[Theta][i],0,\[Pi]/2},{i,n-1}]]
Table[v[n],{n,10}]
Table[\[Pi]^(n/2)/(n/2)!,{n,10}]
0892狸 ◆2VB8wsVUoo
2014/08/24(日) 17:19:38.71>35 名前：KingMathematician ◆LoZDre77j4i1 ：2014/08/24(日) 17:10:09.59
> 私が実際に社会を動かす事を良しとする.
>
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
0893１３２人目の素数さん
2014/08/29(金) 23:39:08.55Table[Integrate[1, Element[Array[x, n], Ball[Table[0, {n}]]]], {n,
10}]
0894１３２人目の素数さん
2014/08/30(土) 14:11:32.46RegionMeasure[ImplicitRegion[(c1^2+c2^2)^(-1/2)>=(c1^2/a1^2+c2^2/a2^2)^(-1/2)+1,{{a1,0,1},{a2,0,1},{c1,-1,1},{c2,-1,1}}]]
Integrate[1,Element[{a1,a2,c1,c2},ImplicitRegion[(c1^2+c2^2)^(-1/2)>=(c1^2/a1^2+c2^2/a2^2)^(-1/2)+1,{{a1,0,1},{a2,0,1},{c1,-1,1},{c2,-1,1}}]]]
orz
■ このスレッドは過去ログ倉庫に格納されています