Mathematica In[1]

このエントリーをはてなブックマークに追加
1ご冗談でしょう?名無しさん
              /: : : : : : : : : : : : : : : : : : : : : : : \
             / : : : : : : : : : : : : : : : : : : : : : : : : : : : : .
           ' : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :.',
          ': : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :.
           : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : ',
          {: : : : : :|: : : : : : : : :/: : : : : : : : : : :ヽ: :} : : : 八
          |: : : : : :|: : : : : : : : i : : ;イ:|;ハ.: : : : :.:.}: ハ_ : : : ミ\
          ', : : i : :|: : : : : : : : | : / j:.|_,.}: : : : :.:.j/´」_:|: :i: :i  ヽ
          ; : :i: :._|; : : : : : : : |:/ ̄|ノ_,,リヘ: : : :圻リノ}:|: :|ハ|
          i : j: :{ i: : : : : : : :代圻~ソヽ  \/`" 小:.:.| !   
          }: ノ:八 i: : : : : : : | `¨´        /:ノ: : :|
          |/: : ::::ヽi::. : : : : : |        , /:::::: : : |
.           /: : : .::::::::::::. : : : :.:.i           /: ::::::: :.八    Mathematica Thread in Physics
         /: : : .::::::::::::::::. : : : : :i、      ´’/ : : |:: : : : :.',
         ,: : : : ::::::::::::::::::::. : : : 八.> ,_ .ィ:::::::: : : |、: : : : : ::.
        ,: {: : :.:::::::::::::::::::::::. : : : :iベ¨7气!}:::::::::: : |.:{ハ : : : : : ,
        {: |: : :::::::::::::::::::ハ::. : : : i. ∨=斗|::::::::::. :.| :i人: : : : : '
          -┴──…''"\∧:::. : :.i,/}  j!ト、:::::::. :.|: :i. \ : : :}
        / -─-ミ       \:,::::. :.i   ノ圦::\:: 人┴-ミヘ: : |
      j___         ∧:::. :.!\,j_j\)::::\: :\ ∧‘:.:|
     〃///////`ヽ \  :.   ∧:::. |\ノ ヽ    \:: \∧ i:|
     {,////////////\  :: .::: ∧:: |  \人   :..:/\:: \ノ
2ご冗談でしょう?名無しさん:2011/11/11(金) 09:08:56.94 ID:???
関連スレ
〓 Mathematica 伍 〓
http://kamome.2ch.net/test/read.cgi/math/1320969748/
3ご冗談でしょう?名無しさん:2011/11/11(金) 10:08:42.99 ID:???
Mekosujica 淫[69]
4ご冗談でしょう?名無しさん:2011/11/11(金) 19:16:10.54 ID:???
え?Mathematicaって日本人には売らないんじゃなかったの?
5ご冗談でしょう?名無しさん:2011/11/11(金) 23:12:35.82 ID:kZMeMUHH
>>4
>Mathematicaって日本人には売らないんじゃなかったの?

日本だけ特別に高額で販売させていただいております
6ご冗談でしょう?名無しさん:2011/11/11(金) 23:27:10.12 ID:???
もうね。
実験データのグラフ化とか、Mathematicaばかり。

テキストファイルから必要な箇所を切り出して、幾つものグラフに書いて
それぞれ関数フィッティングして、フィッティングで得られた値を
集めてきて、別のグラフに書く...とか
出張実験で150ぐらいのグラフをその場でつくるとか、
一回関数を組んじゃうと便利なんだもの。
7ご冗談でしょう?名無しさん:2011/11/12(土) 11:29:54.03 ID:???
お金のない人はMaxima使いましょう
8ご冗談でしょう?名無しさん:2011/11/12(土) 17:30:36.45 ID:43I4XuEY
Adobeのライセンスが変わって年貢を納めなければいけないとか

Matheamticaもそうなりそう
9ご冗談でしょう?名無しさん:2011/11/15(火) 16:09:42.81 ID:Oz/2tF71
ガンマ関数でのガウスの公式をMathematicaをつかってあらわしたいんですが・・・。
初歩的なのわからなくてすいません・・。
10ご冗談でしょう?名無しさん:2011/11/15(火) 20:33:17.18 ID:Gj21l1FN
>>9
もちっと具体的に
11ご冗談でしょう?名無しさん:2011/11/16(水) 14:19:40.44 ID:NEpfedOv
10
Γ(s)=lim n→∞ n!n^s/s(s+1)・・・(s+n)
これなんですが、これをグラフィック化したいのが本音です。
Γ関数は4次元での図になると思うんですが、一次元方面から視点を当てて、
Γ関数がどのように動いているのかを表現したいです!
長々と申し訳ありません・・・。お願いいたします。
12ご冗談でしょう?名無しさん:2011/11/16(水) 14:54:51.36 ID:???
>>11

これじゃいかんの?

参考1:Mathematica Documentation 3.2.10 特殊関数
http://reference.wolfram.com/legacy/v5/TheMathematicaBook/AdvancedMathematicsInMathematica/MathematicalFunctions/3.2.10.ja.html

参考2:EMANの物理学・統計力学・ガンマ関数
(ガンマ関数の定義、複素数への拡張)
http://homepage2.nifty.com/eman/statistic/gamma_func.html
13ご冗談でしょう?名無しさん:2011/11/17(木) 23:02:46.11 ID:FUno0ZiE
電波テロ戦争ですエンジニアさん参加を願います

公安はサリンオウム信者の子供を40歳まで社会から隔離している

オウム信者が地方で現在も潜伏している
それは新興宗教を配下としている公安の仕事だ
で、盗聴機器を開発したら霊魂が寄って呼ぶ来る
電波憑依
スピリチャル、全否定なら江原三輪氏、高橋佳子大川隆法氏は幻聴強制入院です矛盾する日本宗教と精神科
コードレス盗聴
すでに2004年国民の20%は被害<+>エンジニアさん電波戦争しかない<+>中国工作員ふざけるな<+>250〜700台数3万〜7000万円<+>医師も開発絡んだソウル魂インコピー機<+>
盗聴証拠
今年の5月に警視庁防犯課は、被害者のSDカード15分を保持した
有る、国民に出せ!!
*創価は潰せる
犯人は創刊学会幹部キタオカ1962年東北生は、二十代で2人の女性をレイプ殺害して入信した
創価本尊はこれだけで潰せる
*創価幹部は韓国工作員こうのとり学会軍団
創価会員と言えば公明党
<<<<<テロ装置<<東芝部品<<<宗教<<<同和>>>>公安>>>医師>>>魂複写>>>官憲>>>>>日本終<<<Google検索へ
14ご冗談でしょう?名無しさん:2011/11/17(木) 23:05:17.46 ID:???
うるせぇ死ねカス
15ご冗談でしょう?名無しさん:2011/11/29(火) 16:53:43.28 ID:Q3nI8gxP
12<<
ありがとうございます!
また聞きたいのですが、いまCartesianMapを使って図に表したいのですが・・・
現状
Module[{s,x,y,},
s=x+I Y;
Gam[s_]:=Function[s,3!*3^x e(iy*log3)/(s*(s+1)*(s+2)*(s+3))]]

こんな形で表してて、どこをどうすればいいのかわかりません。
勉強している身であまり知識がないのでどうすればいいかのご指導お願いいたします。
16ご冗談でしょう?名無しさん:2011/11/29(火) 19:40:17.70 ID:???
>>15
>Module[{s,x,y,},
>s=x+I Y;
>Gam[s_]:=Function[s,3!*3^x e(iy*log3)/(s*(s+1)*(s+2)*(s+3))]]

どこまでが指数なんだか分母なんだか、わかりにくいのう。
関数の意味もよく分からん。

それはともかく、俺のはver.8なんでCartesianMapでなく
ver.8のParametricPlotで書いていいなら、

Gam[s_] :=
3! 3^(Re[s] Exp[I Im[s] Log[3]])/(s (s + 1) (s + 2) (s + 3));

ParametricPlot[
Through[{Re, Im}[Gam[x + I*y]]], {x, -1, 1}, {y, -2, 2},
PlotStyle -> None]

とすれば、なんかそれらしいのがでてくるけど?
17ご冗談でしょう?名無しさん:2011/11/29(火) 19:41:53.21 ID:???
つか、e(iy*log3) とかそのままタイプして、なんか値を返すのか?
18ご冗談でしょう?名無しさん:2011/12/02(金) 00:01:02.16 ID:Rcph9YoF
>>16.17
ありがとうございます。
全然知識がないものですいません・・・。
ちょっとやってみてまた書き込んでみますね。
Γ関数のガウスの公式をつかっての変換の式で、本当は3のところがnなんです…。
ですが、nでやると値が大きすぎるためn=3でやってみました。
なんかすいません・・・。
19ご冗談でしょう?名無しさん:2011/12/02(金) 09:32:31.40 ID:???
s=x+I y;
Gam[s_]:=....

こうしちゃうとさ
1行目のsと2行目のsは無関係になるよね。

たとえばさ、
s = 1; f[s_] := s^2; f[2]
とすると、4でしょ。
もちろん、?sをすれば、s=1がかえってくる。

f[s_] := とすると、関数定義の外のグローバルなsをブロックして、
関数の定義式の中だけで通用するsを別に用意するからさ。
20ご冗談でしょう?名無しさん:2011/12/02(金) 09:39:49.59 ID:???
あと、
Function[s,....
というのは、純関数の定義だよね。

これと、gam[s_]:=...という関数の定義を混ぜる意味がないように思う。

x^2を計算する関数を定義するときには、
 普通は、f[x_]:= x^2;でいい。
それで、f[2]は4を返す。

もしも、純関数で書くなら、
Function[x, x^2][2];
だね。それで4を返す。

マニュアルの純関数の説明では、
「これまでに示した例では,関数の「名前」を明記することでこれを行ってきた.
純関数と呼ばれる式の構成形式を使うことで,関数名を宣言しなくても,
引数に適用させることができる関数を定義することができる.」
とある。
21ご冗談でしょう?名無しさん:2011/12/02(金) 09:52:48.25 ID:???
最後に、

Mathematicaのマニュアルによると、
 大文字から始まる変数名や関数名は、
 Mathematicaが既に使っていることがあるので、
 ユーザーはなるべく小文字から始まる名前をつかってね。
というようなことが書いてある。

Gam[]よりは、gam[]の方が吉かもしれん。
Mathematicaの定義とかぶらなきゃどうということはないが。
22ご冗談でしょう?名無しさん:2011/12/06(火) 22:11:58.63 ID:4iF9ktds
>> 11,15,18

と謝辞に書いとけよ。こんな感じで。

I also would like to thanks Mr. Joking No Name at 16 and 17
in the thread of the Board of Physics in 2chan
23ご冗談でしょう?名無しさん:2011/12/14(水) 15:32:59.32 ID:5rU3kNFu
mathematicaで、カラビヤウ多様体の描き方を
知ってる方がいらしたら、お教え下さい。
24ご冗談でしょう?名無しさん:2011/12/14(水) 15:55:59.58 ID:???
mekosujicaで、アワビヤラシイ多様体の描き方を
知ってる方がいらしたら、お教え下さい。
25ご冗談でしょう?名無しさん:2011/12/14(水) 17:02:04.96 ID:???
>>23

ググったら30秒くらいで見つかった。
バージョン5.2用とバージョン6用がwolframにあるね。
http://members.wolfram.com/jeffb/visualization/stringtheory.shtml
26ご冗談でしょう?名無しさん :2011/12/14(水) 18:43:50.41 ID:5rU3kNFu
>>25
すみません。カラビヤウを描くためのコードを
知りたいのですが、ご存知ないでしょうか?
27ご冗談でしょう?名無しさん:2011/12/14(水) 21:21:05.47 ID:???
>>26
>すみません。カラビヤウを描くためのコードを
>知りたいのですが、ご存知ないでしょうか?

質問の意味がよくわからんな。
だから、>25のリンク先にあるノートブックのコードはだめなの?

バージョン5.2用
http://members.wolfram.com/jeffb/visualization/notebooks/calabi-yau52.nb

バージョン6用
http://members.wolfram.com/jeffb/visualization/notebooks60/calabi-yau60.nb
28ご冗談でしょう?名無しさん:2011/12/19(月) 01:39:32.17 ID:???
最近Mathematicaを触り始めて凄い楽しいのですが、躓いてしまい力を貸してください
Plot3Dは平面を書くとき、z=ax+by+dといった式に直す必要があると思います
この時(x,y,z)=(1,0,10),(0,1,10)を通るxy平面に対して垂直な面は書けるのでしょうか
zを固定してx、yを動かす方法がわからない状態です。よろしくお願いします
29ご冗談でしょう?名無しさん:2011/12/19(月) 10:58:23.83 ID:???
>>28

おいおい、落ち着けよ。
Mathematicaというより、まずは中学か高校の数学の話だべさ。

まず、条件は、
 a) 求める図形を3次元空間内の平面S ; ax+by+cz+d=0であるとする。
 b) 平面Sは、xy平面に垂直である。
 c) {1,0,10}と{0,1,10}の2点を通る。
でよろしいか?

すると、b) は、
 b') 平面Sは、xy平面に垂直なベクトル{0,0,1}と平行な平面である。
と言い換えられる。
c) + b')より、
 d) 平面Sは、{1,0,10}と{0,1,10}と{1,0,10}+{0,0,1}の3点を通る。
と言える。

つまり、{1,0,10}と{0,1,10}と{1,0,11}が、 ax+by+cz+d=0を満たす。

Solve[{ a 1 + b 0 + c 11 + d == 0, a 1 + b 0 + c 10 + d == 0 ,
a 0 + b 1 + c 10 + d == 0}, {a, b, c, d}]
結果は、
{{b -> a, c -> 0, d -> -a}}

つまり、求める平面の式は、簡単化して、
 x + y + 0 z -1 ==0
こういう陰関数をPlotするなら、
ContourPlot3D[
x + y + 0 z - 1 == 0, {x, -1, 1}, {y, -1, 1}, {z, -1, 1},
AxesLabel -> {x, y, z}, AxesStyle -> Directive[FontSize -> 15]]
かね。ちなみに、ax+by+cz+d=0でc=0が答えなのに、c≠0の式を立てたらあかんわな。
30ご冗談でしょう?名無しさん:2011/12/19(月) 14:45:20.31 ID:UI7+H48W
>>29

中学生はほっとけよw
31ご冗談でしょう?名無しさん:2012/01/04(水) 10:38:18.08 ID:mXW+BRuz
こんにちは、
下記HPの添付5を、mathematica4.2で実行すると、添付6のように計算ができます。
http://www.geocities.jp/dirac_equation/kurukomi/kurikomi.html
しかし、下記HPからmathematica8の評価版をインストールして、計算すると、ちゃんと答えが出ません。どうしたらmathematica8の評価版でも計算出来るでしょうか?
http://www.wolfram.com/mathematica/trial/
なお、添付5は、下記から入手しました。
http://library.wolfram.com/infocenter/MathSource/839/
32ご冗談でしょう?名無しさん:2012/01/04(水) 12:58:58.52 ID:4DvxlNDC
コマンドが改定されたからだよ

