X



トップページ物理
60コメント26KB
Mathematica In[1]
0001ご冗談でしょう?名無しさん垢版2011/11/11(金) 09:08:44.57ID:CDL4iYWQ
              /: : : : : : : : : : : : : : : : : : : : : : : \
             / : : : : : : : : : : : : : : : : : : : : : : : : : : : : .
           ' : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :.',
          ': : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :.
           : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : ',
          {: : : : : :|: : : : : : : : :/: : : : : : : : : : :ヽ: :} : : : 八
          |: : : : : :|: : : : : : : : 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:|
     {,////////////\  :: .::: ∧:: |  \人   :..:/\:: \ノ
0010ご冗談でしょう?名無しさん垢版2011/11/15(火) 20:33:17.18ID:Gj21l1FN
>>9
もちっと具体的に
0011ご冗談でしょう?名無しさん垢版2011/11/16(水) 14:19:40.44ID:NEpfedOv
10
Γ(s)=lim n→∞ n!n^s/s(s+1)・・・(s+n)
これなんですが、これをグラフィック化したいのが本音です。
Γ関数は4次元での図になると思うんですが、一次元方面から視点を当てて、
Γ関数がどのように動いているのかを表現したいです!
長々と申し訳ありません・・・。お願いいたします。
0013ご冗談でしょう?名無しさん垢版2011/11/17(木) 23:02:46.11ID:FUno0ZiE
電波テロ戦争ですエンジニアさん参加を願います

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

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

こんな形で表してて、どこをどうすればいいのかわかりません。
勉強している身であまり知識がないのでどうすればいいかのご指導お願いいたします。
0016ご冗談でしょう?名無しさん垢版2011/11/29(火) 19:40:17.70ID:???
>>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]

