〓 Mathematica 伍 〓
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/11/11(金) 09:02:28.82, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,）)
! | ! .,'-.{ ! !|; |`､.}ﾞ!.! |. ! ヽ.
', ', |Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l ! | ﾞ',
ヽ､', l:!Kノ}. f:_.)iﾞi: ﾘ ! l ル
| l!iヽﾞ- ' , .!__:ﾉ ﾞ ,ﾘ l ﾘ'´
. ',|!!､ r‐┐ ｀ ノ'. /,ｲ ／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣
'i!ﾞ､ヽ､ ﾞｰ' _, ィ,:',:''´ ＜ Mathematica に関する話題はここに書いてね！
ﾞ:､ｨ､jヾｰ::: 'iﾍ .ノ',ﾘ. ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
,､- '´ ヽ､ﾞ､ { ｀＞"､
/＼＼ ', } //｀ヽ
過去スレ
〓Mathematica〓
http://science2.2ch.net/test/read.cgi/math/985023298/
〓 Mathematica 2 〓
http://science3.2ch.net/test/read.cgi/math/1078534285/
〓 Mathematica 3 〓
http://science6.2ch.net/test/read.cgi/math/1121413040/
〓 Mathematica 四 〓
http://kamome.2ch.net/test/read.cgi/math/1197270001/ (鯖飛びでログ消滅)
http://mimizun.com/log/2ch/math/1197270001/
〓 Mathematica 5 〓
http://kamome.2ch.net/test/read.cgi/math/1285859504/
0461１３２人目の素数さん
2013/04/08(月) 15:06:52.38単に線を沢山書いて作っているだけだよ
0462あぼーん
NGNG0463１３２人目の素数さん
2013/04/12(金) 19:26:18.21もし,受験前にこれ知ってるとやる気なくすよな。
正直テクを使わず入力するだけで答えが得られる
2007^2007 % 17
を求めよ
なんて、そのまま入力したら 1 が一瞬で得られる
高木貞次の一変数の微積をしっかり計算できることなんて、
未だに文科省の理系数学の目標なんだな。あほらし
大学入試でPCとMathematica与えて問題解かせるようにすると数学教育も一変するんだがな。
実際そっちの方が研究/開発では実戦的
0464１３２人目の素数さん
2013/04/12(金) 20:33:33.29数学教育は何も変わらない
0465１３２人目の素数さん
2013/04/14(日) 16:29:58.34すくなくとも大学の課題とかはmalabあたりを駆使しないと解けない問題がでる
変わってないのは日本だけ。
チャンコロが作った科挙制度にならって延々同じような試験をいまでもやらかしてる。
資格試験とかもめーんな同類だけど。意味の無い記憶に頼る試験はナンセンス過ぎるんだな
0466１３２人目の素数さん
2013/04/14(日) 19:31:20.23で、アメリカでは
>大学入試でPCとMathematica与えて問題解かせるように
なってるの？
0467１３２人目の素数さん
2013/04/15(月) 18:46:47.43大学入試問題は相当変わるな
0468１３２人目の素数さん
2013/04/17(水) 19:12:45.55とおのむかしからそーなってる
0469１３２人目の素数さん
2013/04/17(水) 20:32:48.860470古賀潤一郎
2013/04/17(水) 23:22:51.070471有能な政治家を貶め、無能な政治家を持ち上げて国力を削ぐカルト教団
2013/04/21(日) 18:41:36.64・沢山の人が偏った意見を一貫して支持する
偏った意見でも、集団の中でその意見が信じられていれば、自分の考え方は間違っているのか、等と思わせる手法
・不利な質問をさせなくしたり、不利な質問には答えない
誰にも質問や反論をさせないことにより、誰もが皆、疑いなど無いんだと信じ込ませる手法
↑マスコミや、カルト教団のネット工作員がやっていること
TVなどが、偏った思想や考え方に染まった人間をよく使ったり、左翼を装った人間にキチガイなフリをさせるのは、視聴者に、自分と違う考え方をする人間が世の中には大勢いるんだなと思わせる効果がある。
.
0472１３２人目の素数さん
2013/04/29(月) 19:29:17.05世界中、日本以外は２６５ドルで買えるんじゃないのか!!!?
0473１３２人目の素数さん
2013/04/29(月) 20:45:13.90日本代理店 J, H, S
0474１３２人目の素数さん
2013/04/30(火) 11:12:18.92どの国に変えても２６５ドルだよ。日本以外は。
しかしどんだけ消費者舐めてんだ？
ふざけるにもほどがあるだろ！
>>473
どこのこと？
0475１３２人目の素数さん
2013/04/30(火) 12:32:26.02いろいろできることを知った。
もうこれでいいや。
0476１３２人目の素数さん
2013/05/01(水) 00:57:32.06あとIPython NotebookはかなりMathematica Notebookに近いよな。
0477１３２人目の素数さん
2013/05/01(水) 20:37:17.270478１３２人目の素数さん
2013/05/02(木) 10:23:23.46日本は中国より信用無いのか？
0479１３２人目の素数さん
2013/05/02(木) 12:26:44.490480１３２人目の素数さん
2013/05/02(木) 13:25:09.430481１３２人目の素数さん
2013/05/02(木) 20:05:02.12高くても買う奴がいるかぎり
値下げはしない
0482１３２人目の素数さん
2013/05/02(木) 20:38:41.37とwolframに聞いてみたんだが、なかなか返事が来ないね。
なんか、wolframはまともな会社に思えなくなってきたな、、。
万一英語版買っても日本語版との差額をふんだくるぞとかいうし。
まるで…
0483あぼーん
NGNG0484１３２人目の素数さん
2013/05/21(火) 13:03:24.93だったら英語版は?と聞くと、今までそうだったから同額設定（？）なんだそうだ。
とにかく日本では高く売るということらしい。
日本語化するだけでhome editionという意図とはかけ離れた価格になるし、
バージョンアップは遅くなるし、ろくな事が無い。
日本語化なんてすべきじゃなかったね。反って買えないよ。
俺もあきらめたよ。無料のcdfplayerとかsageを使う事にする。
sageいれるとmaximaやR、、、とかもいろいろ入るみたいだから面白そうだし。
0485１３２人目の素数さん
2013/05/25(土) 14:47:02.21接続しているようにできれば、US版を購入&アクティベート
可能だよね？
VPNを利用した接続元(IPアドレス)を変更できるフリーのサービス
があるみたいです。
https://www.okayfreedom.com/
カード情報を入力できるほど信用できるかどうか疑問だけど。
誰か試した人いませんか？
0486１３２人目の素数さん
2013/05/29(水) 18:30:38.60難しいと思います。
0487１３２人目の素数さん
2013/06/06(木) 05:51:05.86スタンドアローンの環境で使っている研究所とかあるし。
0488１３２人目の素数さん
2013/06/07(金) 17:30:00.70ただ、その時の情報は記録されてて、起動時にチェックをしてる
アクティベート承認した環境じゃないだろ？と判断したら、
エラーで再要求される
オンラインアクティベートに接続元の地域は関係ないと思う
ネットに繋げないなら、パスワードをクリエイトすればいいだけの話しだし
0489１３２人目の素数さん
2013/06/26(水) 15:33:22.86同じような人いない？
0490１３２人目の素数さん
2013/06/27(木) 11:38:08.21バージョン7だっけ？ それくらいからそんなもんだと思うけど。
なんだか最近なら↓を試してみたら？
http://support.wolfram.co.jp/kb/3274
0491１３２人目の素数さん
2013/06/28(金) 00:52:30.29ver5.2，7，8と使ってきたけど9だけ異常に重い…
リンク先のはまだ試してないからやってみるよ
0492489
2013/06/28(金) 14:27:29.72というか，インストールした直後から重かったので解決できないのは当然だった．
一応，[評価]欄にある「動的更新の実行」のチェックを外すと従来どおりの速度になった．
毎回チェックはずさないといけないけど．
前のバージョンではチェック入ったままだったんだけどなぁ
0493489
2013/06/28(金) 14:52:21.24ネットへの接続を不許可にしたらヘルプが激早になった．
ネットにいちいちアクセスしていたせいだったのか．
0494１３２人目の素数さん
2013/06/28(金) 19:18:06.87プロクシの設定したら治りそうだな
ネット接続しないなら、そのままの方が簡単でオッケーだが
0496１３２人目の素数さん
2013/07/14(日) NY:AN:NY.ANFor文の中の変数iを、回す文章の中のファイル名に入れたいのですが、どう書けばいいのでしょうか？
例えば、Import="フォルダ名//基本名i"という形で、iによって読み込むファイルを変えたり、基本名i=○○○という変数をFor文の中で定義したりしたい場合はどうすればいいのでしょうか。
0497１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANできるならそのとおりにやればいい
0498１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANID:cF5MhH6Y!0499１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANかな。
0500496
2013/07/15(月) NY:AN:NY.AN>>499
すみません、左辺も右辺も、どう書けばいいのかわかりません・・・。
フォルダ内にtime1.datとtime2.datとtime3.datがあり、それぞれをリストの形で読み込みたい、という質問です。
つまり
t1=Import["フォルダ名\\time1.dat,"Table""];
t2=Import["フォルダ名\\time2.dat,"Table""];
t3=Import["フォルダ名\\time3.dat,"Table""];
という作業を（実際はもっと多くの連番ファイルがあるので）Ｆｏｒ文で簡潔におこないたいのです。
教えてくださっているのは、右辺の部分でしょうか？
0501１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANt[i]=Import["フォルダ名\\time"<>i<>".dat","Table"]
]
ていう感じか？動くかどうかは確認していないが。
0502１３２人目の素数さん
2013/07/15(月) NY:AN:NY.AN○For[i=1,i<=3,i++
0503496
2013/07/15(月) NY:AN:NY.ANありがとうございます！
>>502を踏まえて、
For[i=1,i<=3,i++,
t[i]=Import["フォルダ名\\time"<>i<>".dat","Table"]
]
とやってみたり、i++のあとの,を;にしてみたりしたのですが、
StringJoin::string: String expected at position 2 in フォルダ名\time<>1<>.dat. >>
Import::chtype: First argument フォルダ名\time<>1<>.dat is not a valid file, directory, or URL specification. >>
というエラーが出てしまいました。
また、右辺に関してもt[1]やt[2]・・・というように配列に入れていくのではなく、t1やt2という名前の変数として定義したいのですがどうしたらよいのでしょうか。（t0095のように、tのあとに0を付けることもできるようにしたいためです）
0504１３２人目の素数さん
2013/07/16(火) NY:AN:NY.AN0505１３２人目の素数さん
2013/07/16(火) NY:AN:NY.AN0506１３２人目の素数さん
2013/07/16(火) NY:AN:NY.ANまず、499の書いている通り、
<>i<>
でなく
<>ToString[i]<>
でないとだめだと思う。
次に、たとえば、t00095とかにしたいなら
NumberForm[95, 5, NumberPadding -> {"0", ""}]が使えると思う。
よくわからないがTostringしないと"t"とつながらなかった。
最後に
ToExpression["a345=123"]とかつかえば、
a345に123が代入される。
0507１３２人目の素数さん
2013/07/16(火) NY:AN:NY.AN>Tostringしないと"t"とつながらなかった。
Headをみればわかりますよ。
NumberForm[95, 5, NumberPadding -> {"0", ""}]//Head
->NumberForm(*これでは<>:StringJoinは機能しない*)
ToString@NumberForm[95, 5, NumberPadding -> {"0", ""}]//Head
->String(*これなら<>:StringJoinは機能できる*)
0508１３２人目の素数さん
2013/07/16(火) NY:AN:NY.ANドキュメントセンターのグラフのこと？
事前にデフォルトオプションを変更して、
SetOptions[Plot, Exclusions -> None];
Plot[SawtoothWave[x], {x, -2, 2}]
とするか、個別にオプションを指定して
、
Plot[SawtoothWave[x], {x, -2, 2}, Exclusions -> None]
としないとつながらないよ。
0509505
2013/07/16(火) NY:AN:NY.AN0510１３２人目の素数さん
2013/07/17(水) NY:AN:NY.AN当時のやっつけの記述ですので、洗練はされてません。
(*filePath = "DataFolderにPassが通っていないときに記述"*);
fileName =
Table[
StringJoin["AgtABCD_image", StringTake["0" <> ToString[n], -2], ".gif"],
{n, 0, 63}(*データ数:64*)
];
imageList = Map[Import[#(*, Path -> filePath*)] &, fileName];
(*fileNameにリテラル値をいれる場合の例*)
literalName =
Table[
StringJoin["literal_image", StringTake["0" <> ToString[n], -2], ".gif"],
{n, 0, 63}
];
fileName = literalName;
※StringTake["0" <> ToString[n], -2]部は、今だとNumberFormの方がスマート？かも。
0511１３２人目の素数さん
2013/08/17(土) NY:AN:NY.AN0512１３２人目の素数さん
2013/08/25(日) NY:AN:NY.AN0513１３２人目の素数さん
2013/08/25(日) NY:AN:NY.ANWolframもアドビみたいに
ライセンスで金とるだけか
0514１３２人目の素数さん
2013/09/15(日) 19:19:21.980515計算
2013/09/25(水) 22:26:18.76約４時間程度と思います。
自分のPCの速度を見極めたいです。
http://www.geocities.jp/dirac_equation/index3.htm
0516１３２人目の素数さん
2013/09/26(木) 20:54:18.390517計算
2013/09/26(木) 21:03:18.141000秒程度のはずです。
http://www.geocities.jp/dirac_equation/index4.htm
0518計算
2013/09/26(木) 21:26:45.13下記計算を、昨日の２３時ころからやっておりますが、未だに終わりません。
PDF参照
本来なら、約６０秒×２５６／（６０秒×６０分）＝４時間 程度で終わるはずです。
なぜでしょうか？
http://www.geocities.jp/dirac_equation/index3.htm
の計算結果です。
http://www.geocities.jp/dirac_equation/time.pdf
0519１３２人目の素数さん
2013/09/26(木) 21:37:05.960520計算
2013/09/26(木) 22:02:30.42メモリー不足ですと、下記のように62.343750 秒と表示されても、実際はもっと計算時間がかかっているのでしょうか？
x= 0,y= 0
62.343750, Null <
0521１３２人目の素数さん
2013/09/27(金) 15:50:14.81ガンマ行列の性質を使って整理してから g[4〜15]=0 に注意すれば Mathematica に頼るまでもない計算のようだが。
0522計算
2013/09/27(金) 23:31:26.45そうでしょうか？このガンマ行列は、通常の４×４ではありません。256×256が１６個あります。またg[4〜15]=0は、ここで、特別に示しただけです。０でないケースを計算して、どうなるか？を検討するのが目的です。
ところで、
やっぱり、おかしいです。
下記計算を、一昨日の２３時ころからやっておりますが、未だに終わりません。
本来なら、約６０秒×２５６／（６０秒×６０分）＝４時間 程度で終わるはずです。
なぜでしょうか？
http://www.geocities.jp/dirac_equation/index3.htm
0523１３２人目の素数さん
2013/09/27(金) 23:59:17.47D次元（今はD=16）の場合
・γ^μγ^ν+γ^νγ^μ
・γ^μγ_μ
・tr(γ^μγ^ν)
と、それから
・γ^νγ^μγ_ν
・tr(γ^μγ^νγ^ργ^σ)
がそれぞれどうなるか御存知だろうか？それらを使えば、件のトレースはずっと簡単に表される。
>ガンマ行列の性質を使って整理
とはそういうこと。
0524計算
2013/09/28(土) 00:16:20.99手計算はメンドクサイし、間違いも発生します。
すでに、コードは完成しておりますので、 Mathematica に頼ります。
このコードを完成させるのに、むちゃ苦労しました。
0525１３２人目の素数さん
2013/09/28(土) 00:26:35.49「完成している」コードの動作が「おかしい」とはこれ如何に
0526計算
2013/09/28(土) 00:27:21.12256×256の１６個のガンマ行列を使って、下記を得るのが目標です。
256×256の４個のガンマ行列なら、簡単です。
http://www.geocities.jp/dirac_equation/fig.pdf
0527計算
2013/09/28(土) 00:31:26.98そこを質問しているのです。でも多分メモリー不足だと思います。
１つの結果を出すのに、１週間も時間がかかるのに、８GBのメモリーですから、、
0528１３２人目の素数さん
2013/09/28(土) 00:46:13.960529１３２人目の素数さん
2013/09/28(土) 02:05:15.77「コードは完成してる(はずだ)」けど、「正しく実行されてない」と疑ってるってことでしょ。
原因はマシン性能の可能性があるから、もっとパワーのあるマシンで実行してみて欲しいと。
0530137
2013/09/28(土) 09:17:35.90質問１
y = {2, 5, 1, 3, 4}
のような １〜５の整数の５個の乱数のリストを作るには、どのようなコマンドを使用すれば良いでしょうか？
質問２
z = x + x^2 + 4 r + t^2
の場合、式の項数は、４個ですが、どのようなコマンドを使用すれば４個と表示するでしょうか？
0531137
2013/09/28(土) 10:03:09.770532137
2013/09/28(土) 10:47:59.04繰り返し計算をｎ回させ、Xに１から１０の乱数を発生させて、その乱数が２以下の場合に、ｘ＝ｙとして、その値ｙを表示させて、繰り返し計算をそこで、ストップさせることは可能でしょうか？
下記は、失敗例です。どこを修正すれば良いでしょうか？
For[n = 1, n <= 10, n++,
v = RandomSample[Range[10]];
x = v[[n]];
Print["n=", n, ",x=", x];
If[x < 1, y = x & Print[y]];
];
0533137
2013/09/28(土) 13:24:27.210534計算
2013/09/28(土) 13:32:05.57もっと速く計算する方法はないでしょうか？
例えば、マッキントッシュ用のmathematicとパソコンを使用すれば、速く計算可能でしょうか？その場合、どのくらい速くなるのでしょうか？
その他の方法はないでしょうか？
（※；但し、メモリーは８GBで不足してます。購入中です。）
http://www.geocities.jp/dirac_equation/index6.htm
0535１３２人目の素数さん
2013/09/28(土) 20:23:26.69q0とかも、256×256の行列でしょうか？
やっていることは、256×256の行列をいくつか用意して、
それらの間で、何らかの単純な計算をしているということでいいのでしょうか？
逆行列を求めるとかでない限り、そんなに時間がかかるプログラムにはならないと思います。
Cとか、自分でメモリの管理ができる言語で作り直すほうが、手っ取り早いのでは？
0536計算
2013/09/28(土) 20:48:51.29例えば、C言語で、２回微分させるコードを書くのに、５０行必要としました。
Mathematicaなら、１行で済みます。
この計算は、一般相対論の計算より、かなり複雑なので、C言語で書き直すのは、難しすぎて、出来ないです。
メモリの管理なのですが、マックパソコンでは、Mathematicaのメモリー管理が出来るのではないでしょうか？また、windowsよりも、速いのではないでしょうか？
0537１３２人目の素数さん
2013/09/28(土) 21:25:51.91速さはWindowとかMacとかによるものではない
0538計算
2013/09/29(日) 09:30:05.17Windows 8 - 128GB
Windows 8 Professional - 512GB
Windows 8 Enterprise - 512GB
まで、増設できるので、windowsPCのメモリー増設で対応できそうな気がします。
http://bto-pc.jp/btopc-com/repair/ram-max-beginner-pc.html
0539１３２人目の素数さん
2013/09/29(日) 17:31:05.37まず>>523を考えた方がいい。
0540計算
2013/09/29(日) 18:27:47.85当然、523の計算はやっております。ただし、それはここでアップしたコードの計算後です。
このやり方は、教科書の記載通りです。この方法で、４×４ の４個のガンマ行列を使い、以下の正しい答えを導いてます。
http://www.geocities.jp/dirac_equation/fig.pdf
現在、１６GBに増設して、計算中です。
これで、メモリー不足になって計算できない場合は、以下のPCの購入も検討するつもりです。
http://www.dospara.co.jp/5shopping/detail_prime/2/56/3208/0/
0541１３２人目の素数さん
2013/09/29(日) 20:52:48.60>当然、523の計算はやっております。
それを知ってて
http://www.geocities.jp/dirac_equation/index6.htm
のxとyの和をいちいち成分で足すとかどういうマゾ趣味か、とは思うなぁ
0542１３２人目の素数さん
2013/09/29(日) 21:41:14.200543１３２人目の素数さん
2013/09/30(月) 09:10:34.45例えば Σ[k=1,n]k=n(n+1)/2 を知ってたら 1+2+3+4+…+99999+100000 を一つずつ足し算するようなことはせんだろ。計算氏がやってるのはそういうことだ。
0544計算
2013/10/01(火) 21:32:41.36そうかもしれませんが、この方法で、
http://www.geocities.jp/dirac_equation/fig.pdf
の正しい結果を導出しております。
但し、ご指摘通り、遠回りをしている可能性はあります。しかし、この方法は教科書通りです。
ところで、計算ができない理由が解りました。
２５６回の繰り返し計算で、項数が多い式を
y1 = y1 + s1;
で、作って、後で
Print[FullSimplify[ExpandAll[y1]]];
で、一気に整理しようとしたから、メモリーがパンクして計算できませんでした。
１回の繰り返し計算ごとに
Y1=FullSimplify[ExpandAll[y1+s1]];
で、整理するように、すれば、計算できるはずです。（只今、計算中）
0545１３２人目の素数さん
2013/10/01(火) 22:44:55.40もしかして、
>D次元（今はD=16）の場合
>・γ^μγ^ν+γ^νγ^μ
>・γ^μγ_μ
>・tr(γ^μγ^ν)
>と、それから
>・γ^νγ^μγ_ν
>・tr(γ^μγ^νγ^ργ^σ)
>がそれぞれどうなるか
の結果を知らないのか？知ってたら書いてみて。
0546計算
2013/10/01(火) 22:56:02.51y1 = 1024*(-64*k2^2*m^2 - 64*k3^2*m^2 + 64*m^4 - 48*m^2*p0^2 + 48*m^2*p1^2 - 16*k2*m^2*p2 + 48*m^2*p2^2 - 16*k3*m^2*p3 + 48*m^2*p3^2 + 49*k2^2*p0*q0 + 49*k3^2*p0*q0 - 63*m^2*p0*q0 + 49*p0^3*q0 - 49*p0*p1^2*q0 -
49*p0*p2^2*q0 - 49*p0*p3^2*q0 - 49*k2^2*p1*q1 - 49*k3^2*p1*q1 + 63*m^2*p1*q1 - 49*p0^2*p1*q1 + 49*p1^3*q1 + 49*p1*p2^2*q1 + 49*p1*p3^2*q1 + 112*k2*m^2*q2 - 98*k2*p0^2*q2 + 98*k2*p1^2*q2 + 49*k2^2*p2*q2 -
49*k3^2*p2*q2 + 63*m^2*p2*q2 - 49*p0^2*p2*q2 + 49*p1^2*p2*q2 + 98*k2*p2^2*q2 + 49*p2^3*q2 + 98*k2*k3*p3*q2 + 98*k2*p3^2*q2 + 49*p2*p3^2*q2 +
7*(7*k3^2*p3 + p3*(-7*k2^2 + 9*m^2 + 7*(-p0^2 + p1^2 + p2^2 + p3^2)) + 2*k3*(8*m^2 + 7*(-p0^2 + p1^2 + p2*(k2 + p2) + p3^2)))*q3 + k1^2*(-64*m^2 + 49*(p0*q0 + p1*q1 - p2*q2 - p3*q3)) +
k0^2*(64*m^2 + 49*(p0*q0 + p1*q1 + p2*q2 + p3*q3)) + 2*k1*(-8*m^2*(p1 - 7*q1) + 49*((-p0^2)*q1 + p1^2*q1 + (p2*(k2 + p2) + p3*(k3 + p3))*q1 + p1*(k2*q2 + k3*q3))) -
2*k0*(-8*m^2*(p0 - 7*q0) + 49*((-p0^2)*q0 + (p1*(k1 + p1) + p2*(k2 + p2) + p3*(k3 + p3))*q0 + p0*(k1*q1 + k2*q2 + k3*q3))));
纏めると
y1=-24064m^4+26112m^2s+25088m^2u-25088su;
になります。
0547１３２人目の素数さん
2013/10/10(木) 06:53:26.430548１３２人目の素数さん
2013/10/10(木) 20:08:33.09, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/
. | ＼ ∠イ ,ｲイ| ,`-' |
| l＾,人| ` `-' ゝ | このスレは馬と鹿と豚さんばかりね。
| ` -'＼ ｰ' 人
| /(l ＿＿／ ヽ、
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0549１３２人目の素数さん
2013/10/11(金) 10:22:32.09例えば A, B, C を互いに非可換としてその他の数や変数は可換扱いにしたいです．
このとき
Simplify[10*A*20*B-10*B*20*A]
のような計算があったとして，
結果が 200 (A**B-B**A) になるようにしたいのですが…
0550549
2013/10/11(金) 13:02:22.60http://math.ucsd.edu/~ncalg/
0551１３２人目の素数さん
2013/11/08(金) 09:12:06.47例えば｛a,b,c,d,e}（a〜eは数値）といったリストとA（数値）があって、
a,b,c,d,eによる線形結合で、ただし係数は整数、Aに一番近い値を求めたいです。
このときA＝１a+２b+３c+４d+５e
みたいにどういう形なのかも具体的に知りたいです
こういうこ機能のある関数はありますか？
0552１３２人目の素数さん
2013/11/09(土) 01:02:05.66http://reference.wolfram.com/mathematica/ref/LatticeReduce.html
の例題のアプリケーションを見ると良いと思う。
0553551
2013/11/09(土) 13:11:47.000554１３２人目の素数さん
2013/11/10(日) 01:48:52.81>>552のリンク先のLatticeReduceのドキュメントだけど、
アプリケーションの2番目の例題で
自明な初期関係：
a={
{1, 0, 0, 0, -a0},
{0, 1, 0, 0, -a1},
{0, 0, 1, 0, -a2},
{0, 0, 0, 1, -a3},
{0, 0, 0, 0, -a4}
};
ってなってるけど、間違ってるよね？
3つの例を統一的に表現できるコード書いてみた。
t = 3^(1/3);
v = {1, t, t^2, t^3, t^4};
n = 7;
c = Round[10^n v];
A = Join[IdentityMatrix[Length[v]], -Transpose[{c}], 2]
B = LatticeReduce[A]
A.Append[c, 1]
B.Append[c, 1]
B[[;; , ;; -2]] // MatrixForm
r = B[[;; , ;; -2]].v
N[r]
0555１３２人目の素数さん
2013/11/10(日) 01:59:57.67x*π+y*e=123.456
これを上記のコードにあてはめて
v = {Pi, E, 123.456};
n = 10;
c = Round[10^n v];
A = Join[IdentityMatrix[Length[v]], -Transpose[{c}], 2]
B = LatticeReduce[A]
A.Append[c, 1]
B.Append[c, 1]
B[[;; , ;; -2]] // MatrixForm
r = B[[;; , ;; -2]].v
N[r]
を実行してみると
-6747π+3256e+100*123.456≒0
が得られて、変形すると
(6747/100)π-(3256/100)e≒123.456
となるんだけど、これは求めたい解じゃないよね？
つまり、x,yは整数じゃないといけないんだよね？
となるとLatticeReduceで解けなくない？
0556１３２人目の素数さん
2013/11/10(日) 15:06:57.820557551
2013/11/11(月) 01:07:35.82ご指摘の通り係数は整数のものがいいので、LatticeReduceは使うのはやめました
この問題に関しては効率は悪いですがとりあえず別のプログラムを組んで解決しました
でもこういう処理するいい方法がないのか気になりますね
0558551
2013/11/11(月) 01:10:45.510559551
2013/11/11(月) 01:15:59.880560１３２人目の素数さん
2013/11/11(月) 19:19:11.00ttp://twicsy.com/i/AehpGb
ハート形の立体なんておもしろいなぁとか思いつつ、
記載されてる数式を眺めてたんだけど、
「320」ってなんだっ！？と思ったわけ。
明らかに必要ないからね。
んで、いろいろ調べてるうちに以下に辿り着いた。
ttp://demonstrations.wolfram.com/EquationsForValentines/
Taubinっていうのがオリジナルでそれを定数倍したのが
Trott（ウルフラム社員）になってることがわかる。
つまり、本来は同じ結果が出なきゃいけないのに定数倍することによって
じゃっかん滑らかになったからってそれを大々的にアピールするとかｗｗｗ
ウルフラムって馬鹿なの？自社製品を宣伝するんじゃなくて、
自社製品のバグを宣伝してどうすんだよｗｗｗ
■ このスレッドは過去ログ倉庫に格納されています