トップページ⇒math
1001コメント517KB

〓 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/
0757１３２人目の素数さん2014/04/28(月) 17:47:11.44
　　　　　　　　　　__ﾉ)-'´￣￣`ｰ- ､_
　　　　　　　 ,　'´　　_. -‐'''"二ニﾆ=-`ヽ、
　　　　　　／　　 ／:::::; -‐''"　　　　　 　 `ｰノ
　　　　　/　　　／:::::／　　　　　　　　　　　＼
　　　　 /　 　 /::::::/　　　　　　　　　 |　|　|　 |
　　　　 |　　　|:::::/　/　　　　　|　　|　|　|　|　 |
　 　 　 |　　　|::/　/　/　|　　| ||　 |　| ,ﾊ .| ,ﾊ|
　 　 　 |　　　|/　/　/　/|　,ﾊﾉ|　/|ﾉﾚ,ﾆ|ﾙ'　
　　　　 |　　　|　 |　/　/ ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/　　　
.　　　 　|　　　＼ ∠イ　 ,ｲイ|　　　 ,`-' |　　　　　　
　　　　 | 　　　　l＾,人| 　` `-' 　 　 ゝ　 |　　　　　　　　
　 　 　 | 　　　　 ` -'＼　　　 　 　ｰ' 　人　　　　　　　　　　　私は死なないわよ。
　　　　|　　　　　　　 /(l 　　 　＿＿／　 ヽ、　　　　　　　　　　　 でも最近一寸太ったかしら。
　　 　 |　　　　　　　（:::::`‐-､__ 　|::::`､ 　 　 ﾋﾆﾆヽ、　　　　　　 　　　　Windows ver.10 で　　　　
　 　　|　　　　　　／ `‐-､::::::::::`‐-､::::＼　　 /,ﾆﾆ､＼　　　　　　　　　　　　元の痩せた姿にしてよね。
　　　|　　　　　　|::::::::::::::::::|` -､:::::::,ﾍ￣|'､　 ﾋﾆ二、　＼　　　　　　　　　　 　　　
.　　 |　　　　　　/::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ 　 |　'､　　 ＼
　　 |　　　　　 /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ　　ヽ、　 |
　　|　　　　　 |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ 　/:＼__/‐､
　　|　　　　　 |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::|　く:::::::::::::￣|
　　 |　　　　 /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
　　　|　　　　|／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
　　　 |　　　/:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0758１３２人目の素数さん2014/05/02(金) 18:41:01.62
maximaのポインタ的な変数操作
v='a
v::1
をMathematicaでできるでしょうか？
aが未定義のうちは
v=a、v=`a、v=Symbol["a"]
など(いずれか)として
Set[Evaluate[v],1]
でaに値が入りますが
一度aが定義されてしまうと
1=1
になってしまうためうまく行きません
0759１３２人目の素数さん2014/05/02(金) 20:41:06.13
>>758
maximaを知らないので何とも。どういう動作をさせたいのでしょうか？
0760１３２人目の素数さん2014/05/02(金) 20:48:56.76
一時的なわり当ては

/.を使っておこないます
例
x^2+2*x+1/.x->1
こたえ4
0761１３２人目の素数さん2014/05/02(金) 21:52:00.16
>>759
あまり意味のある例ではないと思いますが
Do[Set[Evaluate[v],{}],{v,{`a,`b,`c}}]
(記述の雰囲気だけです。動きません。)
と書いたら
{a,b,c}=Table[{},{3}]
となるような感じで

maximaの場合
for v in ['a,'b,'c] do (v::[]);
と書くと
[a,b,c]:makelist([],3);
になります

(758の2行目は v:'a でした)
0762１３２人目の素数さん2014/05/02(金) 22:12:56.40
>>761
>Do[Set[Evaluate[v],{}],{v,{`a,`b,`c}}]
をやったら
{a,b,c} は {{},{},{}} を返すんじゃないの？
0763１３２人目の素数さん2014/05/02(金) 23:23:17.10
/. を使う案件
0764１３２人目の素数さん2014/05/02(金) 23:39:02.74
>>762
それが動くように見えるのは、
最初に書いた通りa,b,cが未定義の状態だけです

そのあとで
Do[Set[Evaluate[v],{1}],{v,{`a,`b,`c}}]
とやると失敗します
{{},{},{}}={{1},{1},{1}}
になるからです
0765１３２人目の素数さん2014/05/04(日) 23:13:10.13
田中ァ、見てるかぁ？ｗｗｗ
0766１３２人目の素数さん2014/05/12(月) 23:33:12.32
ネタねーのか？
0767１３２人目の素数さん2014/05/13(火) 19:02:25.85
もうすぐ
　　10
が来るぞ
0768１３２人目の素数さん2014/05/24(土) 21:18:49.41
Plot[a*x^2-10b*x,{x,0,1}]みたいに
a,bは定数として、xだけが変数でグラフを書きたいんだけど
どうしたらいい？
上記のだとグラフが一切表示されないんだけど
0769１３２人目の素数さん2014/05/24(土) 21:58:55.70
a, b を定義してないからだろ
0770１３２人目の素数さん2014/05/24(土) 22:02:28.27
>>768
a, bを引数とする関数にすればいいんじゃないの？
それをg(a,b)として、a,bを外から与えれば
グラフがかけるという風にすればいいでしょ。
0771１３２人目の素数さん2014/05/24(土) 22:30:07.77
>>770
？
プログラミングにあまり詳しくないから
使う関数名と簡単な例を用いて教えてくれ
0772１３２人目の素数さん2014/05/24(土) 22:44:19.73
自分でググれ
0773１３２人目の素数さん2014/05/25(日) 15:28:38.23
Manipulate[a, {a, {ggr, ks}}]
0774１３２人目の素数さん2014/05/26(月) 10:16:46.34
>>773
Manipulateが最適解と思うが…

768は釣りだから相手にするな。
0775１３２人目の素数さん2014/05/26(月) 20:22:22.44
>>774
釣りじゃねえよハゲ
0776１３２人目の素数さん2014/05/26(月) 20:33:36.61
数独・ナンプレ解法プログラム

report[board_] := Sow[board]

solutionQ[board_] := And[Not[MemberQ[board, 0, {2}]], report[board]; True]

candidates[board_, position_] := Complement[Range[9],
board[[position[[1]], ;;]],
board[[;; , position[[2]]]],
Flatten[Take[board, Sequence @@ (# - {2, 0} &amp; /@ (3 Ceiling[position/3]))]]]

children[board_] := With[{position = First[Position[board, 0]]}, ReplacePart[board, position -> #] &amp; /@ candidates[board, position]]

descendantSolvableQ[board_] := Or @@ (solvableQ /@ children[board])

solvableQ[board_] := Or[solutionQ[board], descendantSolvableQ[board]]
0777１３２人目の素数さん2014/05/26(月) 20:37:28.97
サンプル

board = {
{1, 0, 0, 0, 0, 7, 0, 9, 0},
{0, 3, 0, 0, 2, 0, 0, 0, 8},
{0, 0, 9, 6, 0, 0, 5, 0, 0},
{0, 0, 5, 3, 0, 0, 9, 0, 0},
{0, 1, 0, 0, 8, 0, 0, 0, 2},
{6, 0, 0, 0, 0, 4, 0, 0, 0},
{3, 0, 0, 0, 0, 0, 0, 1, 0},
{0, 4, 0, 0, 0, 0, 0, 0, 7},
{0, 0, 7, 0, 0, 0, 3, 0, 0}};

AbsoluteTiming[output = Reap[solvableQ[board]]]

{1.070002, {True, {{{{1, 6, 2, 8, 5, 7, 4, 9, 3}, {5, 3, 4, 1, 2, 9,
6, 7, 8}, {7, 8, 9, 6, 4, 3, 5, 2, 1}, {4, 7, 5, 3, 1, 2, 9, 8,
6}, {9, 1, 3, 5, 8, 6, 7, 4, 2}, {6, 2, 8, 7, 9, 4, 1, 3,
5}, {3, 5, 6, 4, 7, 8, 2, 1, 9}, {2, 4, 1, 9, 3, 5, 8, 6,
7}, {8, 9, 7, 2, 6, 1, 3, 5, 4}}}}}}

TableForm[output[[2]]]

（ｒｙ
0778１３２人目の素数さん2014/05/26(月) 20:40:19.99
ちなみにboardの左上の1を0に変えると・・・
0779１３２人目の素数さん2014/05/28(水) 03:43:31.02
初心者なんだけど
まずPECE法が調べてもさっぱりなんだけども
予測子修正子法（PECE）でP(q,Dq)の値を与えた時
(Dq,DDq)を返すにはどうしたらいいかな？
0780１３２人目の素数さん2014/05/28(水) 19:05:05.86
ませまてぃかは
ドラえもんではありません
0781狸 ◆2VB8wsVUoo 2014/05/29(木) 13:04:53.43
狸

>論理性が欠如していようがなんだろうが、芳雄は『後世に語り継がれる秀逸な結果』を残しただろうが
>なのに何が科学者の敵だ
>芳雄は過程はどうであれ結果を残した、お前は科学者をなめるんじゃない
>お前は芳雄を妨害して芳雄の研究成果に悪い影響を与えている、お前こそ研究者の敵だろうが
>今からでもいいから素直になって芳雄に謝ってこい、それぐらいはできるだろうが
>
0782１３２人目の素数さん2014/05/29(木) 17:01:52.44
2^100000000を計算しPDFファイルで書き出したら
6627ページになりびっくりしますた
0783１３２人目の素数さん2014/05/29(木) 18:49:49.56
9^9^9^9
をやってみろ
0784狸 ◆2VB8wsVUoo 2014/05/29(木) 22:37:26.33
狸

>論理性が欠如していようがなんだろうが、芳雄は『後世に語り継がれる秀逸な結果』を残しただろうが
>なのに何が科学者の敵だ
>芳雄は過程はどうであれ結果を残した、お前は科学者をなめるんじゃない
>お前は芳雄を妨害して芳雄の研究成果に悪い影響を与えている、お前こそ研究者の敵だろうが
>今からでもいいから素直になって芳雄に謝ってこい、それぐらいはできるだろうが
>
0785１３２人目の素数さん2014/05/31(土) 15:54:10.35
こんにちは、

下記がエラーになります。
修正方法を教えて下さい。


Label[begin];
h = 0;
For[k = 1, k <= 17, k++,
T1 = k*2;

If[T1 == 8 , h = g + 231;]

If[h == g + 231 , Goto[end]];
Print[h];
Label[begin];
];
Label[end];
0786１３２人目の素数さん2014/05/31(土) 20:35:30.14
>>785
>下記がエラーになります。
Moduleの中に入れないとダメなんじゃないの？よくわからないが、、。

あと、gが定義されていないけど？
0787１３２人目の素数さん2014/06/03(火) 09:49:24.23
＞＞786
有難う御座いました。
解決しませんが、なんとか、このまま使います。
0788１３２人目の素数さん2014/06/03(火) 14:16:48.20
http://library.wolfram.com/howtos/faster/index.ja.html
>Mathematica のプログラムの実行速度を上げるにはどのようにすればよいのでしょうか．
>手続き型プログラミングを避ける
0789１３２人目の素数さん2014/06/03(火) 16:01:43.87
>Mathematica のプログラムの実行速度を上げるにはどのようにすればよいのでしょうか

Mathematicaを使わない
0790１３２人目の素数さん2014/06/04(水) 19:13:42.47
＞Mathematicaを使わない

青海だね
おれなんかのりかえて3年になる。
レベルが上がった気がする。
0791１３２人目の素数さん2014/06/07(土) 02:29:38.18
テスト
0792１３２人目の素数さん2014/06/07(土) 11:01:52.48
Wolframのセミナー講師？　なんかイラツク態度というか
言い回し、口調だよな
客商売したことあんのかね？
0793１３２人目の素数さん2014/06/07(土) 16:31:44.93
どこの会社の誰だ?
0794１３２人目の素数さん2014/06/07(土) 18:03:52.33
http://www.wolfram.com/training/instructors/

これみると日本人は5人しかいない。

高校教師：1人
日本電子計算：2人
ヒューリンクス：2人

セミナー講師って段階で日本電子計算かヒューリンクス
0795１３２人目の素数さん2014/06/07(土) 18:10:56.76
どこのというかWolframの人だよ
セミナーというかMathematica紹介のプレゼン的な内容だったから、>>794に載るような
ことをする人じゃないと思う
講師というか、いわゆるエヴァンジェリストｗなのかな
0796１３２人目の素数さん2014/06/07(土) 19:34:29.66
ttp://www.gdep.jp/column/view/13
0797１３２人目の素数さん2014/06/07(土) 19:35:11.85
日本支社で講師みたいなことする人がいたっけ?
>>794　の誰かだろ
0798１３２人目の素数さん2014/06/10(火) 01:06:33.78
課題で自作関数を1万回実行して平均実行時間を測定しろってのがあるんですけど
Timing関数で時間を図るとしてどういう風にすればいいんでしょうか？
入力値と出力値は決まってます
0799１３２人目の素数さん2014/06/10(火) 18:04:17.92
Help引け
0800１３２人目の素数さん2014/06/12(木) 18:09:41.01
10が出るぞー
今日発表かな。
0801１３２人目の素数さん2014/06/12(木) 18:38:33.03
何が新しいんだ?
0802１３２人目の素数さん2014/06/12(木) 19:48:02.06
無限undo
0803１３２人目の素数さん2014/06/12(木) 20:37:29.26
それくらいか?
0804１３２人目の素数さん2014/06/13(金) 06:02:11.28
面白そうじゃん
https://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.39
Mathematica 10はデルデル詐欺か？
0806１３２人目の素数さん2014/06/27(金) 22:39:56.52
なんで9でいまさらベクトル解析記号的テンソル扱えるようにしたんだろ？
0807１３２人目の素数さん2014/06/28(土) 02:36:08.71
＞＞なんで9でいまさらベクトル解析記号的テンソル扱えるようにしたんだろ？
自己レスだが、単位系換算サポートの上での議論でゲージ原理やら接続、テンソルやらが議論にのぼったんだろうけど
ゲージ原理が時計表示のゲージに曲解されて実装されたのがゲージのような気がする。
0808１３２人目の素数さん2014/07/07(月) 17:06:30.13
フーリエ変換つったら、俺的には F[f] = ∫f(x)e^(-ikx) dx なのだが、
Mathematica では違うのな。

逆とか気持ち悪いぜ……
0809１３２人目の素数さん2014/07/07(月) 22:20:19.41
ｘが整数のときに、Ｏｋと表示する方法を教えて下さい。
下記はどこが悪いでしょうか？


x = 10;
If[x == Integer, Print["Ok"];
0810１３２人目の素数さん2014/07/07(月) 23:00:22.58
>>809
IntegerQでは如何
0811１３２人目の素数さん2014/07/08(火) 09:16:31.32
x = 10;
If[IntegerQ[x] == True, Print["Ok"]];

で出来ました。助かりました。有難う御座いました。
0812１３２人目の素数さん2014/07/09(水) 01:13:15.38
aho
あほ
あっホウ
0813１３２人目の素数さん2014/07/12(土) 01:03:23.37
Version 10 introduces a host of new areas--such as machine
learning, 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.82
お探しのページが見付かりませんでした．

Mathematica の日本語リファレンス、ページが死んでないか？
0815１３２人目の素数さん2014/07/15(火) 19:04:10.29
V10キタコレ
0816１３２人目の素数さん2014/07/17(木) 06:59:13.84
時下ますますご健勝のこととお喜び申し上げます．平素は，
Mathematicaをご利用いただき，誠にありがとうございます．

Mathematica 10日本語版が今週リリースされる予定です
0817１３２人目の素数さん2014/07/17(木) 17:30:30.33
今週はないな
0818１３２人目の素数さん2014/07/17(木) 18:40:35.11
いやもうダウンロードしてインストールしたけど

A 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
http://www.wolfram.com/mathematica/

出たねバージョン１０

なんかおもしろいこと発見したらみんな報告してや
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.29
なんで学生版だけはえーんだよ
0823１３２人目の素数さん2014/07/18(金) 20:57:14.34
うん？会社のもアップされてたぞ。
0824１３２人目の素数さん2014/07/18(金) 21:21:30.75
なんか性能が落ちている部分があるなあ。
0825１３２人目の素数さん2014/07/18(金) 23:46:53.45
>>814
復活してるな。バージョンアップのためのメンテナンスだったのか。
0826１３２人目の素数さん2014/07/19(土) 00:09:14.36
Undo がまともに機能するようになって凄く便利
0827１３２人目の素数さん2014/07/19(土) 08:34:03.13
何で性能落ちてるんだw
https://pbs.twimg.com/media/Bs1SW5JCIAAwGWB.png:large
0828１３２人目の素数さん2014/07/19(土) 11:21:52.79
>>827

より時間がかかってしまうのは、
積分だから今までより多くの変数変換の適用を試みてると考えればまだわかる

ただ、問題は最後の結果。
解ける問題が解けなくなってしまっているのは。。。
ウルフラムに報告しろ！
0829１３２人目の素数さん2014/07/19(土) 11:48:05.18
俺はまだV10もってないんだが、
V10でこの問題解けるようになったか誰か試してみてくれる？

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.07
今走らせてるけど全然結果返ってこないな。結構待ってるんだが。
0831１３２人目の素数さん2014/07/20(日) 14:11:42.16
>>829を計算するプログラム作ってたんだけどバグ見つけた
同じ計算させてるのに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.85
>>831
V9だけど結果同じ
08338322014/07/20(日) 14:48:26.59
× V9
○ V8
0834１３２人目の素数さん2014/07/20(日) 15:09:22.91
V10だけど結果同じ
http://up3.viploader.net/omoshiro/src/vlomoshiro037682.png
0835１３２人目の素数さん2014/07/20(日) 15:31:30.47
>>832-833
確認ありがとう。
「結果同じ」っていうのは
●俺と結果が同じ、つまり、計算結果は両者で違う
のか
●計算結果が両者で同じ
のどちらの意味かな？

>>834
ワオ！正しい方に結果が落ち着いてるね
バグ報告しようかと思ったけどV10でなおってるならいいや
確認ありがとう。

V8とV10は正しく計算されるのにV9だけおかしいのかな。。。

http://up3.viploader.net/omoshiro/src/vlomoshiro037683.png
0836１３２人目の素数さん2014/07/20(日) 15:36:08.02
9.01だけど同じ答だった。
0837１３２人目の素数さん2014/07/20(日) 15:38:48.86
>>836
一回目と二回目で等しい結果が出たという意味ね。
08388322014/07/20(日) 15:44:42.55
>>834
と同じ
0839１３２人目の素数さん2014/07/20(日) 15:49:10.80
>>836-837
確認ありがとう。
俺もV9って書いてたけど正確にはV9.0.1で、
パソコンがおかしくなってんのかと思って再起動してみたけど、
結果は同じ（答えが変わる）だった。
ちなみにWindows版ね。
同じバージョンなのにふるまいが変わるとか一体どうなってんだｗ
08408322014/07/20(日) 15:51:38.14
http://i.imgur.com/kn5ODwv.jpg

いや，よくみたら>>834とは符号逆だな
08418322014/07/20(日) 15:56:15.47
http://i.imgur.com/klIyQok.jpg

9.0.1.0

全然検証してないけど見た目の結果は確かにV8と違うな
0842１３２人目の素数さん2014/07/20(日) 16:24:47.23
>>840-841
V8は間違ってる方に落ち着いてるねｗこれはたちが悪いｗ
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.59
なんか条件による(Ver8,9)

10でQuit[]してスタート->初回が異なる
キャッシュをリセット->初回の値が返る
http://imgur.com/Gxup1AK
0844１３２人目の素数さん2014/07/20(日) 17:23:14.18
>>843
Mathematicaを起動してない状態から起動して
>>831をコピペして実行するとどうなる？
答えが変わるんならV10もV9と同じ振る舞いということになるな

>>834のV10はMacっぽいから、
Macだと両者正しい値がでるのか、
もしくは>>834が他のノートブック開いててそこで計算させてた状態で
コピペ→実行だとキャッシュがらみでうまくいったのかも
08458432014/07/20(日) 17:58:18.29
>>844
10起動 -> 831コピペ&実行 -> 答えは異なる(初回プラス)

なお、
10では仮定Thread[0 < ss]に等号を入れても、8&9のように
ConditionalExpressionは返らず、共にマイナス側が返って
きます。よくわからん。。。
0846１３２人目の素数さん2014/07/20(日) 22:30:51.74
とりあえずn=2に関しては厳密解が計算でいるようになった
だけどせっかく一般的に作ったのに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.27
>>845
10でもなおってないんだ
ま、等号を含めるかどうかで体積は変わらないはずだから
そもそもそれもバグだよねｗ
V8が自信満々に間違った答えを出力してることからも
Integrateについてはあんまり信用しない方がいいね
0848１３２人目の素数さん2014/07/21(月) 02:22:57.43
>>844
> >>834のV10はMacっぽいから、

俺834だけど、起動した直後にコピペして計算結果出しただけだからキャッシュとか関係ないだろう。
0849１３２人目の素数さん2014/07/21(月) 21:28:48.19
俺の考察を「メモ」ページにまとめときました
ttp://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
> {p1^0.5 -

なんで 1/2 じゃなくて 0.5 とか使うのかね
0853１３２人目の素数さん2014/08/15(金) 16:47:49.98
>>851
>どこが間違っているのでしょうか？
自分が何を計算してるのかを全く理解せずに式だけいじろうとしてるところ
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
>>854
>>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に計算結果を追加しました。
■ このスレッドは過去ログ倉庫に格納されています