とすれば、なんかそれらしいのがでてくるけど?
0018ご冗談でしょう?名無しさん垢版2011/12/02(金) 00:01:02.16ID:Rcph9YoF
>>16.17
ありがとうございます。
全然知識がないものですいません・・・。
ちょっとやってみてまた書き込んでみますね。
Γ関数のガウスの公式をつかっての変換の式で、本当は3のところがnなんです…。
ですが、nでやると値が大きすぎるためn=3でやってみました。
なんかすいません・・・。
0019ご冗談でしょう?名無しさん垢版2011/12/02(金) 09:32:31.40ID:???
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を別に用意するからさ。
0020ご冗談でしょう?名無しさん垢版2011/12/02(金) 09:39:49.59ID:???
あと、
Function[s,....
というのは、純関数の定義だよね。

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

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

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

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

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

Gam[]よりは、gam[]の方が吉かもしれん。
Mathematicaの定義とかぶらなきゃどうということはないが。
0022ご冗談でしょう?名無しさん垢版2011/12/06(火) 22:11:58.63ID: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
0023ご冗談でしょう?名無しさん垢版2011/12/14(水) 15:32:59.32ID:5rU3kNFu
mathematicaで、カラビヤウ多様体の描き方を
知ってる方がいらしたら、お教え下さい。
0026ご冗談でしょう?名無しさん 垢版2011/12/14(水) 18:43:50.41ID:5rU3kNFu
>>25
すみません。カラビヤウを描くためのコードを
知りたいのですが、ご存知ないでしょうか?
0027ご冗談でしょう?名無しさん垢版2011/12/14(水) 21:21:05.47ID:???
>>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
0028ご冗談でしょう?名無しさん垢版2011/12/19(月) 01:39:32.17ID:???
最近Mathematicaを触り始めて凄い楽しいのですが、躓いてしまい力を貸してください
Plot3Dは平面を書くとき、z=ax+by+dといった式に直す必要があると思います
この時(x,y,z)=(1,0,10),(0,1,10)を通るxy平面に対して垂直な面は書けるのでしょうか
zを固定してx、yを動かす方法がわからない状態です。よろしくお願いします
0029ご冗談でしょう?名無しさん垢版2011/12/19(月) 10:58:23.83ID:???
>>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の式を立てたらあかんわな。
0030ご冗談でしょう?名無しさん垢版2011/12/19(月) 14:45:20.31ID:UI7+H48W
>>29

中学生はほっとけよw
0031ご冗談でしょう?名無しさん垢版2012/01/04(水) 10:38:18.08ID: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/
0032ご冗談でしょう?名無しさん垢版2012/01/04(水) 12:58:58.52ID:4DvxlNDC
コマンドが改定されたからだよ

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

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

まだ院に入ってないなら、悪いことは言わない。
お前のようなバカはさっさと働きなさい。
0037ご冗談でしょう?名無しさん垢版2012/01/04(水) 19:02:57.60ID: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
0038ご冗談でしょう?名無しさん垢版2012/01/04(水) 19:46:15.66ID:???
>>31
compton2.m
をダウンロードして、mathematicaで開いて実行して、
M=…
M*cj[M]
%//ReleaseHold…

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

ただ、m_eとかが、TamarA`m_eみたいになってしまう。
この変なのをどうすればいいのか俺には分からん、、。
0039ご冗談でしょう?名無しさん垢版2012/01/04(水) 21:53:24.05ID: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では成功しておりません。理由は不明です。
以上、何とかよろしくお願い致します。

0040ご冗談でしょう?名無しさん垢版2012/01/04(水) 22:11:47.65ID: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
0042ご冗談でしょう?名無しさん垢版2012/01/15(日) 23:56:13.69ID:???
もうちょっと問題を切り分けて単純化してくれないと....
分野が違うものにとってはソースを見るだけでも面倒だし。
0043ご冗談でしょう?名無しさん垢版2013/01/20(日) 17:29:58.11ID:CLM/UQiZ
確かに
0046ご冗談でしょう?名無しさん垢版2014/03/19(水) 22:41:03.90ID:???
こんにちは、
下記HPのPDFは、電子の対消滅の計算式です。
(素粒子物理 牧二郎先生、林浩一先生著P134,135抜粋)
この計算を、下記HPのmathematica eμ計算の通り、作成したのですが、エラーが発生して、答が出ません。(1/0で無限になります。)
どこを修正すれば、良いでしょうか?


http://www.geocities.jp/dirac_equation/index10.htm
0047ご冗談でしょう?名無しさん垢版2014/03/20(木) 21:04:02.26ID:1scUXbOD
マルチするならせめてどちらかのスレに誘導しろよ
0049ご冗談でしょう?名無しさん垢版2014/03/21(金) 08:00:57.19ID:???
マルチをするつもりはないのですが、もう1つの方には、「人大杉」と表示され、レスできなくなりました。
現在、Jane Styleという使い慣れていないアプリを使用しているのですが、こちらからは、もう1つの方にレスできず、結果的にマルチになっている状況です。
0053ご冗談でしょう?名無しさん垢版2018/01/31(水) 07:24:00.43ID:co3m1tSm
物理学もおもしろいけどネットで儲かる方法とか
グーグルで検索⇒『羽山のサユレイザ』

IT2NN
0054ご冗談でしょう?名無しさん垢版2018/07/12(木) 20:38:40.20ID:1MdQRTZv
僕の知り合いの知り合いができた在宅ワーク儲かる方法
時間がある方はみてもいいかもしれません
検索してみよう『立木のボボトイテテレ』

UTX
0055ご冗談でしょう?名無しさん垢版2019/10/28(月) 11:12:29.19ID:Jl4CIHxg
福田博造は地獄へ落ちたのか
0056ご冗談でしょう?名無しさん垢版2020/11/27(金) 09:21:25.12ID:PTtk/056
福田博造は地獄へ落ちただろうな
0058ご冗談でしょう?名無しさん垢版2024/03/29(金) 00:15:17.43ID:y9IFEfKV
スノ1人1つスイーツ映画やるのは公休扱いだろと担任に詰め寄り担任じゃはなしにならなかったり裁判まで入力させてるからな
0059ご冗談でしょう?名無しさん垢版2024/03/29(金) 00:42:15.07ID:IYpu4ORF
全てをかけ合う系になったよ
だから全力で舵切れたな
0060ご冗談でしょう?名無しさん垢版2024/03/29(金) 00:50:38.85ID:oa+Hy/Po
なぜ買い向かうのかな
あかんな
レスを投稿する


ニューススポーツなんでも実況