トップページ⇒math
1001コメント317KB

〓 Mathematica 2 〓

レス数が900を超えています。1000を超えると表示できなくなるよ。
0001１３２人目の素数さんNGNG
　　　　　　　　_....._{{ 〃
　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-<
　 　 ,' 　　.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/
0809１３２人目の素数さん05/02/17 15:52:01
>>807
Countを使うなら:
Count[Table[Random[Real, {1, 10}], {20}], n_?(3 <= # <= 8 &)]

Length@Select[...]も悪くはない。
0810１３２人目の素数さん05/02/17 15:53:16
nは要らなかった。_?(...)
0811１３２人目の素数さん05/02/17 16:19:54
>>808
risutoはねぇだろ
081280705/02/17 17:51:15
>>808-811
ありがとうございました。809さんので無事にできました。
0813１３２人目の素数さん05/02/17 19:52:30
>>811

榊原さんの教科書風にしたわ
0814１３２人目の素数さん05/02/18 12:56:39
軽蔑すべき輩だと思う
0815１３２人目の素数さん05/02/18 19:15:37
わざとだろ(w
突っ込まれるの待ってたんじゃね?
0816１３２人目の素数さん05/02/19 16:26:06
質問なんですが
ttp://documents.wolfram.com/v5/Built-inFunctions/NumericalComputation/EquationSolving/FurtherExamples/NBMLImages/NDSolve/NDSolve_79.ja.gif
この画像の3行目の頭にある
u(1,0)[1,t]
はInputFormではどう書くのですか？
0817１３２人目の素数さん05/02/19 16:47:28
Derivative[1, 0][u][1, t]
でした。スマソorz..
0818１３２人目の素数さん05/02/19 20:18:09
Mathematicaからまともな表を表示する方法はありますか？
ちなみにLinkForExcelは可変長リストを返す関数に関して糞なのでダメ。
TeX系はインタラクティブじゃないのでイヤです。
0819１３２人目の素数さん05/02/19 21:42:06
Do[Plot…],{i,0,100}]
って感じでパラメーター変えながら100枚ほどプロットしたのですが、
それぞれのグラフのタイトルとして i=1 とかって表示させるにはどうしたらいいですか？
0820１３２人目の素数さん05/02/19 22:46:20
>>818
これじゃだめかい？だめだろうねぇ。
t = Table[Subscript[a, i, j], {i, 5}, {j, 5}];
DisplayForm[FrameBox[GridBox[t, RowLines -> True, ColumnLines -> True]]]

>>819
Table[Plot[Sin[i x], {x, 0, 2π}, PlotLabel -> "i=" <> ToString[i]], {i, 3}]
0821１３２人目の素数さん05/02/19 23:43:40
>>818
Javaで自作。

Link for Excelが標準的な解なわけだから、
> 可変長リストを返す関数に関して糞
とやらをなんとかして使うのがいいと思うが。
082281905/02/19 23:43:50
>>820
多謝
0823１３２人目の素数さん05/02/20 08:58:22
こんど、PC新しくするんだけど、Mathematicaには
Pentium4とAthlon64のどちらが向いていますか？
今のところ考えているのはPen4 3.4GHz(EM64T)か
Athlon64 3500+くらいのやつです。
0824１３２人目の素数さん05/02/20 11:08:39
Linuxなら64bitのMathematicaがあるよ。Windowsは知らん。
082581805/02/20 18:08:45
>>820-821
ありがとうございます。
とりあえずはLinkForExcel+Excel上のグチャグチャで凌ごうと思います。
0826１３２人目の素数さん05/02/21 10:41:09
>>824
なるほど64bit版ですか。今のところ32bit Windows版を考えてるんですが、安定してるなら
64bit Linux版に乗り換えたいところですね。

どなたかPen4とAthlon64の両方で使っているという方がいましたら、感想などよろしくお願いします。
0827１３２人目の素数さん05/02/22 12:56:28
<< Statistics'Continuous Distributions'

で入れると、どうやっても

Get::noopen: Statistics'Continuousを開くことができません．

でエラー返されちゃうんですが。
バージョンは4です。
0828１３２人目の素数さん05/02/22 13:00:44
>>827
<<Statistics`ContinuousDistributions`

じゃないの？　ContinuousとDistributionsの間に空白が入るのはおかしい。
082982705/02/22 15:01:59
>>828
解決しました。
「'」と「`」が間違っていました。
ありがとうございました。
0830１３２人目の素数さん05/02/22 16:25:36
>>829
バックダッシュは間違えやすいよね
0831１３２人目の素数さん05/02/22 19:15:36
>>830
バッククウォートと言います。

ダッシュってのは横線。'をダッシュというのは日本の中高校数学だけ。
0832１３２人目の素数さん05/02/23 02:32:07
ちなみに'の正式名称は"prime"ね。
0833１３２人目の素数さん05/02/27 19:33:25
上の方と同じような質問なのですが

<< Graphics`ImplicitPlot`
ImplicitPlot[x y - y^2 == 3, {x, -2, 2}]

と打つと
Set::"write": \
"\!\(Options[ImplicitPlot]\)のタグ\!\(ImplicitPlot\)\
はプロテクトされています．"

とでるのですが何がいけなかったでしょうか･･
0834１３２人目の素数さん05/02/27 22:21:26
>>833
俺んとこではでないぞ@Mathematica5
0835１３２人目の素数さん05/02/28 22:28:13
Mathematica link for Excel使っている人いる？
0836１３２人目の素数さん05/03/02 21:39:17
すまへん。LOG[1+a x]にフェルミ分布関数の微分をかけたような
被積分関数を０から無限大まで数値積分したい。
そのとき、Methodはどれを選ぶと早いのでしょうか？実際は
もっと複雑な被積分関数なので、まあまあのせいどで
現実的な時間に終わってほしい。台形？やっぱり？
リカーションとか分割とかナイスな設定あったら教えて。
0837１３２人目の素数さん05/03/03 05:56:09
>>835
持ってるけど使ってないぞ。
0838１３２人目の素数さん05/03/03 11:00:48
age
0839１３２人目の素数さん05/03/03 15:32:07
やっぱり計算が遅い。対数関数と指数関数の積をインテグランドに
したものはやっぱり計算が遅い。前者の対数関数を
ｘの多項式で近似し、後者をガウス関数で近似すればええのんかな？
今のところ、積分区間をゼロから無限大まで引っ張っているけど
どうせ重要な積分区間は限られているし、その狭い区間だけ
計算させるようにすれば格段に短縮されるのかな？もう本当
つらい。
0840１３２人目の素数さん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
ってループさせる方法は？
0841１３２人目の素数さん05/03/03 19:46:19
>>840
式をそう書いたら言いやん。
0842１３２人目の素数さん05/03/11 15:57:38
＞＞376

全く下らん問題やな。数値積分にして、メソッドのところを
モンテカルロにしたら？次にインテグレートで出力結果をえる。
そして二つを比較。
0843１３２人目の素数さん05/03/12 01:59:26
>>842
230日前のレスにレスする君に脱帽。
0844１３２人目の素数さん05/03/12 15:33:27
5.1.1のうpでーたが来たよう
0845１３２人目の素数さん05/03/13 07:05:41
Mathematicaってどこの数学科でも授業で使い方教えてくれるんですか？
0846１３２人目の素数さん05/03/18 22:57:03
円周率の平方根は有理数だそうです。

In[1] := Limit[x - Sqrt[Pi],
x -> 19151686109235253480650512155425821931577/
10805181810206786582137323831443648757793]
Out[1] = 0
0847１３２人目の素数さん05/03/18 23:26:28
>>846
Block[
{$MaxExtraPrecision=100},
Limit[x-Sqrt[Pi], x-> 19151686109235253480650512155425821931577/
10805181810206786582137323831443648757793]
]
0848１３２人目の素数さん05/03/19 03:48:57
>>846
4.2だとそうはならんな。
0849１３２人目の素数さん05/03/19 17:04:04
昨年卒業したので、Ｓｔｕｄｅｎｔ版からＰｒｏ版へのアップグレード
期限がもうすぐなんだけど、卒業証明書って必須なの？実家に置いてきちゃったよ・・・。
0850１３２人目の素数さん05/03/19 23:37:44
>>849
バーシティーウェーブに問い合わせしたら？
0851１３２人目の素数さん2005/03/24(木) 17:06:51
保守上げ
0852１３２人目の素数さん2005/03/24(木) 22:16:56
>>849
後輩か教授に金渡して買ってもらえ。
0853１３２人目の素数さん2005/03/28(月) 07:56:51
interpolation function
良い感じ。
しかし
NIntegrateInterpolatingFunction　意味不明
どうやったら積分収束がはやくなるの
0854１３２人目の素数さん2005/03/29(火) 03:50:46
MATLAB とかMAPLEをサイトライセンスで入れていると、
MATHEMATICAのサイトライセンスの料金を高くするというような
公正な競争を阻害するライセンスの販売システムになっている。
大学ではみんなボイコットして、MAPLEに乗り換えよう。
0855１３２人目の素数さん2005/03/30(水) 22:12:03
Mathematicaが直感的に操作できて（ﾟдﾟ）ｳﾏｰ　
0856１３２人目の素数さん2005/03/31(木) 01:37:06
Mathematica 高すぎ
0857１３２人目の素数さん2005/03/31(木) 14:09:32
Export[ , ,"GIF"]
が動きません。
こんなエラーが出ます。

Export::fmterr: 無効なPNM形式です

"GIF" 以外はJPEGでもBMPでも動くようです。

http://documents.wolfram.com/v4-ja/RefGuide/Export.html

一体なにが問題なのか分るかたお教えください。
0858つづき2005/03/31(木) 14:14:15
http://documents.wolfram.com/mathematica/functions/Export

ここを見ても、正しいことをしているようです。
バージョン5.1です
0859１３２人目の素数さん皇紀2665/04/01(金) 14:51:23
>>818

昔作った関数でこういうのがあった。


TF = TableForm[PaddedForm[#, {6, 5}] & /@ # & /@ #,
Alignments -> Right, TableSpacing -> {0, 0}] &;


Table[…]//TF

みたいな使い方をしてる。

もう忘れたから、
Mathematicaに貼り付けて
「選択項目のヘルプ」でしらべてちょ。
0860１３２人目の素数さん2005/04/02(土) 21:13:23
二分木なデータ構造をMathematicaでやろうとしたらどう実装すればいいですか？
Lisp風なリスト？
0861１３２人目の素数さん2005/04/02(土) 23:30:44
Mathematica自体がLisp風な訳だが。
0862１３２人目の素数さん2005/04/03(日) 11:01:02
>>861
でも、リストの深いところを破壊しようとすると全部書き換える事になって遅くなりませんか？
08638622005/04/03(日) 11:24:11
つうかリストの破壊が基本的にできない気がしてきました。書き換えしかない。
0864１３２人目の素数さん2005/04/03(日) 12:19:28
>>862
二分木のデータを変更すると、そこに至るための上の部分の木も再生成されるが、log(n)だ。
だから、普通Schemeで二分木使う時はset!でデータ書き換えたりしないと思うのだが。

Mathematicaでは「書き換え」はシンボルに対してしかできないだろう。
一部を書き換えると全部が「コピー」または「再生成」される。
どうしても書き換えたいなら、シンボルによる二分木を作っておいて、あとでシンボルに値を
アサインするようにすればいいんじゃないか。
08658622005/04/03(日) 22:19:37
>>864
そうですよね、「コピー」または「再生成」しかできない。
（ちょっとこっちな言語の用語に詳しくないもんで、すいません。）
Unique[]ってのがあるみたいなんで、それでノードを乱打して繋げてみます。
ありがとうございます。

#今日リストの大コピー合戦で実装してみたら、
#Length[Flatten[リスト]] = 10000あたりでもう超ウルトラ遅い。
#ノード数は200ぐらいなのにorz
0866１３２人目の素数さん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]
という解を得るには，どうしたらよいのでしょうか？
0867１３２人目の素数さん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]]
08688662005/04/07(木) 01:36:28
自己解決しました．867さん，ご意見ありがとう
ComplexExpand[Abs[a + b*I]]
でやってくれました
0869１３２人目の素数さん2005/04/08(金) 13:02:41
>>868
> 自己解決

自己解凍みたいw
0870博士 ◆OBpuRHyGa. 2005/04/08(金) 14:58:47
lim f(x)=sin∫f(x)となるようなf(x)をヒルベルト空間を用いて
x→∞
求めよ。バナハ空間でもよい。


0871１３２人目の素数さん2005/04/08(金) 15:17:58
>>870
積分の範囲は？
0872博士 ◆OBpuRHyGa. 2005/04/08(金) 16:12:06
0→∞
0873１３２人目の素数さん2005/04/08(金) 20:41:10
>>870

なんでヒルベルト空間なんているんだよ。
lim[n→∞]f(x)=0
∫[0,∞]f(x)dx=0
にすればいいだけじゃん。
たとえば
f[x]=(-1)^[x]/[1+x/2]
08748732005/04/08(金) 20:45:48
>>873
lim[n→∞]f(x)=0
↓
lim[x→∞]f(x)=0
0875１３２人目の素数さん2005/04/14(木) 15:41:02
>>868

MyRe[H_] :=
1/2*(ComplexExpand[H, TargetFunctions -> {Re, Im}] +
ComplexExpand[Conjugate[H], TargetFunctions -> {Re, Im}])
0876１３２人目の素数さん2005/04/24(日) 10:50:44
保守のため、上昇
0877１３２人目の素数さん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　になります。
0878１３２人目の素数さん2005/04/28(木) 19:14:20
>>877
Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals]でa - I bになるが。
@Mathematica 5.0
0879１３２人目の素数さん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版ですか？
08808782005/04/29(金) 19:29:24
>>879
Win版だけど、バージョンの違いだと思われ。OS変えるだけで
わざわざカーネルのソース書き換えないだろう。
0881１３２人目の素数さん2005/05/01(日) 10:16:29
Simplify[Conjugate[a + I b], a ∈ Reals && b ∈ Reals, ComplexityFunction -> LeafCount]

V4系にはConjugate[..]のほうが単純に見えるんだろう。
0882１３２人目の素数さん2005/05/03(火) 10:00:50
age
0883１３２人目の素数さん2005/05/10(火) 20:25:40
hos
0884１３２人目の素数さん2005/05/11(水) 19:13:01
ちょっと質問なのでが
Mathematicaを使ってアイパターンを描画させたいのですがMathematicaでできるでしょうか？
gnuplot使っていたのですがアイパターンを描けずMathematicaに移行したいと思っています
0885１３２人目の素数さん2005/05/11(水) 21:38:09
gnuplotでできなかったのは、誰にも聞かなかったからで、
Mathematicaでできたとしたら、それは誰かに教えてもらえたからでしょうね。
0886１３２人目の素数さん2005/05/12(木) 17:56:57
プログラム板でレスなかったのでこちらで質問させてください。

CUIで使ってるときに、1行に表示される文字数を変更する方法はありませんか？
例えば、

In[1]:= N[Pi,100]
Out[1]= 3.1415926535897932384626433832795028841971693993751058209749445923078\
> 16406286208998628034825342117068

としたときに、結果を改行せずに表示してほしいのですが。
0887１３２人目の素数さん2005/05/12(木) 18:51:26
Mathematica 5.01のMac版とSolaris版で試したのですが、
In[22]:= 1. 10^(-309) +0.0
Out[22]= 0.
上のようになるのはバグでなく機械精度の問題ですか？

また、この結果を
1. 10^(-309)
となるようにする方法ってありますか？
08888872005/05/12(木) 19:18:22
すいません
In[22]:= 1. 10^(-309) +0.0
Out[22]= 0.
こうなるのは、有効桁数のこと考えたら当たり前ですね。逝ってきまつ
0889１３２人目の素数さん2005/05/12(木) 21:03:13
>>886
無理
0890１３２人目の素数さん2005/05/15(日) 10:27:27
さっきランダムウォークの
コードを書いた。はろうか
0891１３２人目の素数さん2005/05/15(日) 10:57:05
今起きたばかりで
　コンドームと一瞬読めてしまった俺ｶﾞｲﾙ
0892１３２人目の素数さん2005/05/15(日) 16:55:26
迷ってるのなら、貼れ
0893１３２人目の素数さん2005/05/19(木) 22:01:32
>>892

今日また、コードを書いた。
0894１３２人目の素数さん2005/05/20(金) 18:44:08
Ver.5.1　９ラック版ある？
0895１３２人目の素数さん2005/05/20(金) 20:06:58
　　　　　　〃〃　ｶﾞｯ
　　　 ＿__〃　 ∧＿∧
　　 　|　　|　　（・∀・ ,,）
　　 .=|　　|==○○=　 ヽ
　　从|＿_|∩　ヽ　　､(⌒)っ
　 （　｀Д´）/　　 ヽ＿）
　　　　　　/ ←>>894
0896１３２人目の素数さん2005/05/21(土) 21:57:03
>>894
刑務所行って探して来い
0897１３２人目の素数さん2005/05/25(水) 18:10:48
おまいらMathematicaMark 2004の結果教えてくり
やり方は
<<Utilities`Benchmark`
BenchmarkReport[]

おいらはﾁﾝｺﾊﾟｯﾄﾞ T42 PenM 1.7GHz 512MB Mathematica5.1.1で1.59
0898１３２人目の素数さん2005/05/31(火) 12:33:44
３点O(0,0,0)、A(1,0,0)、B(0,1,1)を通る平面上に、
Aを中心とする半径１の円をｙ軸の周りに１回転してできる回転体を
Mathematicaを使って書きたいんですけど
どんなコマンドにしたらいいですか？
0899ChaosicSoul ◆/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 を左から掛ければいいだろう。
0900ChaosicSoul ◆/yaJbLAHGw 2005/05/31(火) 12:51:19
ちなみに正しいスペルは chaotic のようで。
0901８９８2005/05/31(火) 12:55:36
コピペでMathematicaに貼り付けたらできちゃうように
コマンドを作ってもらってもいいでしょうか…
0902ChaosicSoul ◆/yaJbLAHGw 2005/05/31(火) 13:02:29
Re:>>901 今Mathematicaが手元に無いから何とかして自分で作ってくれ。
0903８９８2005/05/31(火) 13:07:28
>902
ありがとうございました。

自分はあんまりMathematicaは使いこなせないんで、
ChaosicSoulさんがMathematicaを使えるようになったり、
他の方でも、できるようならコマンドを書いていただけるとうれしいです。
皆さんお願いします。
0904１３２人目の素数さん2005/05/31(火) 15:26:04
そこまで人に頼るのかよ
0905Zeno ◆5nZQbNmQPs 2005/05/31(火) 17:56:11

Reply-to:>>903
　　　　　　　　　　　　　 君は将来性がない人間だな。
　　　　　　　　　 そんな簡単なグラフは自分で書くようにしろ。

________________________________________________________
セオドア・カジンスキー 「産業社会とその未来」 (aka. マニフェスト) の抜粋
http://www.cnn.com/US/9509/unabomber/09-19/am/excerpt.html
0906１３２人目の素数さん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π}];
0907８９８2005/05/31(火) 21:37:46
>>905
英語ばっかりだ〜
読めん…

>>906
ありがとー！
やってみます。
0908１３２人目の素数さん2005/06/01(水) 00:12:28
学校の宿題を代わりにやってくれるスレはここですか
レス数が900を超えています。1000を超えると表示できなくなるよ。