〓 Mathematica 伍 〓
レス数が950を超えています。1000を超えると書き込みができなくなります。
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/
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
0895１３２人目の素数さん
2014/09/06(土) 07:49:39.29とか
689（x^3+y^7+xy+xy^4+125）
の式から、係数「215」、「689」を抽出するには、どうしたら良いでしょうか？
0896１３２人目の素数さん
2014/09/06(土) 18:16:46.060897１３２人目の素数さん
2014/09/07(日) 20:52:42.26やってみます。
0898１３２人目の素数さん
2014/09/07(日) 21:34:21.14knife.mは最近のバージョンでは使えないんだろ
PlotRangeだと斜めとか切れないし
0899１３２人目の素数さん
2014/09/08(月) 00:34:23.87著者の人は相変わらずMathematica使ってるみたいだけど
興味が折り紙に向かってる模様
http://www.jssst.or.jp/award/detail/id=94
せっかく作ったパッケージなのになんで公開やめちゃったんだろう
新しいバージョンで動かなくなってメンテやめちゃったのかな
0900１３２人目の素数さん
2014/09/14(日) 15:11:45.590901１３２人目の素数さん
2014/09/14(日) 22:00:53.730902１３２人目の素数さん
2014/09/16(火) 18:26:30.32「設問チラ見」で問題が見れる
この問題をMathematicaで解いて使用言語「Wolfram言語」として送ったった。
届いた評価フィードバックがこちら。
【解答内容】
最適解でした。
【コメント】
詳細な説明ありがとうございます。
解の全列挙も良いアピールですね。Wolfram言語の需要がどの程度あるかは予測できかねますが、良いオファーがあることをお祈り申し上げます。
【参考評価値（1～5点）】
(1)提出フォーマット 5.0
(2)入出力 5.0
(3)可読性 5.0
(4)汎用性・拡張性 5.0
(5)判断力・思考力 5.0
総合評価 5.0
10行ぐらいのコードだったんだが、コードの内容には一切触れず。
なんかよくわからんけど合ってるよ、おめでとう、ぐらいの感じだな。
Mathematicaの実行環境がないということで実行すらしてないだろう。
Wolfram Programming Cloudの存在教えてやりゃーよかったかな。
0903１３２人目の素数さん
2014/09/17(水) 18:14:38.85挑戦済み0人になってるけど。
wolfram langageだとカウントもされんのか。
0904１３２人目の素数さん
2014/09/17(水) 22:23:23.03>立体を切断するにはどうすればいいんだ?
>
>knife.mは最近のバージョンでは使えないんだろ
>PlotRangeだと斜めとか切れないし
Ver.10で立体の共通部分とか
できるようになったからやってみそ
0905１３２人目の素数さん
2014/09/17(水) 22:37:28.26＞wolfram langageだとカウントもされんのか。
ワロタｗ
なんか受付期間を変えて何度もやってるみたい
https://codeiq.jp/ace/code_teller/
https://codeiq.jp/ace/code_teller/q986
この問題は言語不問だからWolfram言語でもOK
もちろん言語指定の問題もある
https://codeiq.jp/ace/code_teller/q1076
この問題をC#で解いたらオファーが3個くらいきた
0906１３２人目の素数さん
2014/09/17(水) 22:42:07.67立体の共通部分はまだダメだ
平面だけでした
0907１３２人目の素数さん
2014/09/25(木) 21:03:59.71今は一般人でも5000円もあれば使える時代に何を言ってるのやらｗ
ただしハードが遅いけどね
0908１３２人目の素数さん
2014/09/27(土) 16:08:52.78しかも2社
プログラミング言語総選挙にはC#で応募しといた
https://codeiq.jp/wp/vote/
0909１３２人目の素数さん
2014/10/04(土) 08:03:17.17素粒子論の専門家でmathematicaを使用している人には、どのような先生がいるのでしょうか？
0910１３２人目の素数さん
2014/10/04(土) 13:23:50.720911１３２人目の素数さん
2014/10/04(土) 13:34:58.54有難う御座います。
更に、ご存じないでしょうか？
0912age
2014/10/04(土) 14:25:46.190913１３２人目の素数さん
2014/10/04(土) 14:31:12.170914１３２人目の素数さん
2014/10/22(水) 11:24:33.260915１３２人目の素数さん
2014/10/22(水) 11:35:29.20定積分の計算じゃないの
0916１３２人目の素数さん
2014/10/22(水) 13:41:46.730917１３２人目の素数さん
2014/10/23(木) 06:23:29.80ArcLength
0918１３２人目の素数さん
2014/10/27(月) 21:22:49.19閉曲線のパラメータ表示はわかってるとします。
グラフみればなんとなくどこかわかるけど、グラフが複雑になるとわけわからなくて困ってます。
0919１３２人目の素数さん
2014/11/09(日) 08:42:42.46添付は、電子・ミューオン散乱の計算です。
答えは間違っています。どこを修正すれば、良いでしょうか？
http://www.geocities.jp/dirac_equation/index1.html
0920１３２人目の素数さん
2014/11/09(日) 22:43:38.91y
0921１３２人目の素数さん
2014/11/10(月) 19:56:11.98>>920
ほんとですね。確かに、ｙの部分が、間違っておりました。
更に、電子陽電子の散乱を計算しました。
答えが全く合いません。
原因は、マンデルスタム変数の部分が、コンプトン散乱で使用していた
ものと異なるからだと思っております。
コンプトン散乱の場合は、電子はｍ、光子はｍ＝０
電子陽電子の散乱の場合は、両者ともｍとして計算するためです。
添付にコンプトン散乱の場合のマンデルスタム変数の計算をつけております。
s+t+u=2m^2
が導出出来ません。どこが間違っているのでしょうか？
http://www.geocities.jp/dirac_equation/index1.html
0922１３２人目の素数さん
2014/11/10(月) 20:34:06.13Mathematica文法で {a,b,c,d}^2 はどう計算されるか知ってる？
0923１３２人目の素数さん
2014/11/10(月) 22:35:00.45>>922
{a,b,c,d}.{a,b,c,d}にすれば、良いわけですね。
全く、考えていなかったです。
電子・ミューオン散乱の計算は、相対論的極限で,m=M＝０のときは
答えが合っていますが、そうでない時は、導出できません。
電子、陽電子散乱も、陽電子の質量が考慮されておりませんので
正しい答えが導出出来ません。（多分）
要するに、下記の変換に、ミューオンのＭ（または陽電子のｍ）を反映させて
今度は、s+t+u=4m^2等の関係になるようにすることが必要です。
でも、そのＭ（またはｍ）を、（同じ場所のはずですが）、どこに追加すれば良いかは謎です。
直感的には、対称性からして
p0→(s+m^2+M^2)/(2Sqrt[s]),
p3→(s-m^2+M^2)/(2Sqrt[s]),
t→m^2+M^2-s-u
ではないかと、思っております。
0924１３２人目の素数さん
2014/11/10(月) 22:45:43.70＞ではないかと、思っております。
計算したら、全くはずれてました。(泣)
0925１３２人目の素数さん
2014/11/10(月) 23:16:18.52>>740
0926１３２人目の素数さん
2014/11/11(火) 21:25:34.51>>740
承知致しました。
下記に
１．コンプトン散乱のマンデルスタム変数を満たす確認計算
２．電子・陽電子散乱のマンデルスタム変数を満たす確認計算
３．電子・陽電子散乱の計算
を入れております。
直感的は正しいかもしれません。（一部訂正）
p0→(s+m^2+M^2)/(2Sqrt[s]),
p3→(s-m^2-M^2)/(2Sqrt[s]),
t→２＊（m^2+M^2）-s-u
なぜか、３の計算が合いません。なぜでしょうか？
http://www.geocities.jp/dirac_equation/index1.html
0927１３２人目の素数さん
2014/11/11(火) 21:40:48.10>p0→(s+m^2+M^2)/(2Sqrt[s]),
>p3→(s-m^2-M^2)/(2Sqrt[s]),
この粒子の質量はどうなる？
0928１３２人目の素数さん
2014/11/11(火) 23:42:50.79>j0 -> p3, j1 -> -p3*Sqrt[1 - z^2], j2 -> 0, j3 -> -p3*z
この粒子の質量は？
0929１３２人目の素数さん
2014/11/12(水) 00:48:11.77お返事有難う御座います。
すいません。基本的なことを間違っていました。
この計算は、電子・陽電子ではなく、電子と電子の散乱でした。
＞この粒子の質量はどうなる？
全体的には４元運動量として、保存されます。質量だけに注目すると、ローレンツ変換により、速度に応じて、変化すると思います。
数値では解らないです。
mathematicaの20141112の「電子・電子散乱のマンデルスタム変数の確認 その１」の計算は間違っているでしょうか？いいところまで、満たしております。
１つだけ、“０”になるべきところ、“２ｍ＾２”になります。８０点です。
「電子・電子散乱のマンデルスタム変数の確認 その２」は、ｐとｋ（またqとｋ）で、質量を分担させましたが、結果は、全く、マンデルスタム変数の条件を満たさないです。
０点です。
どこから、考えれば良いでしょうか？
0930１３２人目の素数さん
2014/11/12(水) 00:52:39.68>どこから、考えれば良いでしょうか？
特殊相対論を勉強するところから。
0931１３２人目の素数さん
2014/11/12(水) 19:42:48.56自力も大切ですが、すでに聞きまくりですから(^^;)
0932１３２人目の素数さん
2014/11/12(水) 22:10:32.91お返事有難う御座います。
P0^2-p1^2-p2^2-p^2=m^2
になるということですか。
>p0→(s+m^2+M^2)/(2Sqrt[s]),
>p3→(s-m^2-M^2)/(2Sqrt[s]),
この場合、質量がm^2+M^2になって、間違い。
すると
ｍ→ｍ＋Ｍ
に変更すればよろしいでしょうか？
0933１３２人目の素数さん
2014/11/12(水) 22:49:48.89>>857
0934１３２人目の素数さん
2014/11/13(木) 23:11:44.57>>931
日本語文献は、ほぼ完璧に探しましたが、英文は探しませんでした。
英文の検索しました、まだ完全に探しておりませんが、ほしい式は見つからないです。
>>933
そうですか。
下記ＨＰの本の式（８２．２）は正しいでしょうか？
特に２番目の式が、怪しいのでは？と思っています。
http://www.geocities.jp/dirac_equation/index1.html
0935１３２人目の素数さん
2014/11/13(木) 23:44:41.29間違ってるとしたら、どこがどう間違っていて、正しくはどうだと思うのか書いてみ
0936１３２人目の素数さん
2014/11/15(土) 14:00:23.02>>931
探しましたが、ないです。
下記の最後のプログラムが近いですが、今の計算とは少し違うので、使えないです。
http://www.phy.pmf.unizg.hr/~kkumer/articles/feynman_for_beginners.pdf#search='mandelstam+variables+++mathematica'
もし、参考になるような式があれば、教えて下さい。
>>935
http://en.wikipedia.org/wiki/Mandelstam_variables
を見ると、間違っていないことがわかりました。
ｐの電子とｋの電子が衝突して、ｑとjの状態になって、飛び去る。
ｐのp0^2-p1^2-p2^2-p3^3を計算すると“ｍ＾２”になり
ｋのそれは、”0”
qのそれは、”m＾２”
ｊのそれは、“０”
になります。
ｋが０になるのは、光子質量が０であるからであり、ｋが“ｍ＾２”になるように計算すれば
良いのでしょうか？
いろいろと試しましたが、わからないです。かなり、コンプトン散乱の場合と、変更があるのではないでしょうか？
http://www.geocities.jp/dirac_equation/index1.html
0937１３２人目の素数さん
2014/11/15(土) 20:04:20.80http://www.geocities.jp/dirac_equation/index1.html
の
20141115の計算 その１
で正しい否かを教えて下さい。
0938１３２人目の素数さん
2014/11/15(土) 22:34:19.17http://www.geocities.jp/dirac_equation/index1.html
20041115 電子・電子散乱
どこに、問題があるでしょうか？
0939１３２人目の素数さん
2014/11/16(日) 08:29:42.99答えが合わないはずです。
0940１３２人目の素数さん
2014/11/16(日) 10:57:58.70１、電子・電子散乱のマンデルスタム変数を使った変換式が間違っている
（形からして、重心系の変換式のはず。82.8を代入したら一致するはず）
２、計算結果は正しいが、s,t,uの３種類の変数が混在しているため、見かけ上
異なるように見える
３、その他の箇所で、計算間違いをしている
等々、どこが悪いのか？その原因が解らないので、対策が立てられない状況です。
0941１３２人目の素数さん
2014/11/16(日) 23:27:19.42>どこに、問題があるでしょうか？
sl[j]の定義式より後は全部問題あり。「何一つ正しい式を書いてない」レベル。
0942１３２人目の素数さん
2014/11/17(月) 21:00:48.13承知しました。計算をやり直してみます。
ｐ1,p2,p3,p4とｐ、ｑ、ｋ、ｊが、対応してないかもしれません。
0943１３２人目の素数さん
2014/11/24(月) 10:03:39.52Compton散乱を計算し直しております。
下記ＨＰのCompton散乱は、ＨＰのHeitlerの本を元に計算してます。
プログラムComptonは、式（３９）に式（４）を入れて、式（４０）を導いてます。
プログラムCompton R2は、行列を一部変更したのですが、式（４０）に該当する式に、ｍが残ってしまい、グラフ化しません。グラフ化することは可能でしょうか？
変数、k0,k,u,γに対して、式（４）とγ=k0/uの関係式しかないので、計算は無理でしょうか？
http://www.geocities.jp/dirac_equation/index1.html
0944１３２人目の素数さん
2014/11/24(月) 19:05:30.29ガンマ行列の性質も知らん人にはQEDの計算は無理
0945１３２人目の素数さん
2014/11/25(火) 20:30:33.71すいません。説明不足です。
プログラムCompton R2のガンマ行列もどきは、ガンマ行列になっていません。
そもそも、プログラムCompton R2は、間違っていますが、このような場合、最後に“ｍ”が残ると、変数γだけの式にはならないでしょうか？
変数、k0,k,u,γに対して、式（４）とγ=k0/u（γ=k0/ｍ）の関係式しかないので、計算は無理でしょうか？
0946１３２人目の素数さん
2014/11/26(水) 16:53:01.39最後の式をよく読めばどうしてmが残ってるか分かる
0947１３２人目の素数さん
2014/11/26(水) 22:28:18.93＞最後の式をよく読めばどうしてmが残ってるか分かる
ご指摘有難う御座います。確認します。（現状、よくわからないです）
＞特殊相対論を勉強するところから。
今、酔っています。現状で、仲間（賢い人一杯いると思います）に「特殊相対論」って一言、言ったら、「病院で見てもらえ！」って（目で）言われます（ました）。
でも、ワイスコップの人の心を埋めるのは、「科学」「宗教」「芸術」であるいう言葉を信じます。（確信しています。）実は、数歩、進んでいる。
あくまでも、物理「科学」大好き、（あほですが）よろしくお願いいたします。
0948１３２人目の素数さん
2014/11/27(木) 12:22:31.02すいません。最後の式とは、ハイトラーの本の式（４０）のことでしょうか？
それとも（４３）でしょうか？
0949１３２人目の素数さん
2014/11/27(木) 16:37:23.94貴方が自分で出力して
>最後に“ｍ”が残る
と言ってる式
0950１３２人目の素数さん
2014/11/27(木) 22:21:11.18見て考えたのですが、わかりません。
ｍを消すことは可能でしょうか？
0951１３２人目の素数さん
2014/11/27(木) 23:07:47.44FullSimplify[(m^2)^(3/2)/m^3] はどうなるかをを考えてみ
0952１３２人目の素数さん
2014/11/28(金) 07:17:52.06有難う御座います。考えます。
0953１３２人目の素数さん
2014/11/28(金) 21:41:47.33＞FullSimplify[(m^2)^(3/2)/m^3] はどうなるかをを考えてみ
ボキャブラリーというか、ExpnadAllで展開させることしか、思い浮かばないです。
まず、このFullSimplify[(m^2)^(3/2)/m^3]を、どうしたら、“１”になるのか？わからないです。
±があるので、正確には、これ以上どうしょうも、ないのでしょうか？
（ルートｍ＾２の解は、＋ｍと-ｍの２つが考えられる。）
このコンプトン散乱の場合、±の-は、一応、無視したいです。
その場合でも、どうしたら、良いのか？わからないです。
0954１３２人目の素数さん
2014/11/28(金) 21:44:09.26√ｍ＾２→ｍ に、力ずくで、置き換えれば良いかもしれませんね。
レス数が950を超えています。1000を超えると書き込みができなくなります。