〓 Mathematica 七 〓

1 名前：１３２人目の素数さん mailto:sage [2015/02/27(金) 10:38:04.85 ID:EbhLAWFK.net] 　　　　　　　　_....._{{ 〃 　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､ 　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､ .　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-< 　 　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）) 　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ. 　　　',　',　|Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　 ﾞ', 　　　 ヽ､', l:!Kノ}. 　 　 f:_.)iﾞi: ﾘ !　l　ル　 　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉ ﾞ ,ﾘ l ﾘ'´ 　 　　 .　　　　 ',|!!､ 　 　r‐┐ 　 ｀ ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣ 　 　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　 ＜　 Mathematica に関する話題はここに書いてね！ 　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　 　 　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿ 　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､ 　 　　/＼＼　　　　',　　 }　　 //｀ヽ 過去スレ 〓Mathematica〓 science2.2ch.net/test/read.cgi/math/985023298/ 〓 Mathematica 2 〓 science3.2ch.net/test/read.cgi/math/1078534285/ 〓 Mathematica 3 〓 science6.2ch.net/test/read.cgi/math/1121413040/ 〓 Mathematica 四 〓 kamome.2ch.net/test/read.cgi/math/1197270001/ (鯖飛びでログ消滅) mimizun.com/log/2ch/math/1197270001/ 〓 Mathematica 5 〓 kamome.2ch.net/test/read.cgi/math/1285859504/ 〓 Mathematica 伍 〓 (実質6) wc2014.2ch.net/test/read.cgi/math/1320969748/ 100 名前：１３２人目の素数さん mailto:sage [2015/11/25(水) 15:16:08.17 ID:Qdnl0bDq.net] >>98 前スレより再録 >全然理解してないって昔から指摘されてるのに、全く進歩しないねぇ >「ローレンツ不変量」とか知らんのだろう？ >特殊相対論を勉強するところから。 >デタラメやる前にちっとは勉強したら >特殊相対論の勉強を頑に拒むのはなぜ？ 101 名前：１３２人目の素数さん mailto:sage [2015/11/25(水) 17:30:57.70 ID:UERObkTV.net] お返事有難う御座います。 Simplify[y3 //. {p1 -> 0, p2 -> 0, 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, p0 -> (s + m^2)/(2 Sqrt[s]), p3 -> (s - m^2)/(2 Sqrt[s]), z -> 1 + t/(2 p3^2), t -> 2 m^2 - s - u}] この式（重心系）の意味は、下記HPのP4を見てやっと理解できました。 osksn2.hep.sci.osaka-u.ac.jp/~naga/kogi/konan-class04/ch6-qed.pdf この式(6.26a)を、自力で導出するのは、かなり難しいと思います。 同様に、（知っている実験室系の式）の意味を理解するのは、出来ますが、導出するのは困難です。 従いまして、偏極を指定した断面積を自力で導出するのは、時間の無駄ですので、取りあえず明らめます。 www.kadokawa.co.jp/sp/201312-01/ 下� 102 名前：LHPでの非偏極での計算でOKとします。 http://amonphys.web.fc2.com/amonfc.pdf [] [ここ壊れてます] 103 名前：１３２人目の素数さん mailto:sage [2015/11/27(金) 16:45:33.20 ID:K8C1Fes7.net] 応用問題が解けなきゃ理解したことにならないって高校くらいで教わらないのか 104 名前：１３２人目の素数さん mailto:sage [2015/12/06(日) 02:23:19.03 ID:9ZV2BAqG.net] 初めて書き込みます、宜しくお願いします。 三次関数の解を視覚的に表現してみたくて工夫しました。 切片を変数にしてその解を重ねてみました。 X,Y軸と、奥行きにｉ軸をとっています。 複素数の解も含めれば、どんな三次関数も解が ３個有ると視覚的に分かります。 この表現方法には何の根拠も有りませんが、 色々試してみたら面白いグラフが描けました。 どんな意味が有るのか自分でも分かりませんが、 とりあえず見てやってください。 sol[y_] = Solve[y^2 + y == x^3 - x^2, x] ans[j_, t_] := x /. sol[t][[j]] Do[gr[n] = Graphics3D[Table[Point[{Re[ans[n, t]], Im[ans[n, t]], t}], {t, -2, 0.5, 0.001}]], {n, 1, 3}] Show[gr[1], gr[2], gr[3], ViewPoint -> {-0.5, -5, 0}] 105 名前：１３２人目の素数さん mailto:sage [2015/12/06(日) 04:53:57.45 ID:/Zdpd7nb.net] なんじゃこりゃ〜 106 名前：１３２人目の素数さん mailto:sage [2015/12/06(日) 08:12:18.79 ID:9ZV2BAqG.net] 見ていただいてありがとうございます。 本来は y = ax^3 + bx^2 + cx + d の解を視覚的に見るために作りました。 a=1 b=-5 c=3 d=4 sol[t_]=Solve[a*x^3+b*x^2+c*x+d==t,x] ans[j_,t_]:=x/.sol[t][[j]] m=0.1 Do[gr[n]=Graphics3D[Table[Line[{{Re[ans[n,t]],Im[ans[n,t]],t},{Re[ans[n,t+m]],Im[ans[n,t+m]],t+m}}],{t,-10,10,m}]],{n,1,3}] Show[gr[1],gr[2],gr[3],Axes->True,AspectRatio->Automatic,PlotRange->{{-5,5},{-5,5},{-10,10}},ViewPoint->{2,-5,0}] a,b,c,dは任意の係数を指定します。 mは曲線の滑らかさです。 y軸の任意の値が含まれるx-i複素数平面に 解が必ず３個有ります(重解含む)。 書き替えれば４次、５次関数でも表現できます。 係数には複素数も代入できます。 適当にいじってみていただけたら幸いです。 107 名前：１３２人目の素数さん [2015/12/07(月) 15:12:13.62 ID:eC7hl6uO.net] 面白い曲線だな。ParametricPlot3Dを使った方が 簡単だと思うよ。好みは人それぞれだけどね。 x /. Solve[y^2 + y == x^3 - x^2, x]; Map[{Re[#], Im[#], y} & , %]; ParametricPlot3D[%, {y, -2, 1}, AxesLabel -> {"Re[x]", "Im[x]", "y"}, PlotLabel -> "y^2+y==x^3-x^2"] 108 名前：１３２人目の素数さん mailto:sage [2015/12/07(月) 20:01:25.33 ID:ONtsWkUw.net] ご教授ありがとうございます。 その通りですね。 ここで使用している数式は12/4(金)にNHKＥテレで放送された 数学ミステリー白熱教室で紹介されたものです。 ある数論の問題 y^2 + y == x^3 - x^2 ３次方程式の解で、 素数Pを法とする場合はある調和解析の関数となり、 複素数の場合は幾何学のトーラスとなるという話でした。 ラングランズ・プログラムというものらしいですが、 私には知識不足で理解できませんでした。 とりあえず、グラフの中央にある歪んだ円に大変興味を持ちました。 この辺が何かを意味しているような気がしますが、 私が考えられる次元はこの程度です。 どなたか知識のある方が何かを発見する手掛かりとなれたら本望です。 109 名前：１３２人目の素数さん [2015/12/07(月) 20:25:40.85 ID:Tq8BIhkH.net] 潜望鏡深度まで浮上！ 110 名前：１３２人目の素数さん mailto:sage [2015/12/08(火) 00:14:04.66 ID:KipLHGUT.net] 俺もグラフィックスオブジェクトをShowする方が好きだな Maximaに書き換えるとき見た目を合わせ易いし 111 名前：１３２人目の素数さん mailto:sage [2015/12/08(火) 03:33:17.43 ID:XscFjqkO.net] ただの数字遊びですが、これも面白かったので載せます。 たまねぎのような構造 sol[y_] = Solve[y^3 + y^2 == x^4, x] ans[j_, t_] := x /. sol[t][[j]] Do[gr[n] = Graphics3D[Table[Point[{Re[ans[n, t]], Im[ans[n, t]], t}], {t, -2, 1, 0.001}]], {n, 1, 4}] anim = Table[ 112 名前：Show[gr[1], gr[2], gr[3], gr[4], ViewPoint -> {20*Sin[s], 20*Cos[s], 0}], {s, Pi/24, 2*Pi, Pi/24}] もう良く分からない sol[y_] = Solve[y^4 + y == x^5 - x^3, x] ans[j_, t_] := x /. sol[t][[j]] Do[gr[n] = Graphics3D[Table[Point[{Re[ans[n, t]], Im[ans[n, t]], t}], {t, -2, 1, 0.001}]], {n, 1, 5}] anim = Table[Show[gr[1], gr[2], gr[3], gr[4], gr[5], ViewPoint -> {20*Sin[s], 20*Cos[s], 0}], {s, Pi/24, 2*Pi, Pi/24}] [] [ここ壊れてます] 113 名前：１３２人目の素数さん mailto:sage [2015/12/08(火) 21:05:53.29 ID:fafZe/EG.net] わけワカメw 114 名前：１３２人目の素数さん mailto:sage [2015/12/08(火) 21:20:58.41 ID:XscFjqkO.net] 連投ごめんなさい、 なんか面白くなってしまいました。 連投制限かMathematicaの記述がエラーに引っかかってしまい 他にも面白い曲線があるのですが投稿できずに残念です。 一番面白い？のはこの タイトル<虫>です。 sol[y_] = Solve[y^4 + y == x^4*y^2 - x^3, x] ans[j_, t_] := x /. sol[t][[j]] Do[gr[n] = Graphics3D[Table[Point[{Re[ans[n, t]], Im[ans[n, t]], t}], {t, -2, 1, 0.01}]], {n, 1, 4}] anim = Table[Show[gr[1], gr[2], gr[3], gr[4], ViewPoint -> {20*Sin[s], 20*Cos[s], 0}], {s, Pi/24, 2*Pi, Pi/24}] 115 名前：１３２人目の素数さん mailto:sage [2015/12/12(土) 20:04:05.97 ID:F+My5ir0.net] y^2 + y == x^3 - x^2 に、パラメタpをいれた 即席H-Function:x^3 - p x^2つくってみたがなw x /. Solve[y^2 + y == x^3 - p x^2, x]; Animate[Plot[Evaluate[Re[%]], {y, -1.5, 0.5}, PlotRange -> {-1, 1}, Filling -> {1 -> {{3}, LightPink}, 2 -> {Top, LightPink}}, Axes -> False], {p, 0.1, 0.9}, AnimationRunning -> False] 116 名前：１３２人目の素数さん mailto:sage [2015/12/12(土) 22:59:23.40 ID:mv17/rkp.net] 返信ありがとうございます。 早速ノートブックにコピペしていじってみましたが、 当方Mathematica４を使用しておりまして再現できませんでした。 sol[y_] = Solve[y^2 + y == x^3 - p*x^2, x] p = 1.5 ans[j_, t_] := x /. sol[t][[j]] Do[gr[n] = Graphics3D[Table[Point[{Re[ans[n, t]], Im[ans[n, t]], t}], {t, -2, 1, 0.01}]], {n, 1, 3}] anim = Table[Show[gr[1], gr[2], gr[3], ViewPoint -> {20*Sin[s], 20*Cos[s], 0}], {s, Pi/6, 2*Pi, Pi/6}] このようにpを変化させてできる曲線でイメージはよろしいのでしょうか？ さらに、H-Functionについて不勉強で理解できておりませんのでご教授頂ければ幸いです。 追伸 面白い曲線　タイトル <ケルベロス> sol[y_] = Solve[y^4 + y == x^4*y - x^3, x] ans[j_, t_] := x /. sol[t][[j]] Do[gr[n] = Graphics3D[Table[Point[{Re[ans[n, t]], Im[ans[n, t]], t}], {t, -2, 1, 0.01}]], {n, 1, 4}] anim = Table[Show[gr[1], gr[2], gr[3], gr[4], ViewPoint -> {20*Sin[s], 20*Cos[s], 0}], {s, Pi/24, 2*Pi, Pi/24}] 117 名前：１３２人目の素数さん mailto:sage [2015/12/13(日) 00:25:12.81 ID:axvNzcue.net] ケルベロス？　なんじゃこりゃ H(えっち)-Functionは、 "Mathematica版「おっぱい（曲面）方程式」" でググれば判る様な、無意味な遊びです。 118 名前：１３２人目の素数さん mailto:sage [2015/12/13(日) 02:25:11.92 ID:RhssocmV.net] 紳士の遊びですね？ｗ ケルベロスと名付けたのは、犬みたいな形で顔の部分が3本の曲線だからです。 なんか、>>1 のAAのキャラクターが出てくるアニメに出てきたような記憶が。 こんな曲線もあります <鼓> sol[y_] = Solve[y^4 + y == x^4 - x^2, x] ans[j_, t_] := x /. sol[t][[j]] Do[gr[n] = Graphics3D[Table[Point[{Re[ans[n, t]], Im[ans[n, t]], t}], {t, -2, 1, 0.001}]], {n, 1, 4}] anim = Table[Show[gr[1], gr[2], gr[3], gr[4], ViewPoint -> {20*Sin[s], 20*Cos[s], 0}], {s, Pi/24, 2*Pi, Pi/24}] 119 名前：１３２人目の素数さん [2015/12/19(土) 21:54:10.67 ID:mzCHiIu8.net] 人並みに扱われて勘違いしている田舎者チョンのニケア信教 ◆BkMobJaj1A こと山出しのニケア信教 ◆BkMobJaj1A。 その無知無教養さ加減はハゲの中でも下の下。 大阪大学や関西大学の学生からも馬鹿にされる為体。 St Marco Passion ◆BkMobJaj1A の臆病無比な卑劣さは、最底辺のハゲキモ蛆虫に劣る嫌らしさ！ こんな化け物老去勢オカマの存在を、我が国は赦しておいてはいけない。 早々にこの河馬ヅラをしたアナル変質者を朝鮮半島パクのもとへ強制送還させよう！！！ 120 名前：１３２人目の素数さん [2015/12/25(金) 14:34:39.76 ID:biauouUl.net] ParallelDoの中で添字付共有変数を使いたいんだけどどうすればよいですか？ 以下のように書いたらエラーが出てしまいました： Clear[a]; Do[a[i]=0, {i,1,10}]; (* 初期値 *) SetSharedVariable[a]; ParallelDo[a[i]=i,{i,1,10}] (* 並列計算 *) Do[Print[a[i]],{i,1,10}] 121 名前：１３２人目の素数さん mailto:sage [2015/12/28(月) 00:59:30.78 ID:NFbTjtr3.net] 添字付き変数渡しは、Headのせいか難しいっぽい。 答えにはなっていませんが、ListのPart扱いなら、 a = {}; Do[AppendTo[a, 0], {10}](*初期値*) SetSharedVariable[a]; ParallelDo[a[[i]] = i, {i, 1, 10}] (*並列計算*) Do[Print[a[[i]]], {i, 1, 10}] 122 名前：１３２人目の素数さん mailto:sage [2015/12/28(月) 11:25:53.74 ID:XH7cvZlD.net] >>118 どうもです 私も結局リファレンス通りListで渡しました ParallelDoでかなり高速化できました しかし並列カーネル4個追加で\73,440って高いっす 123 名前：１３２人目の素数さん [2016/02/03(水) 17:25:33.87 ID:REZsa8b9.net] 馬鹿な日本人ほど海外を崇めて、海外はすごいとかほめそやして 挙句の果てに日本人に留学を勧めたり海外の大学を勧めたりするが 実際は外人なんてクズばかりで勉強しないくせに選民意識だけは高いし 人種差別は平気でするしでどうしようもない 留学なんて奨励してる奴はキチガイですわ mathematicaの値段も頭おかしいのか？日本人が優秀だから日本でmathematicaを売ると 技術的にボロ負けしてしまうと危惧してるんだろうか？おそらくそうだろう。 アメリカがすごいのは大学や研究なんかじゃなくて金を荒稼ぎしてるからだからであって 日本が真似すると大やけどどころか沈没だよ。馬鹿な日本人ほどアメリカの真似したがるけど。 日本がすごいのは日本人が頑張ってきたからだよ。 124 名前：１３２人目の素数さん mailto:sage [2016/02/03(水) 20:57:46.08 ID:VcbUmeJE.net] Mathematica は所属してる大学とか組織でライセンス取るものだろ。 125 名前：１３２人目の素数さん mailto:sage [2016/02/04(木) 09:53:50.95 ID:2sg1Gtew.net] Mathematicaで理工学部程度の数学が学べるおすすめの書籍を教えてください。 126 名前：１３２人目の素数さん mailto:sage [2016/02/04(木) 09:56:19.03 ID:2sg1Gtew.net] または理工学部程度の数学でMathematicaを学べる書籍を教えてください。 127 名前：１３２人目の素数さん mailto:sage [2016/02/04(木) 11:25:15.90 ID:A+CLB+hi.net] 対して数出ていないし古いのしかないから 図書館検索で全部借りればいい 128 名前：１３２人目の素数さん mailto:sage [2016/02/04(木) 11:27:18.89 ID:A+CLB+hi.net] 俺のお気に入りはこれ https://www.kinokuniya.co.jp/f/dsg-01-9784000054478 129 名前：１３２人目の素数さん mailto:sage [2016/02/04(木) 11:46:32.82 ID:2sg1Gtew.net] >>125 おもしろそうです 130 名前：１３２人目の素数さん [2016/02/04(木) 18:55:14.58 ID:h6GhpBCQ.net] ませまちか　は　○○○で　○○す物 131 名前：１３２人目の素数さん mailto:sage [2016/02/05(金) 19:02:15.13 ID:VmFpJCfn.net] センター試験や大学入試問題をコマンド一発で解けるのは感動 証明と図形以外の問題で解けないのはないだろって 逆に言えばその二つが出来れば完璧なのに 132 名前：１３２人目の素数さん [2016/02/05(金) 20:55:05.68 ID:4S9VmYC2.net] そんなもの解いてもしょうがない 133 名前：１３２人目の素数さん mailto:sage [2016/02/06(土) 04:07:51.62 ID:CZQRzoj6.net] >>122 おすすめ 『レクチャーズ オン Mathematica』 もくじなどのpdfがある レクチャーズ オン Mathematica・サポートページ www.math.titech.ac.jp/~kawahira/courses/mathematica.html 134 名前：１３２人目の素数さん mailto:sage [2016/02/10(水) 17:25:12.67 ID:xg3fAtil.net] {i,j,k} を要素に持つリストをつくりたいんですがうまい記述の仕方あるでしょうか？ i=0〜Nx j=0〜Ny k=0〜Nz の範囲の整数として { {0,0,0}, {1,0,0}, {2,0,0}, …　…, {3,2,5}, …, {Nx,Ny,Nz} } のような感じです． 135 名前：１３２人目の素数さん mailto:sage [2016/02/10(水) 18:49:52.80 ID:WSI9K ] [ここ壊れてます] 136 名前：NQ8.net mailto: >>131 Flatten[Table[{i, j, k}, {i, 0, 4}, {j, 0, 5}, {k, 0, 6}], 2] とか？ [] [ここ壊れてます] 137 名前：１３２人目の素数さん mailto:sage [2016/02/10(水) 19:08:59.13 ID:xg3fAtil.net] >>132 ありがとうございましゅ！ 138 名前：１３２人目の素数さん mailto:sage [2016/02/10(水) 19:09:16.46 ID:xg3fAtil.net] ×ありがとうございましゅ！ ○ありがとうございます！ 139 名前：１３２人目の素数さん mailto:sage [2016/02/18(木) 23:46:48.85 ID:hcnN01i/.net] 1/0となりますっていうエラーの出る分母の値（0と認識される値?、結果が無限大と認識される値?）ってどれくらいでしょうか。 電子の質量入れて計算すると飛びます。 分母のオーダーが10^(-31-α)とかでもだめですか？ 140 名前：１３２人目の素数さん mailto:sage [2016/02/18(木) 23:47:54.83 ID:cFQO+ZMK.net] >>135 具体例を 141 名前：１３２人目の素数さん mailto:sage [2016/02/19(金) 00:36:33.18 ID:BVdQKtrr.net] ええとすみません。たとえば、 1/Exp[-x] とするとxがどの辺のときに計算が放置されるのでしょうか？ 142 名前：１３２人目の素数さん mailto:sage [2016/02/19(金) 09:45:48.44 ID:sj9/NYrE.net] x=10.^15 は大丈夫 x=10.^15*4 でアンダーフロー 143 名前：１３２人目の素数さん mailto:sage [2016/02/19(金) 19:27:51.70 ID:BVdQKtrr.net] とｎ 電子の質量が入ってくるときとかどうするんだろうね… 144 名前：１３２人目の素数さん mailto:sage [2016/02/19(金) 19:55:53.52 ID:9jEB8QMH.net] >>139 その、電子の質量を入れて計算させると1/0エラーが出るっていう式をそのまま書いてみて 145 名前：１３２人目の素数さん mailto:sage [2016/02/28(日) 18:27:57.65 ID:5GJ+OBvk.net] u'=-uなる微分方程式を初期値hを0から1まで変えながら図にするのに 下のようにしたのですが 上手く行きません。何が悪いのでしょうか Do[{ sol1 = NDSolve[{ u'[t] == -2 u[t], u[0] == h }, u[t], {t, 0, 3}] Plot[Evaluate[u[t] /. %], {t, 0, 3}] }, {h, 0, 1}] 146 名前：１３２人目の素数さん mailto:sage [2016/02/28(日) 19:05:21.41 ID:H7dcneBr.net] Show[Table[ sol1 = NDSolve[{u'[t] == -2 u[t], u[0] == h}, u[t], {t, 0, 3}] ; Plot[Evaluate[u[t] /. sol1], {t, 0, 3}], {h, 0, 1, 0.1}]] or list = Evaluate[ Table[sol1 = NDSolve[{u'[t] == -2 u[t], u[0] == h}, u[t], {t, 0, 3}]; u[t] /. sol1, {h, 0, 1, 0.1}]]; Plot[list, {t, 0, 3}] こんなんみたいのかな 147 名前：１３２人目の素数さん mailto:sage [2016/02/28(日) 19:35:00.37 ID:5GJ+OBvk.net] すいません %→sol1とコンマ忘れでいけましたね 大変ありがとうございました 148 名前：１３２人目の素数さん mailto:sage [2016/03/11(金) 21:14:10.00 ID:LiuUp7W3.net] 質問です 確率変数x,yとzがあり、それぞれ平均値mx,my,mzで標準偏差sx,sy,sz の正規分布をもって変動するものとします。 このとき例えばz x/(x+y)の値がどのように分布するかを簡単にグラフ化 する手法と、その標準偏差を簡単に求める方法はあるでしょうか。 分布図はモンテカルロ法、統計値は一次近似の数式計算とかででしたが、 忘れた頃にやるので、火急の時には厳しい･･･ Ver10以下の範囲で、何かヒントがあればお願いします。 149 名前：１３２人目の素数さん mailto:sage [2016/03/12(土) 21:40:18.28 ID:r9WZX4oj.net] d = TransformedDistribution[ z x/(x + y), {x \[Distributed] NormalDistribution[mx, sx], y \[Distributed] NormalDistribution[my, sy], z \[Distributed] NormalDistribution[mz, sz]}]; example = {mx -> 1, sx -> .2, my -> 2, sy -> .7, mz -> 2.9, sz -> 1.1}; PDF[d /. example, x] d /. example // StandardDeviation (\[Distributed]はESCdistESCと打っても良い。） で確率密度関数と標準偏差が出るはずだけど、式が複雑すぎて 計算しきれていないみたい。まあゼロ割になりかねない式だから、 標準偏差があるかどうか微妙だが。 なのでモンカルを使うと、 data = RandomVariate[d /. example, 100000]; Histogram[ 150 名前：data] data // StandardDeviation という感じ。でもSDは安定しないな。やはり式の形が… スムーズなグラフを描きたいのなら、 Plot[Evaluate@PDF[data // SmoothKernelDistribution, x], {x, -1, 5}] とかやれば良い。 [] [ここ壊れてます] 151 名前：１３２人目の素数さん mailto:sage [2016/03/12(土) 23:42:05.46 ID:o9XKea1M.net] ありがとうございます。 式自体は単純に、基準電圧zを抵抗y,xで分圧した電圧で、 その分布を求める場合ですが、モンカルだよりですかね。 変数毎に乱数発生をさせて、式に代入させていましたが、 いましたが、TransformedDistribution->RandomVariate の流れは参考になりました。 あとSmoothKernelDistributionも(^^)/ 152 名前：１３２人目の素数さん mailto:sage [2016/03/28(月) 20:52:34.19 ID:Y/nPZlRi.net] Maple使ってる人いる？ eliminateって機能があって、 MarhematicaのEliminate関数と同じ、連立方程式からの変数消去なんだけど、 あまりにアホすぎて使い物にならんねコレ 153 名前：１３２人目の素数さん mailto:sage [2016/03/28(月) 23:40:49.43 ID:gs+YwEEW.net] 使い物になるという話かと思ったら、使い物にならないという話かよ。 154 名前：１３２人目の素数さん mailto:sage [2016/03/29(火) 09:21:30.33 ID:MS5DySjF.net] MCのことモンカルって言う人はじめてみたわ 155 名前：１３２人目の素数さん mailto:sage [2016/03/31(木) 12:58:22.72 ID:cpRtmv87.net] Mapleのeliminateってまだ邦訳されてないから、最近追加された機能かな？ 個人的にMathematicaで一番重宝してる変数の整理、消去がMapleだと今市だな。 156 名前：１３２人目の素数さん mailto:sage [2016/03/31(木) 16:47:29.31 ID:lQq9ZfEU.net] ell = 8; m = 0; Integrate[ SphericalHarmonicY[ell, m, \[Theta], 0]* Sin[\[Theta]]^2 (1 + q Sin[\[Theta]]^2)*Sin[\[Theta]], {\[Theta], 0, \[Pi]}] 俺の環境 (ver10.4, Windows7_64bit) でこの計算やらせると落ちてしまう． ell=6だと回る． なんで〜〜（泣） 157 名前：149 mailto:sage [2016/03/31(木) 17:12:04.92 ID:lQq9ZfEU.net] ver10.0だと回る ぐぬぬ… 158 名前：１３２人目の素数さん mailto:sage [2016/04/03(日) 10:21:42.74 ID:+ao9jRpb.net] >>151 割っただろ 159 名前：１３２人目の素数さん mailto:sage [2016/04/03(日) 10:44:28.95 ID:QaspXQZ2.net] やばそうな割り算無さそうなんだけどなぁ． ver 10.3 (Windows7_64bit) でも落ちる． ver 8 (student ed., Windows XP pro) ではちゃんと計算される． 10.1と10.2はアップデートしてしまったから不明． ver 10.3しか確かめてないけどMacProでも落ちるからハードが原因じゃない気がする． はぁ〜とりあえずver 10.0入ってる少し古い計算機使います． 160 名前：１３２人目の素数さん [2016/04/04(月) 12:24:52.18 ID:8UdNWihF.net] 10.4だけど確かに落ちる。展開するとどんな形になるのかと思って、FunctionExpandを適用してみた。 ell = 8; m = 0; FullSimplify[FunctionExpand[SphericalHarmonicY[ell, m, theta, 0] Sin[theta]^2 (1 + q Sin[theta]^2) Sin[theta]]] 三角関数の組み合わせ。qに具体的な値を入れると積分できるか？ Table[Integrate[%, {theta, 0, Pi}], {q, -1, 1, 1/5}] ->{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} 答えは0っぽい。で、この後気まぐれに元の式をもう一度計算してみた。 Integrate[SphericalHarmonicY[ell, m, theta,0] Sin[theta]^2 (1 + q Sin[theta]^2) Sin[theta], {theta, 0, Pi}] ->0 今度は計算できる！不思議！ 161 名前：１３２人目の素数さん mailto:sage [2016/04/06(水) 13:03:16.89 ID:gULLYejW.net] なんだそりゃw 外部からアクセスできないキャッシュの内容によって挙動が変わるのか。 報告しとくかな。 162 名前：１３２人目の素数さん [2016/04/19(火) 20:07:51.97 ID:YSV4FQMK.net] imgur.com/LKDrVCj.jpg ↑Histogramでヒストグラムを描きました。 階級値の数字が300, 400, 500, 600と100間隔で表示されていますが、 300, 350, 400, 450, 500, 550, 600と表示させたいです。 どうすればいいのでしょうか？お願いします。 163 名前：１３２人目の素数さん mailto:sage [2016/04/20(水) 19:34:29.49 ID:EJjD1eVy.net] Ticksで座標軸の目盛りを明示的に指定すれば？ ex.横軸の目盛り範囲[100,800]隔50を指定する sample = RandomVariate[NormalDistribution[500, 100], 1000]; Histogram[sample, {100, 800, 50}, Ticks -> {Range[100, 800, 50], Automatic}] 164 名前：１３２人目の素数さん [2016/04/20(水) 20:43:33.60 ID:V7TU7EwI.net] >>158 imgur.com/XMJ4pgw.jpg ありがとうございました。 でも↑の画像のようになってしまいました。 数字は階級の両端の数字ではなく階級値を表示させたいんです。 165 名前：１３２人目の素数さん [2016/04/20(水) 20:45:48.41 ID:V7TU7EwI.net] >>158 ありがとうございました。 Ticks -> {Range[125, 775, 50], Automatic}] みたいにやればできるんですね。 ありがとうございました。 166 名前：１３２人目の素数さん mailto:sage [2016/04/25(月) 22:30:29.53 ID:nOxB4NtD.net] >>125 その本良かったんですか。買ったけど読んでない。いまではmathematicaは やめてmapleを使うことにした。 167 名前：１３２人目の素数さん [2016/04/26(火) 07:33:53.58 ID:Bqjxy5C+.net] Mapleのほうが使いやすいよね。 168 名前：１３２人目の素数さん mailto:sage [2016/05/10(火) 17:58:44.15 ID:3vtuhWyp.net] A = { {1,2}, {3,4} } // MatrixForm A.A (*これだと行列計算はしてくれない*) これ一度 MarixFormでラッピングしてしまった後から、コード上でアンラップする方法って無いんですかね？ 実用性云々(ラッピングしなきゃいいじゃん) ではなく単純に疑問です。 A = ToExpression @ StringSplit[#," "]& @ StringSplit[#,"\n\n"]& @ ToString[A] 文字列にして読み直ししたら一応いけたんですが、改行は2回スペースも2回続きそんなのをあてにしていいものやら... もっと単純に一発でできてもよさそう。 169 名前：１３２人目の素数さん mailto:sage [2016/05/10(火) 21:58:53.04 ID:N5W3xcsN.net] >>163 A[[1]] 170 名前：１３２人目の素数さん mailto:sage [2016/05/10(火) 23:35:12.79 ID:ISuglzw8.net] 行列を代入したい、でもMatrixFormで確認もしたい、という時は、 (A={{1,2},{3,4}})//MatrixForm みたいにやるな。 171 名前：１３２人目の素数さん mailto:sage [2016/05/11(水) 01:10:37.40 ID:xWSeBPqz.net] >>164, >>165 ありがとうございます。 >>164 「〜Form の類いは、関数というよりはデータ構造構築子で、元データを1個だけ格納してる」 って理解で合ってるでしょうか この辺りの詳細って WolframのWEBドキュメントに載ってます？ 172 名前：１３２人目の素数さん mailto:sage [2016/05/11(水) 16:10:36.87 ID:fSvrr8qW.net] >>166 Part 関数を参照 173 名前：１３２人目の素数さん mailto:sage [2016/05/11(水) 19:31:46.17 ID:xWSeBPqz.net] >>167 ありがとうございます。 174 名前：１３２人目の素数さん mailto:sage [2016/05/11(水) 23:05:47.89 ID:QU9FYAS5.net] 強引だけどA = { {1,2}, {3,4} } // MatrixFormで A.A を計算するだけならMatrixFormでThreadもあり A.A // Thread[#, MatrixForm] & 175 名前：１３２人目の素数さん mailto:sage [2016/05/12(木) 10:39:58.32 ID:QuIsVXVk.net] プログラミングで普通に思い浮かべるThreadとは関係無くて、 そのまんま「糸」でデータを縫い合わせる関数なのね。 基本的なのにまだまだ知らない関数があるなあ。 176 名前：１３２人目の素数さん [2016/06/15(水) 11:57:22.83 ID:uO5svwl0.net] 梅津千春＝立岡咲樹　かかってこいヨーイチ 177 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:37:21.73 ID:KWnWf/3r.net] ￥ 178 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:37:47.07 ID:KWnWf/3r.net] ￥ 179 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:38:14.77 ID:KWnWf/3r.net] ￥ 180 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:38:44.15 ID:KWnWf/3r.net] ￥ 181 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:39:12.19 ID:KWnWf/3r.net] ￥ 182 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:39:33.61 ID:KWnWf/3r.net] ￥ 183 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:40:03.74 ID:KWnWf/3r.net] ￥ 184 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:40:56.9 ] [ここ壊れてます] 185 名前：9 ID:KWnWf/3r.net mailto: ￥ [] [ここ壊れてます] 186 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:41:17.39 ID:KWnWf/3r.net] ￥ 187 名前：￥ ◆2VB8wsVUoo mailto:sage [2016/06/15(水) 18:41:37.11 ID:KWnWf/3r.net] ￥ 188 名前：￥ ◆2VB8wsVUoo mailto:age [2016/06/16(木) 18:30:38.78 ID:PgonDtoO.net] ￥ 189 名前：￥ ◆2VB8wsVUoo mailto:age [2016/06/18(土) 02:40:26.20 ID:3GXqstO/.net] ￥ 190 名前：１３２人目の素数さん [2016/06/18(土) 21:01:21.49 ID:bG4TQCiY.net] 大類昌俊 関西すうがく徒のつどい 関東すうがく徒のつどい https://s32.postimg.org/r0ygqs3dx/image.jpg https://s32.postimg.org/ar8anvspx/image.jpg https://s32.postimg.org/ql33m4579/cropped_6acae72b0e380001645f78ed66f39459.jpg https://s31.postimg.org/f0wok1s6j/Sommer201942_20_Brutaler20_Partisanenkampf2.jpg https://s31.postimg.org/6itpm717f/attach24_787637.jpg https://s31.postimg.org/mc9p02tmj/Sommer201942_20_Brutaler20_Partisanenkampf.jpg 191 名前：１３２人目の素数さん mailto:sage [2016/07/08(金) 17:35:43.32 ID:cG8o2Rlw.net] a = 1.234567890; N[a, 10] の出力が 1.23457 になってしまうのだけれど， N関数って機械精度の数を渡されると常に6桁表示になるの？ 192 名前：１３２人目の素数さん mailto:sage [2016/07/08(金) 17:56:49.57 ID:cG8o2Rlw.net] Nのせいじゃなくて標準出力が機械精度に対して6桁表示なのね… InputForm[a]もしくはNumberForm[a,10]で表示された 193 名前：１３２人目の素数さん mailto:sage [2016/07/09(土) 07:07:00.48 ID:fDkSF8u6.net] 10.4.0 for Mac OS X x86 (64-bit) (March 1, 2016) なんだけど、 フォントのサイズ変更したいから www.math.titech.ac.jp/~kawahira/courses/mathematica/01.pdf にある > メニューから「編集」→「環境設定...」* → 「詳細」→ 「オプションインスペクタを開く」→ 「書式設定」→ 「フォント設定」→ 「FontSize(サイズ)」とすすみ,これをデフォルトの 12 から適宜増や す.「適用」ボタンを押すとフォントサイズが大きくなる. を行おうとしたが、「適用」ボタンがグレーのままで押しても適用出来ない。 どうしたら良いんだコレ？ i.imgur.com/BEHcB9D.png 194 名前：１３２人目の素数さん mailto:sage [2016/07/15(金) 12:32:01.25 ID:C6cVoJv7.net] list1 = { {x1, y1, z1}, {x2, y2, z2} }; list2 = { {v1, v2}, {w1, w2} }; こういうタイプの2つのリストの各要素を結合して， { {x1, y1, z1, {v1, v2} }, {x2, y2, z2, {w1, w2} } } というリストを作るために FlattenAt[#, 1] & /@ Partition[Riffle[list1, list2], 2] と記述してるんですが，より高速な方法があれば教えてください 195 名前：186 mailto:sage [2016/07/16(土) 00:47:18.27 ID:xDKztG0Z.net] Clear[result]; result = Transpose@list1; result = AppendTo[result, list2]; result = Transpose@result; で速度が2倍ぐらいになりました 196 名前：１３２人目の素数さん mailto:sage [2016/07/16(土) 06:13:38.45 ID:GrLp8R2o.net] 速度はわからんけど、 MapThread[Append[#1, #2] &, {list1, list2}] かなあ。 197 名前：186 mailto:sage [2016/07/16(土) 18:15:35.46 ID:xDKztG0Z.net] >>190 更に速くなりました 記述もスマートですしこれを使いたいと思います ありがとうございました 198 名前：１３２人目の素数さん mailto:sage [2016/07/20(水) 21:06:00.17 ID:K8aqavNL.net] ver10.4 で並列カーネルの状態ダイアログから並列カーネルを全部閉じると ↓みたいにエラーが発生するんだけど同じひといる？ 1つずつ閉じても最後の1個を閉じるところで同じエラーがでる i.imgur.com/kCogne9.png 199 名前：１３２人目の素数さん mailto:sage [2016/08/22(月) 22:06:54.62 ID:+pVfiDmq.net] TEST 200 名前：１３２人目の素数さん [2016/08/23(火) 06:38:21.34 ID:2MEC67az.net] 11のアイコンしょぼすぎ

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef