〓 Mathematica 七 〓
, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,）)
! | ! .,'-.{ ! !|; |`､.}ﾞ!.! |. ! ヽ.
', ', |Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l ! | ﾞ',
ヽ､', l:!Kノ}. f:_.)iﾞi: ﾘ ! l ル
| l!iヽﾞ- ' , .!__:ﾉ ﾞ ,ﾘ l ﾘ'´
. ',|!!､ r‐┐ ｀ ノ'. /,ｲ ／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣
'i!ﾞ､ヽ､ ﾞｰ' _, ィ,:',:''´ ＜ Mathematica に関する話題はここに書いてね！
ﾞ:､ｨ､jヾｰ::: 'iﾍ .ノ',ﾘ. ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
,､- '´ ヽ､ﾞ､ { ｀＞"､
/＼＼ ', } //｀ヽ
過去スレ
〓Mathematica〓
http://science2.2ch.sc/test/read.cgi/math/985023298/
〓 Mathematica 2 〓
http://science3.2ch.sc/test/read.cgi/math/1078534285/
〓 Mathematica 3 〓
http://science6.2ch.sc/test/read.cgi/math/1121413040/
〓 Mathematica 四 〓
http://kamome.2ch.sc/test/read.cgi/math/1197270001/ (鯖飛びでログ消滅)
http://mimizun.com/log/2ch/math/1197270001/
〓 Mathematica 5 〓
http://kamome.2ch.sc/test/read.cgi/math/1285859504/
〓 Mathematica 伍 〓 (実質6)
http://wc2014.2ch.sc/test/read.cgi/math/1320969748/
>どんなプロセスを経てるのかはブラックボックスだが。
Traceを使えば途中のプロセスを確認できます(見方はヘルプ参照)。
ex.
Simplify[Sum[
n Binomial[m, n] FactorialPower[m, n] FactorialPower[l - m,
m - n], {n, 1, m}]/FactorialPower[l, m] // FunctionExpand,
Assumptions -> Element[m, Integers]] // Trace
両方共Linuxだしいけそうなんだけど
が10000を超えたみたいだね
このスレに投稿したことある人っているの？
wolframaplha（Web版）で、
１．Mathematicaのコマンドを実行するにはどうしたらいいのでしょうか？
やろうとすると、認識できませんでしたとか何とか、拒否されます。
２．変数名に、たとえばi0など使うと、ベッセル級数であると勝手に解釈されてしまいます。
この解釈をやめさせる（エスケープする）には、どうしたらよいでしょうか？
変数名のi0が、ベッセル級数？<なら変名数変えろよ。
おら心折れて、やる気△1000％なんで↓↓↓対応よろ。
図解が便利だって言ってたけど結局研究進められず
田舎に帰ってったわ 無駄な買い物だったね
これってなんで出るんだ？
そうモアレ
これって回避する方法ってなんかある？
なんかうまくいきません
ノートブックといわれるファイルに書くのでしょうか？
それともテキストエディタを使ってコードを書いて、ノートブックから読み込むのでしょうか？
1+1
Enter
とか入力したらどうなるん？
つーかそもそもどういう環境？
それは門外不出
遅すぎますか？
10.0.0.2ってMACしかなくね？
アメリカに遅れること約1ヶ月か。
plot[....]
でぐらふがかけたのに（いまでも）
いまでは
ss := plot[....];
ではぐらふがかけないね
ss
出かける。
いつからこうなったんだろう
コマンドのマセマティカは動くのですが、WINDOW（X)画面はINITIALZEで灯ったままです。
そのまえは動いていたのですが、同じような経験のかたはいませんか？
やすいから まず CARD（１５００円）をかって再インストールするんだな
あとは誰かよく知っている奴に教えても懶惰な
あ、訂正します：
Raspberry Pi2が届きました。
RaspbianをインストールしMathematicaを使ってみましたが問題がありました。
最初は起動できるのですが、一度終了して、また起動しようとすると、
「Initializing Kernels ...」というメッセージのところでフリーズ
してしまうんですよね。
調べたところまだ修正されていないバグだそうです。
対策は、WolframというGUIではないMathematicaを起動して、
PacletUpdate["CloudObject"]
を評価すると以後、問題なく使えるようになるようです。
ちなみに、Raspberry Pi2の性能でもほとんど速度的に問題なく
Mathematicaを使えるようです。
SDカードにも注意が必要です。
Amazon.co.jpで推奨のTranscendの32GBのmicroSDカードを最初に買った
のですが、エラーが出て起動できなくなりました。調べたら他にも同じ
人がいて、どうも相性が良くないようです。相性というか、Raspberry Pi2
かmicroSDカードのどちらかがmicroSDカードの仕様を満たしていないか、
microSDカードの仕様自体に問題があるかですよね。ひどい話です。
SAMSUNGの32GBのmicroSDカードを新たに買いましたが、全く問題なく
使用できています。
家でもMathematicaが安く使えるっていいですよね。
下記の計算で、z1=z2となりますが、前に係数epが付いただけで、
y1=y2、y1=y3となりません。
なぜでしょうか？
ep =.;
a =.;
z1 = Exp[-Log[a]]
z2 = Exp[Log[1/a]]
z1 - z2
y1 = Exp[ep/2*(Log[1/a])]
y2 = Exp[ep/2*(-Log[a])]
y3 = Exp[-ep/2*(Log[a])]
FullSimplify[y1 - y2]
FullSimplify[y1 - y3]
すいませんが、以下を教えて下さい。
質問１、
a^(k)=Exp[k*Log[a]]
ですが、（下記の計算参照）
y1 = (1/a)^(k/2)を、変換して
y2 = Exp[k/2*(Log[1/a])]
か、もしくは
y3 = Exp[k/2*(-Log[a])]
に書き換えることは可能でしょうか？
可能でしたら、
(1/a)^(k/2)→ Exp[k/2*(Log[1/a])]
(1/a)^(k/2)→ Exp[k/2*(-Log[a])]
に書き換える方法を知りたいです。
x1 = a^(k)
x2 = Exp[k*Log[a]]
x1 - x2
y1 = (1/a)^(k/2)
y2 = Exp[k/2*(Log[1/a])]
y3 = Exp[k/2*(-Log[a])]
下記の計算で、z1=z2となりますが、前に係数kが付いただけで、
y1=y2、y1=y3となりません。
なぜでしょうか？
k =.;
a =.;
z1 = Exp[-Log[a]]
z2 = Exp[Log[1/a]]
z1 - z2
y1 = Exp[k/2*(Log[1/a])]
y2 = Exp[k/2*(-Log[a])]
y3 = Exp[-k/2*(Log[a])]
FullSimplify[y1 - y2]
FullSimplify[y1 - y3]
多価性
多価関数の分岐問題を無視するなら、FullSimplifyでなく
PowerExpandを強行作用させる。
>>43
組込関数PowerのProtect属性を一度外してオーバーロード
(再定義)する。
その際には、Exp[k Log[a]]がa^kに自動的に簡約されない
ように、HoldFormなどのラッパーをかませる。
またaがEだと無限地獄に陥るのでUnsameQなどの判定も必要。
ご回答有難う御座います。
＞PowerExpandを強行作用させる。
解りました。
>>43
やってみます。
言ってみたものの結構むずい、再帰終端があまいですが単純なサンプル例。
Unprotect[Power, Log];
Power[a_/;a=!=E, b_ /;b=!=-1]:=(b Log[a])//HoldForm[Power[E, #]] &
Log[Power[a_,-1]]:=-Log[a]
Protect[Power, Log];
組込関数をいじっているので、副作用に注意が必要です。
(*Test*)
(1/a)^(k/2) (* -> Exp[-1/2 k Log[a] *)
% - Exp[k/2 Log[1/a]](* -> 0 *)
Exp[-Log[a]] - Exp[Log[1/a]](* -> 0 *)
ご回答有難う御座います。
自分でやって、直ぐに諦めました。
ご回答を見て、「やっぱ難しい。」ことが確認出来ました。
副作用に注意の上、応用して、使用させて頂きます。
助かりました。
コンピューターとか数学の真の姿を見た気分だったw。当時は。
早速コピーして持って帰ったｗ
商用だなんてえれー自信過剰なやつだと思ったもんだ。
機能にもビビったが何よりも驚いたのは値段だな。adobe製品より高いじゃんって
スノレパ切り捨てかよ
６０００円もあればラズパイで使えるのに
NEXTみたいだｗ
高校生が宿題やるくらいだったら余裕って程度のマシンスペックだし個人用ならアリかもな
2Bでどうにかまともに操作できるスペック
BまではWolfram(CUI)でないと苦痛なレベル
Import[
"IMGA.avi", {{
"BitDepth",
"ColorSpace",
"Duration",
"FrameCount",
"FrameRate",
"ImageSize",
"VideoEncoding"
}}]
Import["IMGA.avi", {"Frames", 1}] <-この結果が欲しい
(* out *)
{8, RGBColor, 7.173776069290165`, 172, 23.976215362548828`, {1280,
720}, "MJPG"}
Import::fmterr: データをvideo形式でインポートすることができません． >>
こんな感じで読み込めないんだよ
これは文法的問題なのか動画のエンコードの問題なのか
一応Mathematica9で読み込めるMotionJPEGでエンコードしたんだが
誰か解決出来る人いない？
ちなみに，H261,H263などのコーデックは試した.
10.2で手持ちのaviファイルで実行してみたけど、2つとも正常に動作した。
なので、エンコードの問題だね。
下記HPのP１０ 上から６行目計算をしたいのですが、入射・散乱光子の偏極の部分（下記コード参照）を、このように計算したらいいのか？解りません。
教えて下さい。
http://members3.jcom.home.ne.jp/nososnd/qed/comp.pdf
(*コンプトン散乱～クライン・仁科の公式 P１０ 上から６行目計算*)
s =.;
u =.;
(*ガンマ行列*)
g[0] = {{1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, -1, 0}, {0, 0, 0, -1}};
g[1] = {{0, 0, 0, 1}, {0, 0, 1, 0}, {0, -1, 0, 0}, {-1, 0, 0, 0}};
g[2] = {{0, 0, 0, -I}, {0, 0, I, 0}, {0, I, 0, 0}, {-I, 0, 0, 0}};
g[3] = {{0, 0, 1, 0}, {0, 0, 0, -1}, {-1, 0, 0, 0}, {0, 1, 0, 0}};
(*入射・散乱電子・光子*)
sl[k] = g[0]*k0 + g[1]*(-k1) + g[2]*(-k2) + g[3]*(-k3);
sl[k'] = g[0]*k0' + g[1]*(-k1') + g[2]*(-k2') + g[3]*(-k3');
sl[pi] = g[0]*pi0 + g[1]*(-pi1) + g[2]*(-pi2) + g[3]*(-pi3);
sl[pf] = g[0]*pf0 + g[1]*(-pf1) + g[2]*(-pf2) + g[3]*(-pf3);
e4 = IdentityMatrix[4];
ms = m*e4;
(*入射・散乱光子の偏極 ここの部分が解りません*)
(*sl[epsilon]=g[0]*epsilon0+g[1]*(-epsilon1)+g[2]*(-epsilon2)+g[3]*(-epsilon3);
sl[epsilon']=g[0]*epsilon0'+g[1]*(-epsilon1')+g[2]*(-epsilon2')+g[3]*(-epsilon3');*)
sl[epsilon] = e4;
sl[epsilon'] = e4;
(*トレース部分の計算*)
y1 = 0; y2 = 0;
For[x = 0, x <= 3, x++,
For[y = 0, y <= 3, y++,
s1 = Tr[sl[epsilon'].sl[k].sl[epsilon].(sl[pi] + ms).sl[epsilon].sl[k].sl[epsilon'].(sl[pf] + ms)];
y1 = y1 + s1;
]];
(*計算結果の整理 入射・散乱光子の偏極でのこの部分も解りません*)
y1 = y1 //. {pi1 -> 0, pi2 -> 0, k0 -> pi3, k1 -> 0, k2 -> 0, k3 -> -pi3,
pf0 -> pi0, pf1 -> pi3*Sqrt[1 - z^2], pf2 -> 0, pf3 -> pi3*z, k' 0 -> pi3,
k' 1 -> -pi3*Sqrt[1 - z^2], k' 2 -> 0, k' 3 -> -pi3*z, pi0 -> (s + m^2)/(2 Sqrt[s]),
pi3 -> (s - m^2)/(2 Sqrt[s]), z -> 1 + t/(2 pi3^2), t -> 2 m^2 - s - u, epsilon0 -> 1,
epsilon1 -> 1, epsilon2 -> 1, epsilon3 -> 1, epsilon0' -> 1, epsilon1' -> 1, epsilon2' -> 1,
epsilon3' -> 1};
s = 2*k*pi + m^2;
u = -2*k'*pi + m^2;
Print["(*計算結果*)"];
Print[Simplify[y1]];
Print["(*正しい計算結果*)"];
8*(pi*k)*(pi*k')
すいません。最後の正しい計算結果は、以下となるはずです。
Print["(*正しい計算結果*)"];
8*(pi*k)*((pi*k') + 2*(epsilon'*k)^2)
＞ベクトルとスカラーの区別もつかんのか
この部分のことでしょうか？
(*入射・散乱光子の偏極 ここの部分が解りません*)
(*sl[epsilon]=g[0]*epsilon0+g[1]*(-epsilon1)+g[2]*(-epsilon2)+g[3]*(-epsilon3);
sl[epsilon']=g[0]*epsilon0'+g[1]*(-epsilon1')+g[2]*(-epsilon2')+g[3]*(-epsilon3');*)
sl[epsilon] = e4;
sl[epsilon'] = e4;
ベクトルにしました。
http://www43.tok2.com/home/iq188/
epsilon0 -> 1,
epsilon1 -> 1, epsilon2 -> 1, epsilon3 -> 1, epsilon0' -> 1, epsilon1' -> 1, epsilon2' -> 1,
epsilon3' -> 1
この部分を、直せば、計算できるのでしょうか？
入射光子は
k0 -> pi3,
k1 -> 0,
k2 -> 0,
k3 -> -pi3,
で、これに直交するので、たぶん
epsilon0 -> 0,
epsilon1 -> 1,
epsilon2 -> 1,
epsilon3 -> 0,
となると思います。
では、散乱光子に
k' 0 -> pi3,
k' 1 -> -pi3*Sqrt[1 - z^2],
k' 2 -> 0,
k' 3 -> -pi3*z,
に直交する
epsilon0' -> ?,
epsilon1' -> ?,
epsilon2' -> ?,
epsilon3' -> ?;
値は、あるのでしょうか？考え方が間違っているでしょうか？
>>69のpdfを最後まで読めば
Mathematica9でAVIのデータって対応コーデックでエンコードしてるはずなのに
読み込めない事あるよね
PDFのP13とP14からε(1)、ε‘(1)（=εd[1]）, ε(2)、ε‘(2)は以下のように思います。
P９のε・ε＝－１になるような“ε”が、具体的にどのような行列になるのか？解りません。
ε[1] = {0, 1, 0, 0};
ε[2] = {0, 0, 1, 0};
εd[1] = {0, Cos[theta], 0, Sin[theta]};
εd[2] = {0, 0, 1, 0};
>ε[1] = {0, 1, 0, 0};
>ε[2] = {0, 0, 1, 0};
>εd[1] = {0, Cos[theta], 0, Sin[theta]};
>εd[2] = {0, 0, 1, 0};
各々について
>P９のε・ε＝－１
の意味での “ε[1]・ε[1]” 等を計算してみ
お返事有難う御座います。
各々について
>P９のε・ε＝－１
の意味での “ε[1]・ε[1]” 等を計算してみ
考えたのですが、ε=(0,ε)の定義というか意味が解らないです。
pi=(m,0)なら、電子の静止質量ｍとし、pi=(m,0,0,0,)とわかるのですが、
ε=(0,ε)なら、ε= {0, 1, 0, 0}を素直に入れると、ε= {0,0,1, 0, 0}になってしまいます。
またεは、スカラーや行列ではなく、ベクトルなら、ε・ε＝－１になるのは,、下記のように虚数しかないと思います。
しかし、そうするとε[2]・ εd[2] =-1になってしまいます。
ε[1] = {0, I, 0, 0};
ε[2] = {0, 0, I, 0};
εd[1] = {0, I*Cos[theta], 0, I*Sin[theta]};
εd[2] = {0, 0, I, 0};
太字は３次元空間のベクトルなんじゃないの > >>69のpdf
>>太字は３次元空間のベクトルなんじゃないの > >>69のpdf
その通りです。解りました。
ε・ε＝－１
ε'・ε'＝－１
ε(1)・ε'(1)＝CosΘ
ε(1)・ε'(2)＝ε(2)・ε'(1)＝0
ε(2)・ε'(2)＝1←この式は、（－１）の間違いでは無いでしょうか？PDFのP14の一番上の式
そうしますと、以下のmathematica programで、計算できるのですが、如何でしょうか？
ep = (1/Sqrt[2])*{0, I, I, 0};
epdash = (I/Sqrt[2])*{0, -Cos[theta], 1, Sin[theta]};
Simplify[ep.ep]
Simplify[epdash.epdash]
(*分解*)
ep1 = {0, I, 0, 0};
ep2 = {0, 0, I, 0};
epdash1 = {0, -I*Cos[theta], 0, I*Sin[theta]};
epdash2 = {0, 0, I, 0};
(ep1.epdash1)^2 + (ep2.epdash2)^2 + (ep1.epdash2)^2 + (ep2.epdash1)^2
Simplify[ep1.ep1]
Simplify[ep2.ep2]
Simplify[ep1.epdash1]
Simplify[ep2.epdash2]
Simplify[ep1.epdash2]
Simplify[epdash2.ep1]
修正しました。これで如何でしょうか？
g = {{1, 0, 0, 0}, {0, -1, 0, 0}, {0, 0, -1, 0}, {0, 0, 0, -1}};
ep = (1/Sqrt[2])*{0, 1, 1, 0};
epdash = (1/Sqrt[2])*{0, -Cos[theta], 1, Sin[theta]};
Simplify[ep.g.ep]
Simplify[epdash.g.epdash]
(*分解*)
ep1 = {0, 1, 0, 0};
ep2 = {0, 0, 1, 0};
epdash1 = {0, -Cos[theta], 0, Sin[theta]};
epdash2 = {0, 0, 1, 0};
Simplify[(ep1.g.epdash1)^2 + (ep2.g.epdash2)^2 + (ep1.g.epdash2)^2 + \
(ep2.g.epdash1)^2]
Simplify[ep1.g.ep1]
Simplify[ep2.g.ep2]
Simplify[ep1.g.epdash1]
Simplify[ep2.g.epdash2]
Simplify[ep1.g.epdash2]
Simplify[epdash2.g.ep1]
Print["(*４元ベクトル*)"];
g4 = {{1, 0, 0, 0}, {0, -1, 0, 0}, {0, 0, -1, 0}, {0, 0, 0, -1}};
ep = (1/Sqrt[2])*{0, 1, 1, 0};
epdash = (1/Sqrt[2])*{0, Cos[theta], 1, Sin[theta]};
Simplify[ep.g4.ep]
Simplify[epdash.g4.epdash]
Print["(*３次元空間ベクトル分解*)"];
g3 = {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}};
ep1 = {1, 0, 0};
ep2 = {0, 1, 0};
epdash1 = {Cos[theta], 0, Sin[theta]};
epdash2 = {0, 1, 0};
Simplify[(ep1.g3.epdash1)^2 + (ep2.g3.epdash2)^2 + (ep1.g3.epdash2)^2 \
+ (ep2.g3.epdash1)^2]
Simplify[ep1.g3.ep1]
Simplify[ep2.g3.ep2]
Simplify[ep1.g3.epdash1]
Simplify[ep2.g3.epdash2]
Simplify[ep1.g3.epdash2]
Simplify[epdash2.g3.ep1]
epdash = (1/Sqrt[2])*{0, Cos[theta], 1, Sin[theta]};
としますと、kdashとepdashの内積が、"0"になりません。
下記のどこが間違っているでしょうか？
k0 = pi3;
k1 = 0;
k2 = 0;
k3 = -pi3;
kdash0 = pi3;
kdash1 = -pi3*Sqrt[1 - z^2];
kdash2 = 0;
kdash3 = -pi3*z;
episilon0 = 0;
episilon1 = (1/Sqrt[2])*1;
episilon2 = (1/Sqrt[2])*1;
episilon3 = 0;
(*この部分は正しいでしょうか？*)
episilondash0 = -1;
episilondash1 = Sqrt[1 - z^2];
episilondash2 = 1;
episilondash3 = z;
(*-------------------------*)
Simplify[episilondash0*episilondash0 - episilondash1*episilondash1 - episilondash2*episilondash2 - episilondash3*episilondash3]
Simplify[episilon0*episilon0 - episilon1*episilon1 - episilon2*episilon2 - episilon3*episilon3]
Simplify[kdash0*kdash0 - kdash1*kdash1 - kdash2*kdash2 - kdash3*kdash3]
Simplify[episilondash0*kdash0 - episilondash1*kdash1 - episilondash2*kdash2 - episilondash3*kdash3]
Simplify[episilon1*episilondash1 + episilon2*episilondash2 + episilon3*episilondash3]
この図の通り、計算したのですが、やはりダメです。
s = 2*k*pi + m^2;
u = -2*kdash*pi + m^2;
z = 1 + t/(2 *pi3^2);
t = 2 m^2 - s - u;
pi0 = (s + m^2)/(2 Sqrt[s]);
pi1 = 0;
pi2 = 0;
pi3 = (s - m^2)/(2 Sqrt[s]);
pf0 = pi0;
pf1 = pi3*Sqrt[1 - z^2];
pf2 = 0;
pf3 = pi3*z;
k0 = pi3;
k1 = 0;
k2 = 0;
k3 = -pi3;
kdash0 = pi3;
kdash1 = -pi3*Sqrt[1 - z^2];
kdash2 = 0;
kdash3 = -pi3*z;
episilon0 = 0;
episilon1 = (1/Sqrt[2])*1;
episilon2 = (1/Sqrt[2])*1;
episilon3 = 0;
(*この部分は正しいでしょうか？*)
episilondash0 = 0;
episilondash1 = (1/Sqrt[2])*Sqrt[1 - z^2];
episilondash2 = (1/Sqrt[2])*1;
episilondash3 = (1/Sqrt[2])*z;
(*----------------------*)
Print["(*この値は"０"になる*)"];
Print[Simplify[episilon0*k0 - episilon1*k1 - episilon2*k2 - episilon3*k3]];
Print[Simplify[episilondash0*kdash0 - episilondash1*kdash1 - episilondash2*kdash2 - episilondash3*kdash3]];
(*Print[Simplify[episilondash1*kdash1+episilondash2*kdash2+episilondash3*kdash3]];*)
Print[Simplify[episilon0*pi0 - episilon1*pi1 - episilon2*pi2 - episilon3*pi3]];
Print[Simplify[episilondash0*pi0 - episilondash1*pi1 - episilondash2*pi2 - episilondash3*pi3]];
Print[Simplify[k0*k0 - k1*k1 - k2*k2 - k3*k3]];
Print[Simplify[kdash0*kdash0 - kdash1*kdash1 - kdash2*kdash2 - kdash3*kdash3]];
Print["(*この値は"-1"になる*)"];
Print[Simplify[episilon0*episilon0 - episilon1*episilon1 - episilon2*episilon2 - episilon3*episilon3]];
Print[Simplify[episilondash0*episilondash0 - episilondash1*episilondash1 - episilondash2*episilondash2 - episilondash3*episilondash3]];
>kdash1 = -pi3*Sqrt[1 - z^2];
>kdash2 = 0;
>kdash3 = -pi3*z;
と
>episilondash1 = Sqrt[1 - z^2];
>episilondash2 = 1;
>episilondash3 = z;
が直交するはずだと思ってるの？
＞が直交するはずだと思ってるの？
その通りだと思い、計算していました。
下記HPの「mathematica」をクリック願います。
計算結果が表示します。
http://www43.tok2.com/home/iq188/
いま、PDFのP8の３行目
http://members3.jcom.home.ne.jp/nososnd/qed/comp.pdf
”p・? = p・?′ = 0（実験室系なので）になるので ”
という文から、そもそも以下（重心系）を、実験室系に変更しないと計算（ｐiε＝pfε'=0）が合わないのでは？
と考えてます。従いまして、現在
s = 2*k*pi + m^2;
u = -2*kdash*pi + m^2;
z = 1 + t/(2 *pi3^2);
t = 2 m^2 - s - u;
pi0 = (s + m^2)/(2 Sqrt[s]);
pi1 = 0;
pi2 = 0;
pi3 = (s - m^2)/(2 Sqrt[s]);
pf0 = pi0;←（この形は、重心系）
pf1 = pi3*Sqrt[1 - z^2];←（この形は、重心系）
pf2 = 0;
pf3 = pi3*z;←（この形は、重心系）
k0 = pi3;
k1 = 0;
k2 = 0;
k3 = -pi3;
kdash0 = pi3;
kdash1 = -pi3*Sqrt[1 - z^2];
kdash2 = 0;
kdash3 = -pi3*z;
http://osksn2.hep.sci.osaka-u.ac.jp/~naga/kogi/konan-class04/ch6-qed.pdf
のP4を元に、以前教えて頂いた重心系を、実験室系に変更したいのですが、解らず困っています。
”p・イプシロン = p・イプシロン′ = 0（実験室系なので）になるので ”
pi0 = (s + m^2)/(2 Sqrt[s]);←（この形は、重心系）
pi1 = 0;
pi2 = 0;
pi3 = (s - m^2)/(2 Sqrt[s]);←（この形は、重心系）
pf0 = pi0;
pf1 = pi3*Sqrt[1 - z^2];
pf2 = 0;
pf3 = pi3*z;
正しいでしょうか？
s =.;
u =.;
z =.;
t =.;
(*u=-2*k*q+m^2;
t=2 m^2-s-u;
z=1+t/(2 *k3^2);
w=1+u/(2 *k3^2);*)
k0 = (s - m^2)/ Sqrt[s];
k1 = 0;
k2 = 0;
k3 = -(s - m^2)/ Sqrt[s];
p0 = m;
p1 = 0;
p2 = 0;
p3 = 0;
j0 = (s - m^2)/(2 Sqrt[s]);
j1 = -(s - m^2)/(2 Sqrt[s])*Sqrt[1 - w^2];
j2 = 0;
j3 = (s - m^2)/(2 Sqrt[s])*w;
q0 = (s + m^2)/(2 Sqrt[s]);
q1 = (s - m^2)/(2 Sqrt[s])*Sqrt[1 - z^2];
q2 = 0;
q3 = (s - m^2)/(2 Sqrt[s])*z;
Print["(計算確認*)"];
p = p0*p0 - p1*p1 - p2*p2 - p3*p3;
q = q0*q0 - q1*q1 - q2*q2 - q3*q3;
k = k0*k0 - k1*k1 - k2*k2 - k3*k3;
j = j0*j0 - j1*j1 - j2*j2 - j3*j3;
Print["(*０になる*)"];
Simplify[(p + k)^2 - (j + q)^2]
Simplify[(p - q)^2 - (k - j)^2]
Simplify[(p - j)^2 - (q - k)^2]
Simplify[k]
Simplify[j]
Print["(*mになる*)"];
Simplify[p]
Simplify[q]
s = m^2;
Print["(*０になる*)"];
Simplify[p + k - j - q]
Simplify[k3 - (j3 + q3)](*z軸のつり合い*)
Simplify[j1 - q1](*ｙ軸のつり合い*)
sってどんな量？
http://members3.jcom.home.ne.jp/nososnd/field/mand.pdf
のP２の下から６行目の式です。
結局、
Simplify[j*q]
Simplify[p*k]
は、”０”ですから、ｓ＝ｍ＾２ではないでしょうか？
s=m^2 なら
>k0 = (s - m^2)/ Sqrt[s];
>k1 = 0;
>k2 = 0;
>k3 = -(s - m^2)/ Sqrt[s];
>j0 = (s - m^2)/(2 Sqrt[s]);
>j1 = -(s - m^2)/(2 Sqrt[s])*Sqrt[1 - w^2];
>j2 = 0;
>j3 = (s - m^2)/(2 Sqrt[s])*w;
は全部ゼロになるが、そういう計算をやりたいの？
＞は全部ゼロになるが、そういう計算をやりたいの？
違います。
すいません、入口が解らないです。
http://members3.jcom.home.ne.jp/nososnd/qed/comp.pdf
のP４の下から、８行目に、pi^2=pf^2=m^2,k^2=0とあります。
k^2= k0*k0 - k1*k1 - k2*k2 - k3*k3
は、実験室系でも、”0”になるのでしょうか？
最初は、
k0 = Sqrt[s] - m;
k1 = 0;
k2 = 0;
k3 = (Sqrt[s] - m);
p0 = m;
p1 = 0;
p2 = 0;
p3 = 0;
でよろしいでしょうか？
k1 = 0;
k2 = 0;
k3 = (Sqrt[s] - m);
p0 = m;
p1 = 0;
p2 = 0;
p3 = 0;
s = Simplify[2*(p*k) + m^2]
で、s=m^2になってしまいます。
すると、k0 = ｗとか、新しい変数を使う必要があるのでしょうか？
すると、実験室系では、重心系のように
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}]
とか、纏められるのでしょうか？
前スレより再録
>全然理解してないって昔から指摘されてるのに、全く進歩しないねぇ
>「ローレンツ不変量」とか知らんのだろう？
>特殊相対論を勉強するところから。
>デタラメやる前にちっとは勉強したら
>特殊相対論の勉強を頑に拒むのはなぜ？
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を見てやっと理解できました。
http://osksn2.hep.sci.osaka-u.ac.jp/~naga/kogi/konan-class04/ch6-qed.pdf
この式(6.26a)を、自力で導出するのは、かなり難しいと思います。
同様に、（知っている実験室系の式）の意味を理解するのは、出来ますが、導出するのは困難です。
従いまして、偏極を指定した断面積を自力で導出するのは、時間の無駄ですので、取りあえず明らめます。
http://www.kadokawa.co.jp/sp/201312-01/
下記HPでの非偏極での計算でOKとします。
http://amonphys.web.fc2.com/amonfc.pdf
三次関数の解を視覚的に表現してみたくて工夫しました。
切片を変数にしてその解を重ねてみました。
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}]
本来は 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複素数平面に
解が必ず３個有ります(重解含む)。
書き替えれば４次、５次関数でも表現できます。
係数には複素数も代入できます。
適当にいじってみていただけたら幸いです。
簡単だと思うよ。好みは人それぞれだけどね。
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"]
その通りですね。
ここで使用している数式は12/4(金)にNHKＥテレで放送された
数学ミステリー白熱教室で紹介されたものです。
ある数論の問題 y^2 + y == x^3 - x^2 ３次方程式の解で、
素数Pを法とする場合はある調和解析の関数となり、
複素数の場合は幾何学のトーラスとなるという話でした。
ラングランズ・プログラムというものらしいですが、
私には知識不足で理解できませんでした。
とりあえず、グラフの中央にある歪んだ円に大変興味を持ちました。
この辺が何かを意味しているような気がしますが、
私が考えられる次元はこの程度です。
どなたか知識のある方が何かを発見する手掛かりとなれたら本望です。
Maximaに書き換えるとき見た目を合わせ易いし
たまねぎのような構造
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[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}]
なんか面白くなってしまいました。
連投制限か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}]
即席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]
早速ノートブックにコピペしていじってみましたが、
当方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}]
H(えっち)-Functionは、
"Mathematica版「おっぱい（曲面）方程式」"
でググれば判る様な、無意味な遊びです。
ケルベロスと名付けたのは、犬みたいな形で顔の部分が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}]
こと山出しのニケア信教 ◆BkMobJaj1A。
その無知無教養さ加減はハゲの中でも下の下。
大阪大学や関西大学の学生からも馬鹿にされる為体。
St Marco Passion ◆BkMobJaj1A の臆病無比な卑劣さは、最底辺のハゲキモ蛆虫に劣る嫌らしさ！
こんな化け物老去勢オカマの存在を、我が国は赦しておいてはいけない。
早々にこの河馬ヅラをしたアナル変質者を朝鮮半島パクのもとへ強制送還させよう！！！
以下のように書いたらエラーが出てしまいました：
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}]
答えにはなっていませんが、ListのPart扱いなら、
a = {};
Do[AppendTo[a, 0], {10}](*初期値*)
SetSharedVariable[a];
ParallelDo[a[[i]] = i, {i, 1, 10}] (*並列計算*)
Do[Print[a[[i]]], {i, 1, 10}]
どうもです
私も結局リファレンス通りListで渡しました
ParallelDoでかなり高速化できました
しかし並列カーネル4個追加で\73,440って高いっす
挙句の果てに日本人に留学を勧めたり海外の大学を勧めたりするが
実際は外人なんてクズばかりで勉強しないくせに選民意識だけは高いし
人種差別は平気でするしでどうしようもない
留学なんて奨励してる奴はキチガイですわ
mathematicaの値段も頭おかしいのか？日本人が優秀だから日本でmathematicaを売ると
技術的にボロ負けしてしまうと危惧してるんだろうか？おそらくそうだろう。
アメリカがすごいのは大学や研究なんかじゃなくて金を荒稼ぎしてるからだからであって
日本が真似すると大やけどどころか沈没だよ。馬鹿な日本人ほどアメリカの真似したがるけど。
日本がすごいのは日本人が頑張ってきたからだよ。
図書館検索で全部借りればいい
https://www.kinokuniya.co.jp/f/dsg-01-9784000054478
おもしろそうです
証明と図形以外の問題で解けないのはないだろって
逆に言えばその二つが出来れば完璧なのに
おすすめ
『レクチャーズ オン Mathematica』
もくじなどのpdfがある
レクチャーズ オン Mathematica・サポートページ
http://www.math.titech.ac.jp/~kawahira/courses/mathematica.html
i=0～Nx
j=0～Ny
k=0～Nz
の範囲の整数として
{
{0,0,0},
{1,0,0},
{2,0,0},
… …,
{3,2,5},
…,
{Nx,Ny,Nz}
}
のような感じです．
Flatten[Table[{i, j, k}, {i, 0, 4}, {j, 0, 5}, {k, 0, 6}], 2]
とか？
ありがとうございましゅ！
○ありがとうございます！
電子の質量入れて計算すると飛びます。
分母のオーダーが10^(-31-α)とかでもだめですか？
具体例を
1/Exp[-x]
とするとxがどの辺のときに計算が放置されるのでしょうか？
x=10.^15*4 でアンダーフロー
電子の質量が入ってくるときとかどうするんだろうね…
その、電子の質量を入れて計算させると1/0エラーが出るっていう式をそのまま書いてみて
下のようにしたのですが
上手く行きません。何が悪いのでしょうか
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}]
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}]
こんなんみたいのかな
%→sol1とコンマ忘れでいけましたね
大変ありがとうございました
確率変数x,yとzがあり、それぞれ平均値mx,my,mzで標準偏差sx,sy,sz
の正規分布をもって変動するものとします。
このとき例えばz x/(x+y)の値がどのように分布するかを簡単にグラフ化
する手法と、その標準偏差を簡単に求める方法はあるでしょうか。
分布図はモンテカルロ法、統計値は一次近似の数式計算とかででしたが、
忘れた頃にやるので、火急の時には厳しい･･･
Ver10以下の範囲で、何かヒントがあればお願いします。
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[data]
data // StandardDeviation
という感じ。でもSDは安定しないな。やはり式の形が…
スムーズなグラフを描きたいのなら、
Plot[Evaluate@PDF[data // SmoothKernelDistribution, x], {x, -1, 5}]
とかやれば良い。
式自体は単純に、基準電圧zを抵抗y,xで分圧した電圧で、
その分布を求める場合ですが、モンカルだよりですかね。
変数毎に乱数発生をさせて、式に代入させていましたが、
いましたが、TransformedDistribution->RandomVariate
の流れは参考になりました。
あとSmoothKernelDistributionも(^^)/
eliminateって機能があって、
MarhematicaのEliminate関数と同じ、連立方程式からの変数消去なんだけど、
あまりにアホすぎて使い物にならんねコレ
個人的にMathematicaで一番重宝してる変数の整理、消去がMapleだと今市だな。
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だと回る．
なんで～～（泣）
ぐぬぬ…
割っただろ
ver 10.3 (Windows7_64bit) でも落ちる．
ver 8 (student ed., Windows XP pro) ではちゃんと計算される．
10.1と10.2はアップデートしてしまったから不明．
ver 10.3しか確かめてないけどMacProでも落ちるからハードが原因じゃない気がする．
はぁ～とりあえずver 10.0入ってる少し古い計算機使います．
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
今度は計算できる！不思議！
外部からアクセスできないキャッシュの内容によって挙動が変わるのか。
報告しとくかな。
↑Histogramでヒストグラムを描きました。
階級値の数字が300, 400, 500, 600と100間隔で表示されていますが、
300, 350, 400, 450, 500, 550, 600と表示させたいです。
どうすればいいのでしょうか？お願いします。
ex.横軸の目盛り範囲[100,800]隔50を指定する
sample = RandomVariate[NormalDistribution[500, 100], 1000];
Histogram[sample, {100, 800, 50},
Ticks -> {Range[100, 800, 50], Automatic}]
http://imgur.com/XMJ4pgw.jpg
ありがとうございました。
でも↑の画像のようになってしまいました。
数字は階級の両端の数字ではなく階級値を表示させたいんです。
ありがとうございました。
Ticks -> {Range[125, 775, 50], Automatic}]
みたいにやればできるんですね。
ありがとうございました。
その本良かったんですか。買ったけど読んでない。いまではmathematicaは
やめてmapleを使うことにした。
A.A (*これだと行列計算はしてくれない*)
これ一度 MarixFormでラッピングしてしまった後から、コード上でアンラップする方法って無いんですかね？
実用性云々(ラッピングしなきゃいいじゃん) ではなく単純に疑問です。
A = ToExpression @ StringSplit[#," "]& @ StringSplit[#,"\n\n"]& @ ToString[A]
文字列にして読み直ししたら一応いけたんですが、改行は2回スペースも2回続きそんなのをあてにしていいものやら...
もっと単純に一発でできてもよさそう。
A[[1]]
(A={{1,2},{3,4}})//MatrixForm
みたいにやるな。
>>164
「～Form の類いは、関数というよりはデータ構造構築子で、元データを1個だけ格納してる」 って理解で合ってるでしょうか
この辺りの詳細って WolframのWEBドキュメントに載ってます？
Part 関数を参照
A.A を計算するだけならMatrixFormでThreadもあり
A.A // Thread[#, MatrixForm] &
そのまんま「糸」でデータを縫い合わせる関数なのね。
基本的なのにまだまだ知らない関数があるなあ。
関西すうがく徒のつどい
関東すうがく徒のつどい
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
N[a, 10]
の出力が 1.23457 になってしまうのだけれど，
N関数って機械精度の数を渡されると常に6桁表示になるの？
InputForm[a]もしくはNumberForm[a,10]で表示された
フォントのサイズ変更したいから http://www.math.titech.ac.jp/~kawahira/courses/mathematica/01.pdf にある
> メニューから「編集」→「環境設定...」* → 「詳細」→ 「オプションインスペクタを開く」→ 「書式設定」→ 「フォント設定」→ 「FontSize(サイズ)」とすすみ,これをデフォルトの 12 から適宜増や す.「適用」ボタンを押すとフォントサイズが大きくなる.
を行おうとしたが、「適用」ボタンがグレーのままで押しても適用出来ない。
どうしたら良いんだコレ？
http://i.imgur.com/BEHcB9D.png
{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]
と記述してるんですが，より高速な方法があれば教えてください
result = Transpose@list1;
result = AppendTo[result, list2];
result = Transpose@result;
で速度が2倍ぐらいになりました
MapThread[Append[#1, #2] &, {list1, list2}]
かなあ。
更に速くなりました
記述もスマートですしこれを使いたいと思います
ありがとうございました
↓みたいにエラーが発生するんだけど同じひといる？
1つずつ閉じても最後の1個を閉じるところで同じエラーがでる
http://i.imgur.com/kCogne9.png
>1 ：名無しさん ：2006/04/30(日) 01:41:01 ID:KPnB.CH2
> 迷惑かしらん
>
>5568 ：名無しさん：2016/08/17(水) 18:26:13 ID:???
> うるさい
>
>5569 ：kmath1107★：2016/08/17(水) 21:46:32 ID:???
> 人への念の盗み見による介入を阻め。
>
>5571 ：名無しさん：2016/08/17(水) 23:39:07 ID:???
> うるさい
>
>5576 ：kmath1107★ ：2016/08/18(木) 20:58:14 ID:???
> 人への念の盗み見による介入を阻め。
>
>5577 ：名無しさん ：2016/08/18(木) 21:05:02 ID:???
> >>5575
> うるさい
>
> >>5576
> 賛同致します
>
>5578 ：kmath1107★ ：2016/08/19(金) 08:46:22 ID:???
> 人への念の盗み見による介入を阻め。
>
> Re:>>5577 人への念の盗み見による介入が無くなれば世の不和が無くなるだろう.
>
>5582 ：￥ ◆2VB8wsVUoo ：2016/08/19(金) 08:53:36 ID:???
> 芳雄が理想とし、自ら体現する大学教授とは？
> ０．自分が『お教授である』という利点を徹底活用して、偉そうに振舞う。
> １．年寄りや権威には擦り寄って顔色を窺い、ラクして損しない様にスル。
> ２．難しい分野や困難な研究テーマは徹底して避けて、努力を最小化する。
> ３．高い学歴とか権威を効率的に利用して、自分を飾って偉く見せ掛ける。
> ４．他人に見える部分だけを巧みに繕ってメッキし、人格者のフリをする。
> ５．相手のオツムの質を窺い、シッタカだけで見識がある様に見せ掛ける。
> ６．自分よりも優秀な人間は絶対に敵に回さないでヘラヘラと仲良くする。
> ７．自分から見てダメオツムな野郎は、上から目線で威圧して屈服させる。
> ８．大して中身が無いカラッポ知識を針小棒大に騒ぎ立て、蘊蓄を傾ける。
> ９．自分の大脳が働いてない低能ぶりは、口先で適当に誤魔化して逃げる。
>
> ￥
>
>>151と>>192治ってるね
なんとか安定したような気がする。
>1 ：名無しさん ：2006/04/30(日) 01:41:01 ID:KPnB.CH2
> 迷惑かしらん
>
>5568 ：名無しさん：2016/08/17(水) 18:26:13 ID:???
> うるさい
>
>5569 ：kmath1107★：2016/08/17(水) 21:46:32 ID:???
> 人への念の盗み見による介入を阻め。
>
>5571 ：名無しさん：2016/08/17(水) 23:39:07 ID:???
> うるさい
>
>5576 ：kmath1107★ ：2016/08/18(木) 20:58:14 ID:???
> 人への念の盗み見による介入を阻め。
>
>5577 ：名無しさん ：2016/08/18(木) 21:05:02 ID:???
> >>5575
> うるさい
>
> >>5576
> 賛同致します
>
>5578 ：kmath1107★ ：2016/08/19(金) 08:46:22 ID:???
> 人への念の盗み見による介入を阻め。
>
> Re:>>5577 人への念の盗み見による介入が無くなれば世の不和が無くなるだろう.
>
>5582 ：￥ ◆2VB8wsVUoo ：2016/08/19(金) 08:53:36 ID:???
> 芳雄が理想とし、自ら体現する大学教授とは？
> ０．自分が『お教授である』という利点を徹底活用して、偉そうに振舞う。
> １．年寄りや権威には擦り寄って顔色を窺い、ラクして損しない様にスル。
> ２．難しい分野や困難な研究テーマは徹底して避けて、努力を最小化する。
> ３．高い学歴とか権威を効率的に利用して、自分を飾って偉く見せ掛ける。
> ４．他人に見える部分だけを巧みに繕ってメッキし、人格者のフリをする。
> ５．相手のオツムの質を窺い、シッタカだけで見識がある様に見せ掛ける。
> ６．自分よりも優秀な人間は絶対に敵に回さないでヘラヘラと仲良くする。
> ７．自分から見てダメオツムな野郎は、上から目線で威圧して屈服させる。
> ８．大して中身が無いカラッポ知識を針小棒大に騒ぎ立て、蘊蓄を傾ける。
> ９．自分の大脳が働いてない低能ぶりは、口先で適当に誤魔化して逃げる。
>
> ￥
>
Plot[f[x], {x, 0, 0.001}]
とすると、グラフがぺったんこになってしまいます。
これをちゃんと描くにはどうしたらいいのでしょうか？
>1 ：名無しさん ：2006/04/30(日) 01:41:01 ID:KPnB.CH2
> 迷惑かしらん
>
>5568 ：名無しさん：2016/08/17(水) 18:26:13 ID:???
> うるさい
>
>5569 ：kmath1107★：2016/08/17(水) 21:46:32 ID:???
> 人への念の盗み見による介入を阻め。
>
>5571 ：名無しさん：2016/08/17(水) 23:39:07 ID:???
> うるさい
>
>5576 ：kmath1107★ ：2016/08/18(木) 20:58:14 ID:???
> 人への念の盗み見による介入を阻め。
>
>5577 ：名無しさん ：2016/08/18(木) 21:05:02 ID:???
> >>5575
> うるさい
>
> >>5576
> 賛同致します
>
>5578 ：kmath1107★ ：2016/08/19(金) 08:46:22 ID:???
> 人への念の盗み見による介入を阻め。
>
> Re:>>5577 人への念の盗み見による介入が無くなれば世の不和が無くなるだろう.
>
>5582 ：￥ ◆2VB8wsVUoo ：2016/08/19(金) 08:53:36 ID:???
> 芳雄が理想とし、自ら体現する大学教授とは？
> ０．自分が『お教授である』という利点を徹底活用して、偉そうに振舞う。
> １．年寄りや権威には擦り寄って顔色を窺い、ラクして損しない様にスル。
> ２．難しい分野や困難な研究テーマは徹底して避けて、努力を最小化する。
> ３．高い学歴とか権威を効率的に利用して、自分を飾って偉く見せ掛ける。
> ４．他人に見える部分だけを巧みに繕ってメッキし、人格者のフリをする。
> ５．相手のオツムの質を窺い、シッタカだけで見識がある様に見せ掛ける。
> ６．自分よりも優秀な人間は絶対に敵に回さないでヘラヘラと仲良くする。
> ７．自分から見てダメオツムな野郎は、上から目線で威圧して屈服させる。
> ８．大して中身が無いカラッポ知識を針小棒大に騒ぎ立て、蘊蓄を傾ける。
> ９．自分の大脳が働いてない低能ぶりは、口先で適当に誤魔化して逃げる。
>
> ￥
>
入ってた。便利な時代になったなあ。
でも結局無料マイククラフトで遊んでばっかりw
次回、開いたときに、再評価しないといけなくて、面倒なのですが、
みなさんはどうやって計算した結果などを保存していますか？
あと、定義した関数などをマイライブラリにする方法を教えてください。
計算結果格納した変数を用意してmx形式で Export するとか．
次回はそれを Import する．
とjpg画像として保存しようとすると
Export::noopen : Cannot open fig1.jpg.
というエラーが出てしまいます。
どうすれば、いいのでしょうか？
>1 ：名無しさん ：2006/04/30(日) 01:41:01 ID:KPnB.CH2
> 迷惑かしらん
>
>5714 ：名無しさん：2016/09/01(木) 22:40:59 ID:???
> >>5711
> 黙ってろカスが。お前こんなことずっと続けてて父親に申し訳ないと思わないのか。
>
>5718 ：￥ ◆2VB8wsVUoo ：2016/09/02(金) 07:47:45 ID:???
> >>5714
> マジレスしておくが、芳雄にはきちんと罰だけは受けて貰う。あんな酷い事をし
> ておきながら、無傷であの世に逃亡というのは絶対に許されない。死ぬ前に充分
> な精神的苦痛をタップリと味わうべき。あの糞野郎だけは絶対に許されないので。
> 芳雄に対する怒りと憎しみは、馬鹿板に対する怒りとは比べ物にはならんわ。
>
> ￥
>
>5720 ：￥ ◆2VB8wsVUoo ：2016/09/02(金) 08:54:09 ID:???
> >>5714
> 言って置くが、被害を受け始めた高校生の頃から私は芳雄を論理分析し、その欠
> 陥や弱点を精密に理解し、そしてその横暴極まりない無責任な態度に対抗しなが
> ら狙い撃ちにして来た。私は芳雄のせいで甚大な被害を被ったのであり、それを
> 「親が責任を取る」という様ないい加減な逃げ口上で逃亡し、無責任を通す卑怯
> な行為は到底許されない。なのでその報いだけでもきちんと受けさせてやるだけ。
>
> 糞芳雄の野郎、このまま逃げ切りは許さない。自分から言い放った『親としての
> 責任』というものが微塵でも残ってるのであれば、それ相当の行為が自らなされ
> て当然というものだろう。手を切り落とすもよし、足を切り落とすもよし。或い
> は自分で主張した釜ヶ崎に自分で行って、そして労務者にでも殴られて撲殺され
> るのもいいだろう。
>
> とにかく自分で言った事だけは、きちんと自分から実行するべき。知らぬ存ぜぬ
> で、無責任な逃げ切りだけは許されない。ソレこそが芳雄が言う所の卑怯者だか
> らだ。糞芳雄は恥を知るべき。今からでもいいから、尊厳の意味を理解するべき。
>
> ￥
>
でディレクトリを設定したら保存できました。
デフォルトではディレクトリが "/" になっていました。
これを自動的にノートブックが保存されているディレクトリにするにはどうすればいいのでしょうか？
>1 ：名無しさん ：2006/04/30(日) 01:41:01 ID:KPnB.CH2
> 迷惑かしらん
>
>5714 ：名無しさん：2016/09/01(木) 22:40:59 ID:???
> >>5711
> 黙ってろカスが。お前こんなことずっと続けてて父親に申し訳ないと思わないのか。
>
>5718 ：￥ ◆2VB8wsVUoo ：2016/09/02(金) 07:47:45 ID:???
> >>5714
> マジレスしておくが、芳雄にはきちんと罰だけは受けて貰う。あんな酷い事をし
> ておきながら、無傷であの世に逃亡というのは絶対に許されない。死ぬ前に充分
> な精神的苦痛をタップリと味わうべき。あの糞野郎だけは絶対に許されないので。
> 芳雄に対する怒りと憎しみは、馬鹿板に対する怒りとは比べ物にはならんわ。
>
> ￥
>
>5720 ：￥ ◆2VB8wsVUoo ：2016/09/02(金) 08:54:09 ID:???
> >>5714
> 言って置くが、被害を受け始めた高校生の頃から私は芳雄を論理分析し、その欠
> 陥や弱点を精密に理解し、そしてその横暴極まりない無責任な態度に対抗しなが
> ら狙い撃ちにして来た。私は芳雄のせいで甚大な被害を被ったのであり、それを
> 「親が責任を取る」という様ないい加減な逃げ口上で逃亡し、無責任を通す卑怯
> な行為は到底許されない。なのでその報いだけでもきちんと受けさせてやるだけ。
>
> 糞芳雄の野郎、このまま逃げ切りは許さない。自分から言い放った『親としての
> 責任』というものが微塵でも残ってるのであれば、それ相当の行為が自らなされ
> て当然というものだろう。手を切り落とすもよし、足を切り落とすもよし。或い
> は自分で主張した釜ヶ崎に自分で行って、そして労務者にでも殴られて撲殺され
> るのもいいだろう。
>
> とにかく自分で言った事だけは、きちんと自分から実行するべき。知らぬ存ぜぬ
> で、無責任な逃げ切りだけは許されない。ソレこそが芳雄が言う所の卑怯者だか
> らだ。糞芳雄は恥を知るべき。今からでもいいから、尊厳の意味を理解するべき。
>
> ￥
>
俺は毎回
SetDirectory@NotebookDirectory[];
って書いてるけどこれが面倒なの？
は、
1/2 + sqrt(3)/2*i
1/2 - sqrt(3)/2*i
のどちらを表わすのでしょうか？
何か決まりがあるのでしょうか？
いちいち、実部と虚部を取って確認するのが面倒です。
>414 名前：１３２人目の素数さん ：2016/09/15(木) 12:12:20.88 ID:zW1tnm+Q
> 痴漢でもなんでも、業績を上げていれば世間はほっておかない。
> 哲也がヤケッパチになったのは、自分がただの有象無象の一人に過ぎないことを思い知ったから。
> 痴漢行為はヤケッパチになった結果としての所業だろう。
>
なら (-1)^(1/3) は Exp[iπ/3] でしょ？
>57 :１３２人目の素数さん 2016/09/18(日) 22:27:13.98 ID:ol2qzlaN
>東京の税金たからないで関西に帰れよ部落民
>
>61 :１３２人目の素数さん 2016/09/18(日) 22:54:42.02 ID:ol2qzlaN
>お前の親でもない赤の他人の東京都からナマポふんだくって一人前ズラしてんじゃねえよ関西部落民
>とっとと痴漢しに関西に帰れよ
>鬼のようにぶっさいくな関東女茨城女よりも四国や関西で痴漢した方が楽しいんだろ？ん？
>64 :１３２人目の素数さん 2016/09/18(日) 23:04:21.99 ID:ol2qzlaN
>関西に帰るのがいやならフランスの旧植民地の大学にでも行きゃあいいだろうがよ
>フランス崇め奉ってる植民地根性なら旧植民地がお似合いだろうがよ
>フランス文化受容強要の美しい現実も見られるだろうしなｗ
>
>66 :１３２人目の素数さん 2016/09/18(日) 23:10:41.77 ID:ol2qzlaN
>親のいいなりの関西のお受験坊ちゃんなんて関東女に相手されるどころかキンタマ握りつぶされそうだしな
>フランス人女には口げんかで負けるし
>フランス植民地選良の黒人女にでも慰めてもらえば？
>
>70 :１３２人目の素数さん 2016/09/18(日) 23:16:27.29 ID:ol2qzlaN
>フランス語でフランスの掲示板で関西のお受験パパの暴虐っぷりとお受験坊やのへなちょこっぷり宣伝して回ればいいのに
>日本人には常識レベルの関西のテストだけ誤魔化すのに特化したお受験親子の実力皆無っぷりなんて自分で実証して回ったって冷笑すらされんだろうに？
>旧関係者のガン無視っぷりにますます逆恨みがつのるだけじゃねぇのｗｗｗｗ
>
>72 :１３２人目の素数さん 2016/09/18(日) 23:20:13.28 ID:ol2qzlaN
>実力ありゃあ旧関係者も可哀想に感じてどっかにアカデミックな引き立てもあらアするかも知らねえが
>宇沢親子をしょっぱくしたコネと実績じゃあ冷笑すらされずガン無視空気だわな日本に限らず
>
>74 :１３２人目の素数さん 2016/09/18(日) 23:22:55.17 ID:ol2qzlaN
>親が教授じゃなかったらアカポスに最初っから引っかからなかった程度の癖に上等だな関西のお受験お坊ちゃん
>
>57 :１３２人目の素数さん 2016/09/18(日) 22:27:13.98 ID:ol2qzlaN
>東京の税金たからないで関西に帰れよ部落民
>
>61 :１３２人目の素数さん 2016/09/18(日) 22:54:42.02 ID:ol2qzlaN
>お前の親でもない赤の他人の東京都からナマポふんだくって一人前ズラしてんじゃねえよ関西部落民
>とっとと痴漢しに関西に帰れよ
>鬼のようにぶっさいくな関東女茨城女よりも四国や関西で痴漢した方が楽しいんだろ？ん？
>64 :１３２人目の素数さん 2016/09/18(日) 23:04:21.99 ID:ol2qzlaN
>関西に帰るのがいやならフランスの旧植民地の大学にでも行きゃあいいだろうがよ
>フランス崇め奉ってる植民地根性なら旧植民地がお似合いだろうがよ
>フランス文化受容強要の美しい現実も見られるだろうしなｗ
>
>66 :１３２人目の素数さん 2016/09/18(日) 23:10:41.77 ID:ol2qzlaN
>親のいいなりの関西のお受験坊ちゃんなんて関東女に相手されるどころかキンタマ握りつぶされそうだしな
>フランス人女には口げんかで負けるし
>フランス植民地選良の黒人女にでも慰めてもらえば？
>
>70 :１３２人目の素数さん 2016/09/18(日) 23:16:27.29 ID:ol2qzlaN
>フランス語でフランスの掲示板で関西のお受験パパの暴虐っぷりとお受験坊やのへなちょこっぷり宣伝して回ればいいのに
>日本人には常識レベルの関西のテストだけ誤魔化すのに特化したお受験親子の実力皆無っぷりなんて自分で実証して回ったって冷笑すらされんだろうに？
>旧関係者のガン無視っぷりにますます逆恨みがつのるだけじゃねぇのｗｗｗｗ
>
>72 :１３２人目の素数さん 2016/09/18(日) 23:20:13.28 ID:ol2qzlaN
>実力ありゃあ旧関係者も可哀想に感じてどっかにアカデミックな引き立てもあらアするかも知らねえが
>宇沢親子をしょっぱくしたコネと実績じゃあ冷笑すらされずガン無視空気だわな日本に限らず
>
>74 :１３２人目の素数さん 2016/09/18(日) 23:22:55.17 ID:ol2qzlaN
>親が教授じゃなかったらアカポスに最初っから引っかからなかった程度の癖に上等だな関西のお受験お坊ちゃん
>
>329 名前：１３２人目の素数さん ：2016/10/26(水) 17:24:53.77 ID:LOz3feeR
> 孫を奪うと言いつつ好みのタイプでムラムラ牌揉み
> 情けないねえ...
>
>331 名前：ニノ ：2016/10/28(金) 04:27:43.07 ID:YMecsJRq
> 女性に痴漢は最低の行為です
> 恥を知ってください
>
ヘルプをみろ
>326 名前：１３２人目の素数さん ：2016/10/26(水) 14:13:06.88 ID:dz3LPMS/
> >>2-3
> 虐待的躾
>
> >>4
> その主張は個人の精神の内に限られた中での狭い正論
> 法的には不当な手続きによる風評悪化工作
>
> 縁組解脱するでもなければ復讐も果たせず、無駄以上のマイナスに終わる人生へ進む増田哲也
>
>327 名前：１３２人目の素数さん ：2016/10/26(水) 14:15:31.32 ID:dz3LPMS/
> 「孫を奪う」
> 手段が
> 「自分の任意を駆使して｣
> …
>
> 警察・検察・弁護士・裁判官「それを仄めかすと言う｣
>
>328 名前：１３２人目の素数さん ：2016/10/26(水) 14:29:49.76 ID:dz3LPMS/
> しかし最後に決めるのは法律ではなく過去判例や世情判例、法律は実は決め切らない
> 口八丁手八丁で判例が左右する
>
> 増田哲也の不毛な戦いは続く
>
>330 名前：１３２人目の素数さん ：2016/10/26(水) 17:54:03.72 ID:dz3LPMS/
> 「ムラムラしたんだから仕方ないがな」みたいな発言
>
> 性感を、絶頂に至らせずに高め強める為に
> ある周波数かつ弱電流高電圧で絶頂越えを感電制止し
> 無限に高まり強まる性感を与えつつ
> 性感を全身に広げる為に口腔鼻腔～肺の内気を
> 窒息感たっぷりに高圧に保ちつつも酸欠しない様に内気を常に酸素還元する
>
> 無限に高まり強まりつつ全身に至る性感の中で
> 余りにも強い性感に死に対する欲望が怒涛の勢いで湧き上がり
> 性感の大海の内に命が尽きる事を願うだろう
> 義父に対する復讐よりも死に焦がれる性欲に堕ちるのだ
> 「ムラムラ」「仕方ない」のだから
>
v11.0.0もv11.0.1もどっちも遅いことを確認した。
比較対象はv10の最後のバージョンのv10.4.1とする。
＝＝＝＝＝ 激遅ぷんぷん丸確認コード ＝＝＝＝＝
まずは適当な画像ファイルを用意する。
ここでは以下のスパイキーとする。
http://www.wolfram.com/homepage/img/header-spikey.png
確認コードはこれを10回インポートするだけ。
a = Table[Import[
2chもバグか？
途中で消えるぞ。
v11にしたんだがImportがくっそ遅くなっててブチ切れ中ｗ
v11.0.0もv11.0.1もどっちも遅いことを確認した。
比較対象はv10の最後のバージョンのv10.4.1とする。
＝＝＝＝＝ 激遅ぷんぷん丸確認コード ＝＝＝＝＝
まずは適当な画像ファイルを用意する。
ここでは以下のスパイキーとする。
http://www.wolfram.com/homepage/img/header-spikey.png
確認コードはこれを10回インポートするだけ。
a = Table[Import["header-spikey.png"], {i, 10}]; // AbsoluteTiming
これが、v10では0.5秒なのに対して、v11では65秒かかる。
ウルフラムって製品リリース前にしっかりテストしないんだよな～。
ユーザーをデバッガーだと勘違いしてる。
ふざけんなよ。
> なんだ？
> 2chもバグか？
> 途中で消えるぞ。
確認画面が出る初回書き込みの時、"でコメントが切れるバグがある。
多分、通常書き込みと違い、レスのエスケープを省略してるんだろう（2ch運営はド阿呆だから）
>326 名前：１３２人目の素数さん ：2016/10/26(水) 14:13:06.88 ID:dz3LPMS/
> >>2-3
> 虐待的躾
>
> >>4
> その主張は個人の精神の内に限られた中での狭い正論
> 法的には不当な手続きによる風評悪化工作
>
> 縁組解脱するでもなければ復讐も果たせず、無駄以上のマイナスに終わる人生へ進む増田哲也
>
>327 名前：１３２人目の素数さん ：2016/10/26(水) 14:15:31.32 ID:dz3LPMS/
> 「孫を奪う」
> 手段が
> 「自分の任意を駆使して｣
> …
>
> 警察・検察・弁護士・裁判官「それを仄めかすと言う｣
>
>328 名前：１３２人目の素数さん ：2016/10/26(水) 14:29:49.76 ID:dz3LPMS/
> しかし最後に決めるのは法律ではなく過去判例や世情判例、法律は実は決め切らない
> 口八丁手八丁で判例が左右する
>
> 増田哲也の不毛な戦いは続く
>
>330 名前：１３２人目の素数さん ：2016/10/26(水) 17:54:03.72 ID:dz3LPMS/
> 「ムラムラしたんだから仕方ないがな」みたいな発言
>
> 性感を、絶頂に至らせずに高め強める為に
> ある周波数かつ弱電流高電圧で絶頂越えを感電制止し
> 無限に高まり強まる性感を与えつつ
> 性感を全身に広げる為に口腔鼻腔～肺の内気を
> 窒息感たっぷりに高圧に保ちつつも酸欠しない様に内気を常に酸素還元する
>
> 無限に高まり強まりつつ全身に至る性感の中で
> 余りにも強い性感に死に対する欲望が怒涛の勢いで湧き上がり
> 性感の大海の内に命が尽きる事を願うだろう
> 義父に対する復讐よりも死に焦がれる性欲に堕ちるのだ
> 「ムラムラ」「仕方ない」のだから
>
SQLインジェクションでもされそうだな
お願いします
Wolfram Alphaの書き方について質問です。グレブナー基底を計算したいのですが、ほぼ同じ文構造にしているのに、片方だけ門前払いのようになってしまいます。... #知恵袋_
http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q10167091570
>前科持ち変質者と絶対出会える掲示板 [無断転載禁止]
>
>1 名前：１３２人目の素数さん 2016/11/16(水) 21:02:24.40 ID:8UX5OsVV
> 変質者前科持ちと気が触れ合える掲示板
>
>11 名前：１３２人目の素数さん ：2016/11/19(土) 08:36:12.59 ID:6KwDBI7h
> 変質者前科持ち＝増田哲也
>
>12 名前：１３２人目の素数さん ：2016/11/19(土) 09:04:39.15 ID:AZB04dZ8
> わざわざ言わんでもええ
>
>13 名前：出会える掲示板 ◆2VB8wsVUoo ：2016/11/19(土) 15:58:01.20 ID:21LrO2+x
> 絶対に…
>
> ケケケ￥
>
>14 名前：１３２人目の素数さん ：2016/11/19(土) 16:31:33.55 ID:6KwDBI7h
> 六十目前で父親逆恨みしたり掲示板逆恨みする根性の腐れっぷりは凄くて困る
>
画像小さすぎ
http://i.imgur.com/4VaTUrE.jpg
http://i.imgur.com/dFvo35f.jpg
黄緑のマーカーの部分しか違わないのになぜですかね
>前科持ち変質者と絶対出会える掲示板 [無断転載禁止]
>
>1 名前：１３２人目の素数さん 2016/11/16(水) 21:02:24.40 ID:8UX5OsVV
> 変質者前科持ちと気が触れ合える掲示板
>
>11 名前：１３２人目の素数さん ：2016/11/19(土) 08:36:12.59 ID:6KwDBI7h
> 変質者前科持ち＝増田哲也
>
>12 名前：１３２人目の素数さん ：2016/11/19(土) 09:04:39.15 ID:AZB04dZ8
> わざわざ言わんでもええ
>
>13 名前：出会える掲示板 ◆2VB8wsVUoo ：2016/11/19(土) 15:58:01.20 ID:21LrO2+x
> 絶対に…
>
> ケケケ￥
>
>14 名前：１３２人目の素数さん ：2016/11/19(土) 16:31:33.55 ID:6KwDBI7h
> 六十目前で父親逆恨みしたり掲示板逆恨みする根性の腐れっぷりは凄くて困る
>
長過ぎて入力できてないんじゃないの
>前科持ち変質者と絶対出会える掲示板 [無断転載禁止]
>
>1 名前：１３２人目の素数さん 2016/11/16(水) 21:02:24.40 ID:8UX5OsVV
> 変質者前科持ちと気が触れ合える掲示板
>
>11 名前：１３２人目の素数さん ：2016/11/19(土) 08:36:12.59 ID:6KwDBI7h
> 変質者前科持ち＝増田哲也
>
>12 名前：１３２人目の素数さん ：2016/11/19(土) 09:04:39.15 ID:AZB04dZ8
> わざわざ言わんでもええ
>
>13 名前：出会える掲示板 ◆2VB8wsVUoo ：2016/11/19(土) 15:58:01.20 ID:21LrO2+x
> 絶対に…
>
> ケケケ￥
>
>14 名前：１３２人目の素数さん ：2016/11/19(土) 16:31:33.55 ID:6KwDBI7h
> 六十目前で父親逆恨みしたり掲示板逆恨みする根性の腐れっぷりは凄くて困る
>
成功したコードと失敗したコードをここに書いてみ。
>前科持ち変質者と絶対出会える掲示板 [無断転載禁止]
>
>1 名前：１３２人目の素数さん 2016/11/16(水) 21:02:24.40 ID:8UX5OsVV
> 変質者前科持ちと気が触れ合える掲示板
>
>11 名前：１３２人目の素数さん ：2016/11/19(土) 08:36:12.59 ID:6KwDBI7h
> 変質者前科持ち＝増田哲也
>
>12 名前：１３２人目の素数さん ：2016/11/19(土) 09:04:39.15 ID:AZB04dZ8
> わざわざ言わんでもええ
>
>13 名前：出会える掲示板 ◆2VB8wsVUoo ：2016/11/19(土) 15:58:01.20 ID:21LrO2+x
> 絶対に…
>
> ケケケ￥
>
>14 名前：１３２人目の素数さん ：2016/11/19(土) 16:31:33.55 ID:6KwDBI7h
> 六十目前で父親逆恨みしたり掲示板逆恨みする根性の腐れっぷりは凄くて困る
>
下の画像はメモ帳のキャプチャでコード書いてあるメモ帳をコピペして計算したから下の画像のとうりです
>前科持ち変質者と絶対出会える掲示板 [無断転載禁止]
>
>1 名前：１３２人目の素数さん 2016/11/16(水) 21:02:24.40 ID:8UX5OsVV
> 変質者前科持ちと気が触れ合える掲示板
>
>11 名前：１３２人目の素数さん ：2016/11/19(土) 08:36:12.59 ID:6KwDBI7h
> 変質者前科持ち＝増田哲也
>
>12 名前：１３２人目の素数さん ：2016/11/19(土) 09:04:39.15 ID:AZB04dZ8
> わざわざ言わんでもええ
>
>13 名前：出会える掲示板 ◆2VB8wsVUoo ：2016/11/19(土) 15:58:01.20 ID:21LrO2+x
> 絶対に…
>
> ケケケ￥
>
>14 名前：１３２人目の素数さん ：2016/11/19(土) 16:31:33.55 ID:6KwDBI7h
> 六十目前で父親逆恨みしたり掲示板逆恨みする根性の腐れっぷりは凄くて困る
>
質問ならテキストで貼れよ
アホなのかな？
>>373
そういうことでしたかすみません
やはり長すぎて入力できていなかったみたいなのですが、入力する方法はお知りでしょうか？
ゆうりょうばんをかうしかないのでしょうか
すみません
変数をxyzwに置き換えたところ、文字数以内に収まりました
ご迷惑をおかけしました
>前科持ち変質者と絶対出会える掲示板 [無断転載禁止]
>
>1 名前：１３２人目の素数さん 2016/11/16(水) 21:02:24.40 ID:8UX5OsVV
> 変質者前科持ちと気が触れ合える掲示板
>
>11 名前：１３２人目の素数さん ：2016/11/19(土) 08:36:12.59 ID:6KwDBI7h
> 変質者前科持ち＝増田哲也
>
>12 名前：１３２人目の素数さん ：2016/11/19(土) 09:04:39.15 ID:AZB04dZ8
> わざわざ言わんでもええ
>
>13 名前：出会える掲示板 ◆2VB8wsVUoo ：2016/11/19(土) 15:58:01.20 ID:21LrO2+x
> 絶対に…
>
> ケケケ￥
>
>14 名前：１３２人目の素数さん ：2016/11/19(土) 16:31:33.55 ID:6KwDBI7h
> 六十目前で父親逆恨みしたり掲示板逆恨みする根性の腐れっぷりは凄くて困る
>
それとも、 Mathematica の基本的な機能は10年以上変わっていませんか？
いまでは手に入りにくい本も多い
ありがとうございます。
見てみようと思います。
ありがとうございます。
そういう「工夫」をしなければならないんですか。
コマンド一つではできませんか？
確かにコマンド一つですね。
でも、高次の項を除去するのに、わざわざ多項式の余りを計算する
という無駄なことをやっているため気持ちが悪いです。
正統的な方法で高次の項を除去しようと思うなら、プログラムを組む
しかないということでしょうか？
グラフ理論の最短路問題を試してみたら、
HighlightGraph[] で最短経路をハイライトできませんでした。
ノードは赤くなるのですが、エッジの色や太さなどが変更できません。
ノードの番号やエッジの番号もマニュアルに載っている方法でやっても表示されなかったです。
別の指定の仕方をすると表示されましたが。
MathematicaとMaple、どっちのほうがいいんですか？
Mapleってアップグレードがないっぽいんですけど、
最新版が出るたびにたけえ金払えってことですか？
メープルは普通のプログラミング言語っぽい。
マセマティカは関数型だからちょっと変わっている。
でもユーザーが多そうなマセマティカかなー。選ぶなら。
でも，Mapleにはアップグレードがないっていうのが気になるんですが…
最新版が出るたびに高い金を払え！ってことですよね？
Table[Coefficient[p[x],x^k] x^k,{k,0,n-1}].Table[x^k,{k,1,n-1}]+p[0]
うまいコマンドがあればいいね。
あと、
Plot3D で描いた曲面 z = f(x, y) 上の曲線 x = g1(t), y = g2(t) を「綺麗に」描く方法を教えてください。
半透明にするにはどうすればいいのでしょうか？
PlotStyleを指定できないようです。
今、たまたまドキュメントを見ていたら、こんな例題がありました：
Normal[(1 + 2 x)^1000 + O[x]^5]
これで、 x^5 以上の項を削除できます。
マセマティカでの O[x]^n とは何を意味するのでしょうか？
ありがとうございます。
本当ですか？
なんか違うような気がするのですが。
Sin[x^2 + x^3] + O[x]^10
↑は、 x^9までの項を求めるという意味ですよね？
とドキュメントに書いてありました。
でもそういう意味だけじゃなくて、
>>433
のように x^9 の項まで求めよと指示するような役割もありますよね。
より一般的な O(x^n) ではないことに注意せよ。
と書かれていますが、 O(x)^n の数学における定義ってなんですか？
ようするに程度問題を表している。
正確には
O[x]^n/x^n-> const as x->0
普通は、面倒だから４２９のように（O[x]^nー＞０，ｘ－＞０）のニュアンスをこめて、
書く人が多いが、やはり定数のほうがいいという人が多い。
（何も書かない人もおおい）
o[x]^n/x^n-> 0 as x->0 があるからね
ありがとうございます。
全然 o(x^n) とは違うんですね。
f(x, y) = x + y - tan(x*y)
(1) (0, 0) の近くで f(x, y) = 0 から y = φ(x) と解けることを示せ。
(2) (1) の φ について、
φ'(0)
φ''(0)
φ^(3)(0)
φ^(4)(0)
φ^(5)(0)
φ^(6)(0)
を Mathematica を使って求めよ。
また、
φ(x) - [-x/(1-x)] = o(x^5) (x → 0)
であることを Mathematica を使って示せ。
↑は、 φ(x) と -x/(1-x) のグラフを -0.99 ≦ x ≦ 0.99 の範囲で描いたものです。
よい精度で近似できていることが分かります。
おしっこを出せと命令されたのにかわりに茶をだしたASKA
与党議員を盗聴している、本省の課長級以上の国家公務員の私邸を盗撮している、と言って
日蓮上人の正法を護持し、日本の真の国主であり平和の王者である池田先生と創価学会をネットで誹謗中傷していた奴ら、
おまえらは、全員ASKAと同じ精神障害ということで警察官に逆らう詐欺師と言う事にこれではっきりとなった。
警察官の３割強はわたくし達創価会員で占められている。
ASKAは必ず再逮捕され、拘置所内で自殺と言う事で間違いなく処理される。
おれたちは絶対負けない。
創価の敵対者達も必ず仏法の力で破滅する事になる。
ASKAの味方をする奴は全員統合失調の薬物中毒者と言う事にもうすでにもうなったんだから、
おとなしく観念しろ。
創価が与党議員を盗聴している、本省の課長級以上の国家公務員の私邸を盗撮している、
という疑惑を流している奴らも同罪ということになる。
池田先生は魂の指導者で日本の国主なのだということは、全世界が一致して認めていることなのだから、
もし仮に公務員や一般国民の生活をご照覧になって指導されたとしてもそれは当たり前のことだ。
このスレにそんな馬鹿は居ないと思うがASKAに肩入れして警官や盗撮の悪口を言ってるやつは
IPアドレスから住居を突き止めて六道地獄を永久に廻る仏罰をくだしてやる。
慰安婦として強制労働を行わせ、心身に消えない傷を負わせた
先日、韓国人元慰安婦のハルモニの方々からお話を聞く機会を頂いた
旧日本軍の兵士達による卑劣で、下賎で、下劣な行為をたくさんされた事を
お聞きし、私は、旧日本軍に対する恥ずかしさ、怒り、悲しさの感情が急に高ぶってしまい
涙と嗚咽は止まらなくなり、思わずその場で土下座をしてしまった
我が国は、国ぐるみで犯罪を行ったのである
私は、1人の日本人として、真の愛国者として、韓国や韓国人ハルモニ達には
一生をかけて未来永劫償うべきであると強く感じた
Raspberry Piで無償で使えるから実は今が一番手軽に使える時代なのに。
(1) x≒0 で
φ(x)＝－x/(1-x)－(1/3)x^6－(4/3)x^7－(10/3)x^8－(20/3)x^9－(59/5)x^10－(99/5)x^11－(497/15)x^12－(284/5)x^13－(31454/315)x^14－…
(2)
φ (0)＝0,
φ '(0)＝-1,
φ "(0)＝-2,
φ^(3) (0)＝-6,
φ^(4) (0)＝-24,
φ^(5) (0)＝-120,
φ^(6) (0)＝-960,
f(x,y) = x+y - tan(x･y)
= x+y - xy - (1/3)(xy)^3 - (2/15)(xy)^5 - …
x≒0 で
φ(x) = - x/(1-x) - (1/3)(x^6)/(1-x)^4 - (2/15)(x^10){1+(17/2)x+(77/2)xx+126x^3+(14129/42)x^4+…}
ウェブの方が情報多いでしょ
基本操作程度なら各大学が公開してるPDFに全部書いてあるし
微積は簡単だろうけど整数問題とかはまだムズイ？
maximaでいくつか入試問題を解いているけど
Mathematicaでも似たような感じなのでは？
maximaってスマホで動くのか
iphoneにはないの？
Mapleの最新版は2016のようですが、2017はいつ出ますか？
アップデートは学生版はないですね。
だけど学生版の価格はすげー安かったような
https://software.univcoop.or.jp/search_result?utf8=%E2%9C%93&unit_page[keyword]=maple&commit=%E6%A4%9C%E7%B4%A2
mapleはとてもいいですよ
5月か６月じゃないですかね
学生版じゃなければアップデートがサポート内に入っているから
無料で2017にアップデートできますよ。学生版は12,960円だから
安いよね。あとhome editionというのもあるけど学生は学生版が
一番やすい
最近、Understanding Maple, Cambridge University Press
と
https://goo.gl/ZpkFCY
が出たよ
というか、このスレはmathematicaスレなので余計なことをこれ以上書きません。
おじゃましました
Maple は普通の手続き型言語ですよね？
Mathematica のほうがプログラミングの勉強にはなるのではないでしょうか？
プログラミング言語として、 Mathematica は実用的なのでしょうか？
関数型言語ってプログラミングしにくくないですか？
括弧の対応がウザイ
見せ方は昔よりましになったが多重になると混乱する
オプションのTrueとかNoneとか間違える
Mathematicaの方言て、どういう意味？
なんかコマンドの名前の多さでいろんなものをカバーしてるようなきがする。
elispみたいなものか
mathematicaを使っていてもプログラミングをしてるという気にならない
IPythonなんてもうMathematicaみたいなもんだし。
きちがい
あと、ご存知ならmaximaにも対応する関数があるか教えて欲しいです
からぬけ出す方法を教えてください。
どうもMATHのバグだといわれるけれど、速度が落ちて困っています。
regarding the debugger tools window.
In 10.0 and 10.1 under Windows once I turn on the debugger
I experience the window persisting even after the menu is unchecked;
it reappears when I close a Notebook window. I suspect this has to do with
the change to individual window menus rather than a menu bar.
Fortunately in my case restarting Mathematica removes the tools window,
however other people have experienced it persisting between sessions; see
Block[{a, c},
{a, _, c}=f[];
{a,c}]
みたいに、関数型言語に良くあるワイルドカード _ みたいなパターンマッチって出来ないのか？
{#[[1]],#[[3]]}&
なんかMathematicaもほんとにオワコン化してきてる印象
￥
http://o.8ch.net/w5kw.png
オワコンと書いてる人がいて驚いた
￥
① 公的年金と生活保護を段階的に廃止して、満18歳以上の日本人に、
ベーシックインカムの導入は必須です。月額約60000円位ならば、廃止すれば
財源的には可能です。ベーシックインカム、でぜひググってみてください。
② 人工子宮は、既に完成しています。独身でも自分の赤ちゃんが欲しい方々へ。
人工子宮、でぜひググってみてください。日本のために、お願い致します。☆☆
Series[Tan[x], {x, 0, 3}]
などと入力すると、出力される
O[x]^4
というような記号の意味を教えてください。
O[x^4] ではなく O[x]^4 と書くのはなぜでしょうか？
http://reference.wolfram.com/language/tutorial/TheRepresentationOfPowerSeries.html
気のせい？
教育関係者とかと私のような個人で欲しいだけの人を差別するのもなんだかなあ
差別なのかなんなのかは知らないけど
今はやや円安で、260円ぐらいになってるけど。
学生版を買えばいいんじゃないか
次のライセンス更新は一般では
機械学習でMathematicaがドラスティックに変わるのはいつ？
ディープラーニングで未解決予想が証明できるようになれるのはいつ？
あと論文の例ぐらいはなんとかなるのではないかと思います。
さきのことはわかりません
Re[x+ iy]
はどうして 単純に x にならないのでしょう？
そこでComplexExpandですよ。
それは関数のリストではないですよね。
ff[x_][i_]:= x^i
FF[y] --->{y,y^2.y^3.y^4,y^5}
あれはDVDかなんかの物理媒体を注文した人の為のものですか？
たとえば他人に譲渡に関して不正をしたら、
訴えることもできるよね。
購入者を知ることはなんにせよ基本だと思うけど。
何か届くとかいうことじゃない、ということでオッケーですか？
今日購入したんですが、マンション名と部屋番号を書き忘れてしまって…
メールで連絡するしかないでしょうか？
そんなもんは別に放っておいてもいいんですか？
どなたかレスください
l = {{1, 2}, {3, 4}}
というリストと
i = {1, 2}
というリストがあったときに、
l[[1, 2]]
の値を取り出すのに、 i を使って取り出す方法はありますか？
g[x_, y_, z_] = x * y * z;
である場合、
argumentCount[f] == 2
argumentCount[g] == 3
となるような関数 argumentCount を作りたいです。
を
Derivative[1, 2, 3][g]
と偏微分するのと同じ結果を
l == {1, 2, 3}
が与えられたときに l を使って計算したいのですが、どうすればいいですか？
ID: mmuQWU7a
基地外臭いｗ
一応できたんですが、
多変数の多項式を次数の昇順に並び替えるにはどうすればいいのでしょうか？
という順序で並べ替えたいです。
多項式の形でそれを構成する単項式の表示順序を変更する方法はありますか？
多変数の Taylor 多項式の計算です。
https://i.imgur.com/ceMoL4L.jpg
sss = Table[{RandomInteger[{0, 2}], RandomInteger[{0, 2}]}, {10}] （＊例＊）
sr = Flatten[Table[{i, # - 1} & /@ Range[(i + j - 1)], {i, 0, 3}], 1]；
Table[Position[sss, sr[[j]]], {j, 1, sr // Length}]//Union
Flatten[%,2]
sr[[%]]
解決策教えていただけないでしょうか...
NMinimize[目的関数、制約条件]、変数]
というとき、たとえば制約条件がλ_1＞0、λ_2＞0、・・・、λ_100＞0のとき、
λ_k＞0,{k,1,100}と表したいのですが、うまくいきません。
制約条件の数が少なければλ_1＞0、λ_2＞0、・・・でよいですが、
1個ずつ入れるには多すぎるのです。
ご教示の程よろしくお願いいたします。
Tableを使って
Table[λ_k, k] >= 0, {k, 1, 100}]}
でござるな
0, (a : Blank[4]) > 0, (a : Blank[5]) > 0, (a : Blank[6]) >
0, (a : Blank[7]) > 0, (a : Blank[8]) > 0, (a : Blank[9]) >
0, (a : Blank[10]) > 0}
λ_k は自然数とする。
このとき解の個数を求めよ
Flatten[Permutations /@ IntegerPartitions[15, {10}], 1] // Length
-> 2002
かな。
素晴らしいですね。
わたしは頭を使うのが苦手で
n100 = 10; A = Array[a, {n100}];
Solve[({(A // Total) == 15, Table[a[ll] > 0, {ll, 1, n100}]}) //
Flatten, A, Integers] // Length
でやりました。
貴君のほうが早く答えが出ますね。
ありがとうございました。
これを止める方法誰かおしえてー
シフトキー周りをエアダスターで掃除
掃除したけどダメでしたぁッ
{g1, g2, x, y, u, v, plr},
g1 = ParametricPlot[s {Cos[phi], Sin[phi]}, {s, r, r + dr}, {psi, phi, phi + dphi}];
u = r {Cos[phi], Sin[phi]};
v = r {Cos[phi + dphi], Sin[phi + dphi]};
x = v - u;
v = (r + dr) {Cos[phi], Sin[phi]};
y = v - u;
plr = Options[g1, PlotRange];
g2 = ParametricPlot[s x + t y, {s, 0, 1}, {t, 0, 1}, plr];
Show[g1, g2]
]
g1とg2を同じPlotRangeで描きたいと思い、↑のコードを書きましたが、
エラーが出てしまいます。
どうすればいいでしょうか？
g2 = ParametricPlot[s x + t y, {s, 0, 1}, {t, 0, 1}, plr];
↓
g2 = ParametricPlot[s x + t y, {s, 0, 1}, {t, 0, 1}];
で動作確認したら
Subscript[x, 1]
を変数として使いたいのですが、
別のところで、
x = 3
などとやるとエラーが出ます。
どうすればいいでしょうか？
例えば
Subscript[x, 1]=123
とすれば
Subscript[x, 1]
は123を返してきます。
しかし、ただのxで参照したら別の変数でしょう？
f[a]=321
とすれば
f[a]
は321を返すんだから説明になってなかったですね。
でも
x_1(ここでは表現が難しい、、）は123を返してきます。
sumList[lis_]:=First[lis//.{x_, y___} -> x + {y}]
sumListの右辺で最初にx+{y}が計算されて{x+y}となるみたいだ。
yはSequenceが入る。
例えば
x+Sequence[a,b,c]は
x+a+b+c
となるようだからlistの要素の総和になる。
を計算すると、
{4, 5}
になりました。
>sumListの右辺で最初にx+{y}が計算されて{x+y}となるみたいだ。
これはなぜなのでしょうか？
.....{x+y+{z,w}}---> x+y+z+w
https://reference.wolfram.com/language/ref/Rule.html
によればRule“->”の右辺は即事に評価すると書いてあるから、
x+{y}はx,yに代入される前に評価されて{x+y}になり
そのあと代入されて、計算が続くということだろう。
いつの間にか定理の自動証明ができるレベルになってたんですね
群の証明なんかが自動的にできてて驚いた
高校1年生の2次関数勉強中。サポートよろ！
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244010380.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244040827.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244240232.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244110703.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244390628.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244540005.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245000222.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245040254.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245180265.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245250755.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245300342.jpg
①井口・千明（東京都葛飾区青戸６－２３－１６）
※盗聴盗撮・嫌がらせつきまとい犯罪者のリーダー的存在／犯罪組織の一員で様々な犯罪行為に手を染めている
低学歴で醜いほどの学歴コンプレックスの塊／超変態で食糞愛好家である／醜悪で不気味な顔つきが特徴的である
②宇野壽倫（東京都葛飾区青戸６－２３－２１ハイツニュー青戸２０２）
※色黒で醜く太っている醜悪黒豚宇野壽倫／低学歴で人間性が醜いだけでなく今後の人生でもう二度と女とセックスをすることができないほど容姿が醜悪である
③色川高志（東京都葛飾区青戸６－２３－２１ハイツニュー青戸１０３）
※色川高志はyoutubeの視聴回数を勝手に短時間に何百何千時には何万回と増やしたり高評価・低評価の数字を一人でいくつも増やしたり減らしたりなどの
youtubeの正常な運営を脅かし信頼性を損なわせるような犯罪的業務妨害行為を行っています
※色川高志は現在、生活保護を不正に受給している犯罪者です／どんどん警察や役所に通報・密告してやってください
【通報先】
◎葛飾区福祉事務所（西生活課）
〒１２４－８５５５
東京都葛飾区立石５－１３－１
℡０３－３６９５－１１１１
④清水（東京都葛飾区青戸６－２３－１９）
※低学歴脱糞老女：清水婆婆 ☆☆低学歴脱糞老女・清水婆婆は高学歴家系を一方的に憎悪している☆☆
清水婆婆はコンプレックスの塊でとにかく底意地が悪い／醜悪な形相で嫌がらせを楽しんでいるまさに悪魔のような老婆である
⑤高添・沼田（東京都葛飾区青戸６－２６－６）
※犯罪首謀者井口・千明の子分／いつも逆らえずに言いなりになっている金魚のフン／親子孫一族そろって低能
⑥高橋（東京都葛飾区青戸６－２３－２３）
※高橋母は夫婦の夜の営み亀甲縛り食い込み緊縛プレイの最中に高橋親父にどさくさに紛れて首を絞められて殺されそうになったことがある
⑦長木義明（東京都葛飾区青戸６－２３－２０） ※日曜日になると風俗店に行っている
関数型で書けばうまくいきますかね？
サンキュー
頑張ってみる
ってのはMathematicaじゃないような気もするけど
自マシンのマルチコアプロセッサー利用なら何もいらない
正規ライセンスあればLightweight Gridも利用できる
https://github.com/WolframResearch/WolframLanguageForJupyter
でnotebook代替もできるぞ！
(deleted an unsolicited ad)
TeXアウトプットもできる
入力はキャラクタベースのTexでやるってこと？
出尽くしってことだろうね
ふうL@Fu_L12345654321
学コン1傑いただきました！
とても嬉しいです！
https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
https://twitter.com/Fu_L12345654321/status/1144528199654633477
(deleted an unsolicited ad)
Threshold[n_] = 10*n;
If[Potential[x] >= threshold[n],
Force[x_, v_] := -k*x-G*v;];
のように、条件にあった時だけIf文内を実行したいのですが、If文の条件判定がうまく機能していないようです。
条件分に変数のある関数が入っていたらまずいのでしょうか。
下のように、条件文の前で関数の値を置き直せば問題ないのでしょうか。
こうすると動くようになったと思いますが、実際のものは煩雑でして正しく動作しているのか分かりません…。
よろしくお願い致します。
Potential[x_] = 0.5*k*x^2;
Threshold[n_] = 10*n;
Potential_value = Potential[x_];
Threshold_value = Threshold[n_];
If[Potential_value >= Threshold_value,
Force[x_, v_] := -k*x-G*v;];
>If[Potential[x] >= threshold[n],
>Force[x_, v_] := -k*x-G*v;];
よくわからないけど、If文の条件が実際に評価されるときに(x,nに）値が入ってるんですか？そうでないと関数Fは未定義のままですよね？
それとThresholdはMathematica自身の関数として定義済みでは？
それと、Potential_value =と_valueをつけると変数にセットできないんじゃありませんか？私の使ってるバージョンではできません。
初心者でトンチンカンな事を聞いてるみたいですみません。
>それとThresholdはMathematica自身の関数として定義済みでは？
これ、そもそも○○[x_]などと定義するのが間違ってるって事でしょうか？
xに値が入るなら○○[x_]と定義する必要がない、と。
>それと、Potential_value =と_valueをつけると変数にセットできないんじゃありませんか？私の使ってるバージョンではできません
たしかにできないみたいです。
すみません。
x,nをループで回して、
Potential = 0.5*k*x^2;
Threshold = 10*n;
If[Potential >= threshold,
Force[v_] := -k*x-G*v;];
というようにすればいいのですかね。
Force[v_] := -k*x-G*v;
と何度も定義し直すのは何か意図があるのですか？
素直にループの外で
Force[v_]:=If[ 0.5 k x^2>=10 n, -k x-G v ,必要なら何か定義...]
と定義しておくのはダメなのですか？
Mathematicaの入門書を読むことをお勧めします。
皆さんどうされてるんでしょう
本読み直します
関数の引数に関数を与えることもできます。
ふうL@Fu_L12345654321
学コン1傑いただきました！
とても嬉しいです！
https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
https://twitter.com/Fu_L12345654321/status/1144528199654633477
(deleted an unsolicited ad)
くだらんものをだすな 資源の浪費だ
↑のθ(～)は何の演算子ですか？
右下見れ
あ、スマホでは出ないんだな。PCで見てくれ
Sum[If[q == p, q, 0], {q, 1, 3}] /. p -> 2
の結果として 2 になること期待したんだけれど，
qが返ってきてしまった…
なして？
だから q だろうね。
Sum[If[q == p, q, 0] /. p -> 2, {q, 1, 3}]
だと２だよ。
Sum は Block を使って変数を置き換えていく。Block は変数を置き換えながら内側の式を評価するが、HoldRest のせいで第2引数の q は評価されずにそのまま残る。
第1引数は評価され置き換わるが、If は第1引数が True または False になるまで開かれないから、結果 >>802 の1行目のように展開され、答えは q となる。
ふうL@Fu_L12345654321
学コン1傑いただきました！
とても嬉しいです！
ｈｔｔｐｓ://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
ｈｔｔｐｓ://twitter.com/Fu_L12345654321/status/1144528199654633477
(deleted an unsolicited ad)
Sum[If[q == p, Evaluate@q, 0], {q, 1, 3}] /. p -> 2
でも2が返るよ。
昔lispやってたので馴染めないんだが、、、。
周囲の長さが一定 2*s である三角形の面積の最大値を求めよ。
ヘロンの公式より、
S = sqrt(s * (s - x) * (s - y) * (s - z))
φ(x, y, z) = x + y + z - 2*s
(grad φ)(x, y, z) = (1, 1, 1) ≠ (0, 0, 0)
0 < x < y + z = 2*s - x
0 < y < z + x = 2*s - y
0 < z < x + y = 2*s - z
だから、
0 < x < s
0 < y < s
0 < z < s
でなければならない。
φ(x, y, z) = 0
0 < x < s
0 < y < s
0 < z < s
という条件下で、
f(x, y, z) = (s - x) * (s - y) * (s - z)
を最大化する (x, y, z) を求めて、面積 S を計算すればよい。
↓は制約条件を満たす点の集合、いわゆる実行可能領域です。
https://imgur.com/XzhCaP7.jpg
s = 1 としました。
↑実行可能領域の正三角形の重心で面積が最大になることが分かりますね。
n m
Σ Σ x[i,j]
i=1 j=1
j≠i
みたいなやつ。意味わかります？
ただ、それはわかってますが、もうちょっとスマートに書きたいということです。
その方が早そうだし
発散するときは
x[i_, j_] := 0 /; i == j;
x[i_, j_] := X[i, j];
Sum[x[i, j], {i, 1, 3}, {j, 1, 3}]
Out=X[1, 2] + X[1, 3] + X[2, 1] + X[2, 3] + X[3, 1] + X[3, 2]
とか？
Mathematicaの本を見ても基本的なことしか書いていないですよね。
f[x_, y_] := x^2 + y^3
g[t_] := {t^2, 3*t+1}
とします。
f[g[t]] とはできないですよね？
f[ｖ_] := v[[1]]^2 + v[[2]]^3
g[t_] := {t^2, 3*t+1}
としないといけないんですか？
1個しか与えてないのだから未定義扱いになるのは当然では？
f@@g[t] で如何
正規行列はユニタリ行列によって対角化できるという定理が書いてあります。
Wolfram Language 12で、 A = U*D*U^* となる U, D を計算する関数を作りました。
orthogonalizationU[A_] := Module[
{eigenvs, U, DIA},
eigenvs = Eigenvalues[A];
DIA = DiagonalMatrix[eigenvs];
eigenvs = DeleteDuplicates[Eigenvalues[A]];
U = Transpose @ ((Flatten[#, 1])& @ (Orthogonalize /@ (NullSpace /@ ((A - #*IdentityMatrix[Length[A]])& /@ eigenvs))));
{U, DIA}
]
以下の定理が書いてありますが、 A = U*T （T が上三角行列）の場合しか証明されていません。
2.5.17 【定理】
任意の n 次正則行列 A はユニタリ行列 U と、対角成分が正の実数であるような上三角行列（下三角行列でもいい） T の積 U*T （T*U でもいい）として
一意的に表わされる。
(2)
A^T = U * T （T は上三角行列）
と分解し、両辺の転置を取ると、
A = T^T * U^T （T^T は下三角行列）
(3)
A の列たちをリバースした行列を A' とする。
A' = U * T （T は上三角行列）
と分解する。
T の列たちをリバースした行列を T' とする。
A = U * T' (T' は左上三角行列）
である。
U の列をリバースした行列を U' とする。
T' の行をリバースした行列を T'' とする。
A = U' * T'' （T'' は下三角行列）
である。
(4)
(3)の手順にしたがって、
A^T = U * T （T は下三角行列）
と分解する。
A = T^T * U^T （T^T は上三角行列）
(2)～(4)のケースの一意性についても
A = U*T （T が上三角行列）のタイプの分解の一意性から自明ですよね。
4つのタイプの分解を行う関数を作りました。
gs[U_, a_] := Append[U, Simplify[normalize[a - Sum[(a.u)*u, {u, U}]]]];
gramschmidt[A_] := Module[{U = {}}, Do[U = gs[U, A[[i]]], {i, 1, Length[A]}]; U];
qrdecom[A_] := Module[{U, T},
U = gramschmidt[A];
T = Table[If[j > k, 0, A[[k]].U[[j]]], {j, 1, Length[A]}, {k, 1, Length[A]}];
{Transpose[U], T}];
qrdecom2[A_] := Module[
{B, U, T},
B = Reverse[A];
U = gramschmidt[B];
T = Table[If[j > k, 0, B[[k]].U[[j]]], {j, 1, Length[B]}, {k, 1, Length[B]}];
{Transpose[Reverse[U]], Reverse[Transpose[Reverse[Transpose[T]]]]}];
qrdecom3[A_] := Module[{B, U, T},
B = Transpose[A];
U = gramschmidt[B];
T = Table[If[j > k, 0, B[[k]].U[[j]]], {j, 1, Length[B]}, {k, 1, Length[B]}];
{Transpose[T], U}];
qrdecom4[A_] := Module[
{B, C, U, T},
B = Transpose[A];
C = Reverse[B];
U = gramschmidt[C];
T = Table[If[j > k, 0, C[[k]].U[[j]]], {j, 1, Length[C]}, {k, 1, Length[C]}];
{Transpose[Reverse[Transpose[Reverse[Transpose[T]]]]], Reverse[U]}];
こいつは数学の本スレで連投を繰り返している荒らしです
レスしないでください
並列処理のリモートカーネルとして使えないものか
ここに入っていればRaspberry Piも要らないね
Mathematica 6のフロントエンドで動けば助かったけど、6.0 FE+12.1 Engineの組み合わせはちょっと無理っぽかった
Macで動いたって言うブログ見たけれど
並列計算用サブカーネルにのみ使って、正規の8カーネル＋Engine 8カーネルで16カーネル並列計算できないか
なんとなくやり方はわかったけど、これは後でやってみる
やりすぎたか
無償のラズパイMathematicaフロントエンドか別に使えるバージョン新しめのMathematicaフロントエンドを持っているなら
でもMathematicaのフロントエンドを使えば高額な有償Mathematicaと同じになっちゃう夢を見た
Wolfram EngineのMathKernelはMathematicaパッケージのMathKernelと全く同じで、特に機能制限されてない
スマホでMathematicaとか未来だ
許諾しているのなら、Jyupiter NotebookをMathematica Notebookに置き換えても構わないよな？
それを許諾しないとは書いてない
ネットワーク経由のリモート計算でならロハでもできるというだけ
それが無理ならおとなしくローカルでJyupiter Notebook
ただJyupiter NotebookとMathematica Notebookの使用感が全く異なっており、ファイル互換性もない
なんとかMathematica NotebookをフロントエンドにしてNotebookを統一したいというのが夢のきっかけだった
何か変わったか？
https://reference.wolfram.com/language/howto/ConnectToARemoteKernel.html.ja
3Dグラフィックが遅いのも入力支援機能がないのもフロントエンドの問題だから解決しない
同じマシンで8 x Workerカーネル動かして計算すること
マスターカーネルを一つにして、Ｎ個のマシンにWorkerカーネルを分散して、8N Workerカーネルで並列計算してパフォーマンス
を見てみたいのだが、うまくいってない
N x 8 Wokerカーネルの並列計算は出来ているが
1 x 8N workerカーネルの並列計算が出来てない
今のところSurface Pro7も使ってN=6
Xeon E5-2643v3 (実6コア/仮想12スレッド) x2
Core i7-8700T（実6コア/仮想12スレッド）x2
Core i7-4790K（実4コア/仮想8スレッド）x1
Core i5-1035G4（実4コア/仮想8スレッド）x1
それぞれのコアで8 Workerカーネル動かして並列計算してる
てことは初代iMacよりも性能上なのか
学生はStudent版買うよりこっちのほうがいいな
バグなのかな？
もっとコアが欲しい
１台のPCで8カーネルの制約を超えた並列計算が可能
並列カーネルの状態表示が壮観
難しい
引数がどうにもうまく渡せない
学術の巨大掲示板群 - アルファ・ラボ ttp://x0000.net
数学 物理学 化学 生物学 天文学 地理地学
IT 電子 工学 言語学 国語 方言 など
PS 連続と離散を統一した！
ttp://x0000.net/topic.aspx?id=3709-0
微分幾何学入門
ttp://x0000.net/topic.aspx?id=3694-0
これでlightweight gridに頼らなくても済みそうだ
・ローカルカーネル x8
・ローカルにあるWolfram engineのカーネル x8
・２つの別マシンのカーネル x8x2
を起動して32カーネル
Windows10オプションのOpenSSH（クラ・サバ）を使用
（フロントエンドや数学的な）サポートはしないけど
ってとこじゃない？。
完全タダは無理じゃない？
自分は正規ライセンスのMathematicaにWolfram Engineを組み合わせてできる範囲の並列分散処理をやってる
ラズパイ買えば、含まれる無償MathematicaとWolfram Engineの組み合わせでできることは拡大するが、ラズパイ
のコストはかかる
また、この方法はバックエンドの処理が重い場合には有効だけれども、複雑なグラフの描画とかフロントエンドの処理
が重い場合にはあまり意味がない
Wolfram EngineをPCで
ってパターンならラズパイ分の値段だけでいける？
いける
それにメリットを感じるかどうかは、やりたいことによる
で、（スレーブカーネルは並列処理のときのみ）
フロントエンド をラズパイ上のMathematicaで、 マスターカーネル ― スレーブカーネル をPC上のWolfram Engineで
分散処理することができる
また、フロントエンド ― マスターカーネル をラズパイ上のMathematicaで、スレーブカーネルをPC上のWolfram Engineで
分散処理させることもできる
並列処理計算なんてしない人が多いんだろうけど
その場合はマスターカーネルをPCに任せるという話になる。マスターカーネルでの処理が重ければ分散処理の意義が出てくる
だろう
どうやってやるの？
これ、ParallelTableとかでシームレスに32並列行けるの？
フロントエンド をラズパイ上のMathematicaで、 マスターカーネル―スレーブカーネル をPC上のWolfram Engineで分散処理は、
https://community.wolfram.com/groups/-/m/t/94917
https://support.wolfram.com/ja/12495?src=mathematica
http://www3.econ.fukuoka-u.ac.jp/usage/doku.php/windows
を基に設定した。これらを理解したうえで設定すればできた。
フロントエンド―マスターカーネル をラズパイ上のMathematicaで、スレーブカーネルをPC上のWolfram Engineで分散処理は、
Windows10のオプション機能であるOpenSSHを使えば実現できる
行ける
ありがとうございます。できるかどうか分かりませんが、とりあえず読んでみようと思います。
コネクトする方法が書かれているが、自分は前者で実現した。
PC側でカーネルをリスンモードで起動するバッチファイルを作成し、先にカーネルを起動する。フロントエンドのカーネル環境設定
でフロントエンドからPCに接続するようにMLOpenの引数を設定したうえで、マスターカーネルを起動する。
LinkNameは自分で適当にポートを２つ決め打ちして両方で合わせればいい。
やろうと思えば、PC上のカーネルをOpenSSHなどでリモート起動することもできるはず。
残念
解説お願いします
seedList =
With[{vertices = .1 N[PolyhedronData["Icosahedron"][[1, 1]]]},
Join[Map[{#, 2} &, vertices],
Map[{# + {1, 1, 1}, -2} &, vertices]]];