と見ないでレス
33ご冗談でしょう?名無しさん:2012/01/04(水) 13:42:36.71 ID:mXW+BRuz
31です。
そうですか。1度、お試し頂ければ幸いです。
また、評価版だからかもしれませんが、mathematica8は、起動用アイコンをデスクトップに持ってきて使用すると、途中で停止します。そんなことはないでしょうか?vistaを使用しております。

34ご冗談でしょう?名無しさん:2012/01/04(水) 16:27:16.20 ID:4DvxlNDC
ショートカットをデスクトップに作るといいよ
35ご冗談でしょう?名無しさん:2012/01/04(水) 16:32:48.92 ID:P4R0322E
>>33

まだ院に入ってないなら、悪いことは言わない。
お前のようなバカはさっさと働きなさい。
36ご冗談でしょう?名無しさん:2012/01/04(水) 16:39:36.20 ID:???
坂の上のメコスジ
37ご冗談でしょう?名無しさん:2012/01/04(水) 19:02:57.60 ID:mXW+BRuz
>お前のようなバカはさっさと働きなさい。
お返事有り難うございます。
バカだから、質問しております。賢いのなら教えて下さい。秀才なら簡単なことでしょう。
本当に困っております。お願いします。
不明点は、以下です。

下記HPの添付5のミュウ粒子の崩壊計算は、mathematica4.2で実行すると、添付6のようにちゃんと計算ができます。
http://www.geocities.jp/dirac_equation/kurukomi/kurikomi.html
しかし、添付5のコンプトン散乱の計算が出来ません。どこをどのように修正すれば良いでしょうか?
なお、添付5は、下記から入手しました。
http://library.wolfram.com/infocenter/MathSource/839/

なお、コンプトン散乱の式は下記HPの添付7 式(6.103)を使っています。
http://www.geocities.jp/dirac_equation/kurukomi/kurikomi.html
38ご冗談でしょう?名無しさん:2012/01/04(水) 19:46:15.66 ID:???
>>31
compton2.m
をダウンロードして、mathematicaで開いて実行して、
M=…
M*cj[M]
%//ReleaseHold…

とその通りやっていけば大体出てくるみたいだけど?

ただ、m_eとかが、TamarA`m_eみたいになってしまう。
この変なのをどうすればいいのか俺には分からん、、。
39ご冗談でしょう?名無しさん:2012/01/04(水) 21:53:24.05 ID:mXW+BRuz
お返事有り難うございます。
実行して頂いて、本当に有り難うございます。
>ただ、m_eとかが、TamarA`m_eみたいになってしまう。
>この変なのをどうすればいいのか俺には分からん、、。
その通りです。最初、私がやった際にも、TamarAが現れて困りました。
まず、このTamarAが出ないようにする必要があります。
TamarAGuide.nb
に書かれている通りの結果になる必要があります。
まず、TamarA.mを実行して、その後、TamarAGuide.nbを、実行する必要があります。
もしくは、compton2.m の上の塊を実行して、その後、
M=…
M*cj[M]
%//ReleaseHold…
を、実行してもOKのはずです。
ミュウの崩壊だけは、必ずできるはずです。但し、私はmathematica4を使用しております。
はmathematica8では成功しておりません。理由は不明です。
以上、何とかよろしくお願い致します。

