1 名前:132人目の素数さん mailto:sage [04/03/06 09:51] _....._{{ 〃 , - ' ,..、、.ヾ{{フ'⌒`ヽ、 / ,:', -‐‐` ´ '´⌒ヽ ヾ:、 . ,' ,'´ ,ィ ,ィ ,' , `ヽ', ',-< ,' .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/
835 名前:132人目の素数さん mailto:sage [05/02/28 22:28:13 ] Mathematica link for Excel使っている人いる?
836 名前:132人目の素数さん mailto:sage [05/03/02 21:39:17 ] すまへん。LOG[1+a x]にフェルミ分布関数の微分をかけたような 被積分関数を0から無限大まで数値積分したい。 そのとき、Methodはどれを選ぶと早いのでしょうか?実際は もっと複雑な被積分関数なので、まあまあのせいどで 現実的な時間に終わってほしい。台形?やっぱり? リカーションとか分割とかナイスな設定あったら教えて。
837 名前:132人目の素数さん mailto:sage [05/03/03 05:56:09 ] >>835 持ってるけど使ってないぞ。
838 名前:132人目の素数さん [05/03/03 11:00:48 ] age
839 名前:132人目の素数さん [05/03/03 15:32:07 ] やっぱり計算が遅い。対数関数と指数関数の積をインテグランドに したものはやっぱり計算が遅い。前者の対数関数を xの多項式で近似し、後者をガウス関数で近似すればええのんかな? 今のところ、積分区間をゼロから無限大まで引っ張っているけど どうせ重要な積分区間は限られているし、その狭い区間だけ 計算させるようにすれば格段に短縮されるのかな?もう本当 つらい。
840 名前:132人目の素数さん [05/03/03 16:59:32 ] Do[式, {i, min_i, max_i, step_i}] → step_i おきに i を min_i から max_i まで変えながら式を繰り返す. てなこっちゃけど、結果を対数表示させたいので、min_i*step_i^n ってループさせる方法は?
841 名前:132人目の素数さん mailto:sage [05/03/03 19:46:19 ] >>840 式をそう書いたら言いやん。
842 名前:132人目の素数さん [05/03/11 15:57:38 ] >>376 全く下らん問題やな。数値積分にして、メソッドのところを モンテカルロにしたら?次にインテグレートで出力結果をえる。 そして二つを比較。
843 名前:132人目の素数さん mailto:sage [05/03/12 01:59:26 ] >>842 230日前のレスにレスする君に脱帽。
844 名前:132人目の素数さん mailto:sage [05/03/12 15:33:27 ] 5.1.1のうpでーたが来たよう
845 名前:132人目の素数さん [05/03/13 07:05:41 ] Mathematicaってどこの数学科でも授業で使い方教えてくれるんですか?
846 名前:132人目の素数さん [05/03/18 22:57:03 ] 円周率の平方根は有理数だそうです。 In[1] := Limit[x - Sqrt[Pi], x -> 19151686109235253480650512155425821931577/ 10805181810206786582137323831443648757793] Out[1] = 0
847 名前:132人目の素数さん mailto:sage [05/03/18 23:26:28 ] >>846 Block[ {$MaxExtraPrecision=100}, Limit[x-Sqrt[Pi], x-> 19151686109235253480650512155425821931577/ 10805181810206786582137323831443648757793] ]
848 名前:132人目の素数さん mailto:sage [05/03/19 03:48:57 ] >>846 4.2だとそうはならんな。
849 名前:132人目の素数さん mailto:sage [05/03/19 17:04:04 ] 昨年卒業したので、Student版からPro版へのアップグレード 期限がもうすぐなんだけど、卒業証明書って必須なの?実家に置いてきちゃったよ・・・。
850 名前:132人目の素数さん mailto:sage [05/03/19 23:37:44 ] >>849 バーシティーウェーブに問い合わせしたら?
851 名前:132人目の素数さん [2005/03/24(木) 17:06:51 ] 保守上げ
852 名前:132人目の素数さん mailto:sage [2005/03/24(木) 22:16:56 ] >>849 後輩か教授に金渡して買ってもらえ。
853 名前:132人目の素数さん [2005/03/28(月) 07:56:51 ] interpolation function 良い感じ。 しかし NIntegrateInterpolatingFunction 意味不明 どうやったら積分収束がはやくなるの
854 名前:132人目の素数さん [2005/03/29(火) 03:50:46 ] MATLAB とかMAPLEをサイトライセンスで入れていると、 MATHEMATICAのサイトライセンスの料金を高くするというような 公正な競争を阻害するライセンスの販売システムになっている。 大学ではみんなボイコットして、MAPLEに乗り換えよう。
855 名前:132人目の素数さん mailto:sage [2005/03/30(水) 22:12:03 ] Mathematicaが直感的に操作できて(゚д゚)ウマー
856 名前:132人目の素数さん mailto:sage [2005/03/31(木) 01:37:06 ] Mathematica 高すぎ
857 名前:132人目の素数さん mailto:sage [2005/03/31(木) 14:09:32 ] Export[ , ,"GIF"] が動きません。 こんなエラーが出ます。 Export::fmterr: 無効なPNM形式です "GIF" 以外はJPEGでもBMPでも動くようです。 documents.wolfram.com/v4-ja/RefGuide/Export.html 一体なにが問題なのか分るかたお教えください。
858 名前:つづき mailto:sage [2005/03/31(木) 14:14:15 ] documents.wolfram.com/mathematica/functions/Export ここを見ても、正しいことをしているようです。 バージョン5.1です
859 名前:132人目の素数さん [皇紀2665/04/01(金) 14:51:23 ] >>818 昔作った関数でこういうのがあった。 TF = TableForm[PaddedForm[#, {6, 5}] & /@ # & /@ #, Alignments -> Right, TableSpacing -> {0, 0}] &; Table[…]//TF みたいな使い方をしてる。 もう忘れたから、 Mathematicaに貼り付けて 「選択項目のヘルプ」でしらべてちょ。
860 名前:132人目の素数さん [2005/04/02(土) 21:13:23 ] 二分木なデータ構造をMathematicaでやろうとしたらどう実装すればいいですか? Lisp風なリスト?
861 名前:132人目の素数さん mailto:sage [2005/04/02(土) 23:30:44 ] Mathematica自体がLisp風な訳だが。
862 名前:132人目の素数さん [2005/04/03(日) 11:01:02 ] >>861 でも、リストの深いところを破壊しようとすると全部書き換える事になって遅くなりませんか?
863 名前:862 mailto:sage [2005/04/03(日) 11:24:11 ] つうかリストの破壊が基本的にできない気がしてきました。書き換えしかない。
864 名前:132人目の素数さん mailto:sage [2005/04/03(日) 12:19:28 ] >>862 二分木のデータを変更すると、そこに至るための上の部分の木も再生成されるが、log(n)だ。 だから、普通Schemeで二分木使う時はset!でデータ書き換えたりしないと思うのだが。 Mathematicaでは「書き換え」はシンボルに対してしかできないだろう。 一部を書き換えると全部が「コピー」または「再生成」される。 どうしても書き換えたいなら、シンボルによる二分木を作っておいて、あとでシンボルに値を アサインするようにすればいいんじゃないか。
865 名前:862 mailto:sage [2005/04/03(日) 22:19:37 ] >>864 そうですよね、「コピー」または「再生成」しかできない。 (ちょっとこっちな言語の用語に詳しくないもんで、すいません。) Unique[]ってのがあるみたいなんで、それでノードを乱打して繋げてみます。 ありがとうございます。 #今日リストの大コピー合戦で実装してみたら、 #Length[Flatten[リスト]] = 10000あたりでもう超ウルトラ遅い。 #ノード数は200ぐらいなのにorz
866 名前:132人目の素数さん mailto:sage [2005/04/06(水) 23:34:24 ] 教えてください Simplify[Abs[a + b*I], {a ∈ Reals, b ∈ Reals}] とすると, Abs[a + I b] という答えが返ってきて困っています.aもbも実数と言ってるんだから,もっと簡単になるだろ! Sqrt[a^2 + b^2] という解を得るには,どうしたらよいのでしょうか?
867 名前:132人目の素数さん mailto:sage [2005/04/07(木) 00:35:10 ] >>866 機械にとっては、次のどちらが簡単だと考える方が自然かね? Abs[Plus[a, Times[Complex[0, 1], b]]] Power[Plus[Power[a, 2], Power[b, 2]], Rational[1, 2]]
868 名前:866 mailto:sage [2005/04/07(木) 01:36:28 ] 自己解決しました.867さん,ご意見ありがとう ComplexExpand[Abs[a + b*I]] でやってくれました
869 名前:132人目の素数さん [2005/04/08(金) 13:02:41 ] >>868 > 自己解決 自己解凍みたいw
870 名前:博士 ◆OBpuRHyGa. mailto:もっちーAGE [2005/04/08(金) 14:58:47 ] lim f(x)=sin∫f(x)となるようなf(x)をヒルベルト空間を用いて x→∞ 求めよ。バナハ空間でもよい。
871 名前:132人目の素数さん [2005/04/08(金) 15:17:58 ] >>870 積分の範囲は?
872 名前:博士 ◆OBpuRHyGa. mailto:もっちーAGE [2005/04/08(金) 16:12:06 ] 0→∞
873 名前:132人目の素数さん mailto:sage [2005/04/08(金) 20:41:10 ] >>870 なんでヒルベルト空間なんているんだよ。 lim[n→∞]f(x)=0 ∫[0,∞]f(x)dx=0 にすればいいだけじゃん。 たとえば f[x]=(-1)^[x]/[1+x/2]
874 名前:873 mailto:sage [2005/04/08(金) 20:45:48 ] >>873 lim[n→∞]f(x)=0 ↓ lim[x→∞]f(x)=0
875 名前:132人目の素数さん [2005/04/14(木) 15:41:02 ] >>868 MyRe[H_] := 1/2*(ComplexExpand[H, TargetFunctions -> {Re, Im}] + ComplexExpand[Conjugate[H], TargetFunctions -> {Re, Im}])
876 名前:132人目の素数さん [2005/04/24(日) 10:50:44 ] 保守のため、上昇
877 名前:132人目の素数さん [2005/04/28(木) 00:33:04 ] Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals]で Conjugate[a + I b]となり, a-I bにならないのですがどうしてでしょう? ちなみに、 Simplify[Conjugate[I b], a ∈ Reals && b ∈ Reals]は -i b になるし、 Simplify[Conjugate[a+I ], a ∈ Reals && b ∈ Reals]は -i+a になります。
878 名前:132人目の素数さん mailto:sage [2005/04/28(木) 19:14:20 ] >>877 Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals]でa - I bになるが。 @Mathematica 5.0
879 名前:132人目の素数さん [2005/04/28(木) 19:55:03 ] >>878 Mathematica4.0.1 mac版 では Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals]は Conjugate[a + I b] です。878さんのはwindows版ですか?
880 名前:878 mailto:sage [2005/04/29(金) 19:29:24 ] >>879 Win版だけど、バージョンの違いだと思われ。OS変えるだけで わざわざカーネルのソース書き換えないだろう。
881 名前:132人目の素数さん mailto:sage [2005/05/01(日) 10:16:29 ] Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals, ComplexityFunction -> LeafCount] V4系にはConjugate[..]のほうが単純に見えるんだろう。
882 名前:132人目の素数さん [2005/05/03(火) 10:00:50 ] age
883 名前:132人目の素数さん [2005/05/10(火) 20:25:40 ] hos
884 名前:132人目の素数さん [2005/05/11(水) 19:13:01 ] ちょっと質問なのでが Mathematicaを使ってアイパターンを描画させたいのですがMathematicaでできるでしょうか? gnuplot使っていたのですがアイパターンを描けずMathematicaに移行したいと思っています
885 名前:132人目の素数さん mailto:sage [2005/05/11(水) 21:38:09 ] gnuplotでできなかったのは、誰にも聞かなかったからで、 Mathematicaでできたとしたら、それは誰かに教えてもらえたからでしょうね。
886 名前:132人目の素数さん [2005/05/12(木) 17:56:57 ] プログラム板でレスなかったのでこちらで質問させてください。 CUIで使ってるときに、1行に表示される文字数を変更する方法はありませんか? 例えば、 In[1]:= N[Pi,100] Out[1]= 3.1415926535897932384626433832795028841971693993751058209749445923078\ > 16406286208998628034825342117068 としたときに、結果を改行せずに表示してほしいのですが。
887 名前:132人目の素数さん mailto:sage [2005/05/12(木) 18:51:26 ] Mathematica 5.01のMac版とSolaris版で試したのですが、 In[22]:= 1. 10^(-309) +0.0 Out[22]= 0. 上のようになるのはバグでなく機械精度の問題ですか? また、この結果を 1. 10^(-309) となるようにする方法ってありますか?
888 名前:887 mailto:sage [2005/05/12(木) 19:18:22 ] すいません In[22]:= 1. 10^(-309) +0.0 Out[22]= 0. こうなるのは、有効桁数のこと考えたら当たり前ですね。逝ってきまつ
889 名前:132人目の素数さん mailto:sage [2005/05/12(木) 21:03:13 ] >>886 無理
890 名前:132人目の素数さん [2005/05/15(日) 10:27:27 ] さっきランダムウォークの コードを書いた。はろうか
891 名前:132人目の素数さん [2005/05/15(日) 10:57:05 ] 今起きたばかりで コンドームと一瞬読めてしまった俺ガイル
892 名前:132人目の素数さん mailto:sage [2005/05/15(日) 16:55:26 ] 迷ってるのなら、貼れ
893 名前:132人目の素数さん [2005/05/19(木) 22:01:32 ] >>892 今日また、コードを書いた。
894 名前:132人目の素数さん [2005/05/20(金) 18:44:08 ] Ver.5.1 9ラック版ある?
895 名前:132人目の素数さん mailto:sage [2005/05/20(金) 20:06:58 ] 〃〃 ガッ ___〃 ∧_∧ | | (・∀・ ,,) .=| |==○○= ヽ 从|__|∩ ヽ 、(⌒)っ ( `Д´)/ ヽ_) / ←>>894
896 名前:132人目の素数さん mailto:sage [2005/05/21(土) 21:57:03 ] >>894 刑務所行って探して来い
897 名前:132人目の素数さん [2005/05/25(水) 18:10:48 ] おまいらMathematicaMark 2004の結果教えてくり やり方は <<Utilities`Benchmark` BenchmarkReport[] おいらはチンコパッド T42 PenM 1.7GHz 512MB Mathematica5.1.1で1.59
898 名前:132人目の素数さん mailto:sage [2005/05/31(火) 12:33:44 ] 3点O(0,0,0)、A(1,0,0)、B(0,1,1)を通る平面上に、 Aを中心とする半径1の円をy軸の周りに1回転してできる回転体を Mathematicaを使って書きたいんですけど どんなコマンドにしたらいいですか?
899 名前:ChaosicSoul ◆/yaJbLAHGw [2005/05/31(火) 12:50:54 ] Re:>>898 パラメータ表示で曲面を描くコマンドがあったはずだ。一個のパラメータで円を作って、もう一個のパラメータは回転用に使えばいいだろう。 円のパラメータ表示は、(1+cos(t),√(2)sin(t)/2,√(2)sin(t)/2)でいいだろう。 これに回転用の行列を掛ければ完成だ。パラメータ表示を縦ベクトルにして、 ((cos(s),0,-sin(s)),(0,1,0),(sin(s),0,cos(s)))^T を左から掛ければいいだろう。
900 名前:ChaosicSoul ◆/yaJbLAHGw [2005/05/31(火) 12:51:19 ] ちなみに正しいスペルは chaotic のようで。
901 名前:898 mailto:sage [2005/05/31(火) 12:55:36 ] コピペでMathematicaに貼り付けたらできちゃうように コマンドを作ってもらってもいいでしょうか…
902 名前:ChaosicSoul ◆/yaJbLAHGw [2005/05/31(火) 13:02:29 ] Re:>>901 今Mathematicaが手元に無いから何とかして自分で作ってくれ。
903 名前:898 mailto:sage [2005/05/31(火) 13:07:28 ] >902 ありがとうございました。 自分はあんまりMathematicaは使いこなせないんで、 ChaosicSoulさんがMathematicaを使えるようになったり、 他の方でも、できるようならコマンドを書いていただけるとうれしいです。 皆さんお願いします。
904 名前:132人目の素数さん mailto:sage [2005/05/31(火) 15:26:04 ] そこまで人に頼るのかよ
905 名前:Zeno ◆5nZQbNmQPs [2005/05/31(火) 17:56:11 ] Reply-to:>>903 君は将来性がない人間だな。 そんな簡単なグラフは自分で書くようにしろ。 ________________________________________________________ セオドア・カジンスキー 「産業社会とその未来」 (aka. マニフェスト) の抜粋 www.cnn.com/US/9509/unabomber/09-19/am/excerpt.html
906 名前:132人目の素数さん mailto:sage [2005/05/31(火) 18:25:23 ] >>903 何故かこっち(Methematica 5.0)では途中で固まるが、これであってるはず ParametricPlot3D[Transpose[{{Cos[s], 0, -Sin[s]}, {0, 1, 0}, {Sin[s], 0, Cos[s]}}].{1 + Cos[t], Sin[t]/Sqrt[2], Sin[t]/Sqrt[2]}, {t, 0, 2π}, {s, 0, 2π}];
907 名前:898 mailto:sage [2005/05/31(火) 21:37:46 ] >>905 英語ばっかりだ〜 読めん… >>906 ありがとー! やってみます。
908 名前:132人目の素数さん mailto:sage [2005/06/01(水) 00:12:28 ] 学校の宿題を代わりにやってくれるスレはここですか
909 名前:ねえねえ [2005/06/01(水) 09:57:48 ] science3.2ch.net/test/read.cgi/math/1051275230/l50 >>297 誰か私にちからを…
910 名前:ChaosicSoul ◆/yaJbLAHGw [2005/06/01(水) 10:37:09 ] Re:>>909 Cのファイル入出力でバイナリモードを使えばできる。何が改行コードかに注意。
911 名前:ねえねえ [2005/06/01(水) 11:14:48 ] 承知!!
912 名前:ねえねえ [2005/06/01(水) 11:27:30 ] >>910 ありがとー 出来そうです。
913 名前:132人目の素数さん [2005/06/02(木) 13:40:42 ] u, v を実数, k を定数としたとき、 -5u^{4}v + 10u^{2}v^{3}-v^{5} = k のグラフをUV平面に書きたいのですが, 方法が分かりません。 とりあえず k=1 の時を書こうとしましたが、 Plot を使うには v について解く必要がありそうですし, Solve とか試してみましたがエラーで全然動きません。 どなたか教えて下さい。
914 名前:913 mailto:sage [2005/06/02(木) 14:31:55 ] Plot[{{Root[1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 1]}, {Root[ 1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 2]}, {Root[ 1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 3]}, {Root[ 1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 4]}, {Root[ 1 + 5*u^4*#1 - 10*u^2*#1^3 + #1^5 &, 5]}}, {u, -10, 10}] とかしてみたらそれらしいのが出ました(k=1)。 何か違う気もしますが,もう少しいじってみます。
915 名前:132人目の素数さん mailto:sage [2005/06/02(木) 15:06:21 ] >>913 ImplicitPlotで。 UVP[k_] := ImplicitPlot[ - 5 u^4 v + 10 u^2 v^3 - v^5 == k, {u, -10, 10}, {v, -10, 10}] UVP[1]
916 名前:915 mailto:sage [2005/06/02(木) 15:07:04 ] <<Graphics`ImplicitPlot` ね。
917 名前:132人目の素数さん mailto:sage [2005/06/02(木) 22:11:37 ] >>915 うまいね、UVPという関数を作る感じか。
918 名前:918 [2005/06/02(木) 23:00:40 ] 9=1+8 9-1=8
919 名前:919 mailto:sage [2005/06/02(木) 23:02:02 ] クイック ワイパー
920 名前:132人目の素数さん [2005/06/07(火) 11:13:05 ] 今日、大学生協行ったら5.2へのフリーアップグレードつきのパッケージが出てた。 俺が2週間前に買ったときは無かったのに ┌─┐ |も.| |う | │来│ │ね│ │え .| │よ .| │ !!.│ ショボーン バカ └─┤ プンプン ヽ(´・ω・`) ノヽ(`Д´)ノ (`Д´)ノ ( `Д) | ̄ ̄ ̄|─| ̄ ̄ ̄|─| ̄ ̄ ̄|─□( ヽ┐U 〜 〜  ̄◎ ̄ . ̄◎ ̄  ̄◎ ̄ ◎−>┘◎
921 名前:132人目の素数さん [2005/06/07(火) 16:51:27 ] ここは大学で出た宿題がわかんないんだけど聞いてもいいのかな?
922 名前:132人目の素数さん mailto:sage [2005/06/07(火) 17:21:59 ] 5.2ってまだ出てないよね いつ出るんだろ
923 名前:132人目の素数さん mailto:sage [2005/06/07(火) 18:05:57 ] >>921 いいんジャマイカ?
924 名前:921 [2005/06/07(火) 23:07:59 ] やったね^^ありがと 定点A(a,0)と楕円(x^2/2)+y^2=1上の動点Pとの距離APの最小値bを求めよ また、aを横軸、bを縦軸にとって、aとbの関係をグラフに示せ っていう問題なんですが、条件を置いてからのそこから先の展開ができません。
925 名前:132人目の素数さん mailto:sage [2005/06/08(水) 02:15:02 ] スレ違いだから簡単に Pを(√2 cosθ,sinθ)と置いて距離APの2乗を導き、微分=0を調べればいいんでは つか高校レベルじゃない?arccosとか出てくるから大学なのか
926 名前:132人目の素数さん mailto:sage [2005/06/09(木) 18:57:38 ] 5.1.1がリリースされたみたい。 5.1自体つい半年前に出たばかり。 ただいま、某サイトからダウソ中。
927 名前:132人目の素数さん mailto:sage [2005/06/09(木) 20:22:38 ] 5.1.1は前からあるよ。 つかワレ厨は氏ね!
928 名前:132人目の素数さん mailto:sage [2005/06/09(木) 20:47:05 ] ウルフラムのホムペからも金払えばダウソ出来るよ。
929 名前:132人目の素数さん mailto:sage [2005/06/10(金) 02:35:20 ] 既知
930 名前:132人目の素数さん [2005/06/15(水) 18:34:10 ] m
931 名前:132人目の素数さん [2005/06/17(金) 02:36:27 ] LiveGraphics3Dで3Dグラフィクスを作ってるんだが、Mathematicaでちゃんと PlotRange使ってプロット範囲指定してるのにWriteLiveFormで変換してHTMLで見ると PlotRangeで指定したことがまるっきり無視されるんだが、どうすればいいんだ?
932 名前:132人目の素数さん mailto:sage [2005/06/17(金) 08:12:56 ] jisaku.u-satellite.net/uploader/files/869.gif
933 名前:132人目の素数さん [2005/06/17(金) 12:26:16 ] インテルオワットル
934 名前:132人目の素数さん mailto:sage [2005/06/17(金) 20:17:49 ] ソンナコトワカットル
935 名前:132人目の素数さん [2005/06/18(土) 12:28:34 ] matlabにしてもmathematicaにしても、ここまで差が出るとは。 理系の研究室にopteronが売れまくるのも当然だな