40ご冗談でしょう?名無しさん:2012/01/04(水) 22:11:47.65 ID:mXW+BRuz
下記HPの添付8にtamarAの説明
添付9に、μ粒子の崩壊、μ粒子の生成、コンプトン散乱の式、計算結果を記載しました。添付9をご覧になるとわかりますが、μ粒子の崩壊、μ粒子の生成は正しい結果が導かれていますが、コンプトン散乱は、むちゃくちゃな酷い結果です。
M1=lg[ub[p],cj[eps[n]]*I*e*gm[up[n]],I*(sl[p]+sl[k]+m)/((p+k)^2-m^2),I*e*gm[up[m]]*eps[m],u[p]]
を、どのように修正すべきでしょうか?

http://www.geocities.jp/dirac_equation/kurukomi/kurikomi.html
41ご冗談でしょう?名無しさん:2012/01/15(日) 11:29:54.22 ID:???
ご機嫌如何でしょうか?
42ご冗談でしょう?名無しさん:2012/01/15(日) 23:56:13.69 ID:???
もうちょっと問題を切り分けて単純化してくれないと....
分野が違うものにとってはソースを見るだけでも面倒だし。
43ご冗談でしょう?名無しさん:2013/01/20(日) 17:29:58.11 ID:CLM/UQiZ
確かに
44ご冗談でしょう?名無しさん:2013/04/14(日) 06:12:42.43 ID:???
で、コンプトンは片づいたのかな
45ご冗談でしょう?名無しさん:2014/03/19(水) 22:33:19.85 ID:???
以下を教えてください。
46ご冗談でしょう?名無しさん:2014/03/19(水) 22:41:03.90 ID:???
こんにちは、
下記HPのPDFは、電子の対消滅の計算式です。
(素粒子物理 牧二郎先生、林浩一先生著P134,135抜粋)
この計算を、下記HPのmathematica eμ計算の通り、作成したのですが、エラーが発生して、答が出ません。(1/0で無限になります。)
どこを修正すれば、良いでしょうか?


http://www.geocities.jp/dirac_equation/index10.htm
47ご冗談でしょう?名無しさん:2014/03/20(木) 21:04:02.26 ID:1scUXbOD
マルチするならせめてどちらかのスレに誘導しろよ
48ご冗談でしょう?名無しさん:2014/03/20(木) 23:41:58.99 ID:???
このスレの住民って数学板のスレも見てるんじゃね
逆もしかり

マルチしても意味ないんじゃないかと
49ご冗談でしょう?名無しさん:2014/03/21(金) 08:00:57.19 ID:???
マルチをするつもりはないのですが、もう1つの方には、「人大杉」と表示され、レスできなくなりました。
現在、Jane Styleという使い慣れていないアプリを使用しているのですが、こちらからは、もう1つの方にレスできず、結果的にマルチになっている状況です。
50ご冗談でしょう?名無しさん:2014/03/21(金) 08:51:31.05 ID:???
再度、mathematica 伍を、Jane Styleで探すと、見つかり、レスも出来そうです。
51ご冗談でしょう?名無しさん
>>49
最近鯖がai.2ch.netに変わった