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

〓 Mathematica 3 〓

■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん2005/07/15(金) 16:37:20
　　　　　　　　_....._{{ 〃
　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-<
　 　 ,' 　　.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/
0072GiantLeaves ◆6fN.Sojv5w 2005/08/22(月) 21:41:33
talk:>>71 Flatとかいう関数が無かった？あるいはPart関数とリストに要素を追加する関数を組み合わせるか。
0073GiantLeaves ◆6fN.Sojv5w 2005/08/22(月) 21:44:02
talk:>>71 http://documents.wolfram.co.jp/v5/Tour/MathematicaAsAProgrammingLanguage.ja.html
0074１３２人目の素数さん2005/08/23(火) 10:13:35
ありがとうございました
0075１３２人目の素数さん2005/08/25(木) 13:57:06
基本的なことかもしれませんが・・・
x=n(n=const)
のグラフはどのように描けばよいのでしょうか？
Plot[x = 3, {x, 0, 5}]
としてみたのですがエラーが出てy=3が描画されてしまいます
0076１３２人目の素数さん2005/08/26(金) 21:54:12
ParametricPlot[{x=3,y=t},{t,0,5}]
こんな感じで、いかがでしょう?
0077１３２人目の素数さん2005/08/30(火) 19:20:33
うまいね
0078１３２人目の素数さん2005/08/31(水) 00:29:31
age
0079１３２人目の素数さん2005/09/08(木) 21:52:47
sage
0080１３２人目の素数さん2005/09/09(金) 00:25:59
age
0081１３２人目の素数さん2005/09/09(金) 02:34:09
ﾈﾀくれ
0082１３２人目の素数さん2005/09/09(金) 12:12:43
3次元で r = θ の極座標を描きたい。
0083１３２人目の素数さん2005/09/10(土) 01:23:23
>>82
SphericalPlot3D?
0084１３２人目の素数さん2005/09/16(金) 12:53:36
分散共分散行列を1発で求める関数はありますでしょうか？
0085１３２人目の素数さん2005/09/16(金) 16:44:11
CovarianceMatrix or CovarianceMatrixMLE
0086１３２人目の素数さん2005/09/18(日) 16:31:55
質問させてください
1×n行列を転置してn×1行列にしたいのですが,
Transpose[matrix]
を使ってもエラーが出てできませんでした
どうすればよいのでしょうか？
0087１３２人目の素数さん2005/09/18(日) 18:36:47
In[1]:=Transpose[{{1,2,3}}]
Out[1]={{1},{2},{3}}
0088１３２人目の素数さん2005/09/18(日) 22:15:06
>>87
ありがとうございました
参考になりました

もう1つ質問させてください
リストの正数や負数の数をカウントしたいのですがどうすればよいのでしょうか？
Count[list , _>0]やCount[list , # > 0 &]
としてみたのですがうまくいきませんでした
0089１３２人目の素数さん2005/09/19(月) 08:06:03
Count[{0, 1, -2, 3, 4, -5}, _?Positive]
Count[{0, 1, -2, 3, 4, -5}, _?(# > 0 &)]
0090１３２人目の素数さん2005/09/19(月) 19:48:04
ある計算とその結果を出力させる関数1を別の関数2から実行するとき,
一時的に関数1の出力をなくすようなことはできますか？
0091１３２人目の素数さん2005/09/20(火) 09:04:00
age
0092１３２人目の素数さん2005/09/30(金) 02:57:30
mage
0093１３２人目の素数さん2005/09/30(金) 13:25:36
agenai
0094１３２人目の素数さん2005/10/02(日) 14:20:05
今4.1を使っているのですが、5.2にバージョンアップするだけの
価値はありますか？
0095１３２人目の素数さん2005/10/02(日) 15:29:26
公式ウェブページ見て、便利そうと思わなければ特に変えなくてもいいと思う。

地味なところでは計算精度・速度の向上もあるが、100*100の行列の演算くらいの大規模計算で無い限り差が分からないと思う。
0096942005/10/03(月) 17:56:08
>>95
ありがとうございます。
あまり利点がなさそうなのでうｐは当分見送ることにします。
0097１３２人目の素数さん2005/10/04(火) 00:42:27
>>96
取り敢ず体験版使ってみれば？
日本語化されてるし。
0098942005/10/05(水) 04:17:31
>>97
使ってみました。
特に面白そうなのはなかったんですがReduceだけは便利そうですね。
Sin(x)==Cos(x)みたいなのが解けるっていうのにはちょっと惹かれました。
0099１３２人目の素数さん2005/10/07(金) 12:38:29
リストをランダムに並び替えるにはどのようにしたら良いでしょうか？
アドバイスよろしくお願いします
0100１３２人目の素数さん2005/10/07(金) 16:18:44
NormalPValueは平均0分散1のときのp値を求めることができるようですが，
任意の平均，分散の正規分布のp値を求めたいときはどうすればいいですか？
0101１３２人目の素数さん2005/10/15(土) 06:41:18
　
0102１３２人目の素数さん2005/10/18(火) 17:11:12
AxesLabelで軸のラベルを書くとy軸のラベルはy軸と直角に表示されますが，
これを回転してy軸と水平に表示させる方法はありますか？
0103１３２人目の素数さん2005/10/19(水) 18:54:43
ListPlotで、離散的な関数を図示しているのですが、
Frame-> Trueの際にTicksを有効にしたいのですが、
なにか名案があったらお教え下さい。
0104１３２人目の素数さん2005/10/20(木) 20:00:43
age
0105１３２人目の素数さん2005/10/23(日) 23:52:15
N-joy the noiZ


sR=441*4;siZe=sR*75;timeLine=Range[siZe];t=n/siZe;
a=0.41+0.59*Exp[-3*(2t-1)^2]*(Sin[Pi*((6*Abs[t-0.5])!)])^2; mL=5;mR=11;x0=0.98;
yL=Cos[mL*ArcCos[x]];fL=Compile[{x,n},Evaluate[x+a*(yL-x)]];
dL=FoldList[fL,x0,timeLine];
ListPlot[dL,PlotStyle->{RGBColor[0.05,0.2,0.85]}];
yR=Cos[mR*ArcCos[x]];fR=Compile[{x,n},Evaluate[x+a*(yR-x)]];
dR=FoldList[fR,x0,timeLine];
ListPlot[dR,PlotStyle->{RGBColor[0.85,0.2,0.1]}];
ListPlay[{dL,0.5*dR},SampleRate->sR,SampleDepth->16]
0106１３２人目の素数さん2005/10/25(火) 02:26:22
ﾜﾚﾏﾃｨｶはどこで拾えますか？（MathID とか聞いてこないやつ）
0107１３２人目の素数さん2005/10/25(火) 18:34:26
氏ね
0108１３２人目の素数さん2005/10/25(火) 19:42:44
反応するなよ
0109１３２人目の素数さん2005/10/25(火) 19:57:35
みんな市んじまえ
0110１３２人目の素数さん2005/10/26(水) 10:35:52
>>109
どうした。悩みがあるなら聞くぞ？
0111１３２人目の素数さん2005/10/26(水) 20:57:03
>>99>>100>>102>>103
に答えてくれないかな、いや、ください。
0112１３２人目の素数さん2005/10/29(土) 16:59:19
x^2 + y^2 = r^2
を
x = rcosθ
y = rsinθ
みたいに媒介変数表示に変換したいです。
上記のような単純な円だといいのですが、
任意の楕円
ax^2 + bxy + cy^2 + dx + ey = 1
とかから媒介変数に変換するにはどうしたらいいんでしょうか？
そもそもmathematicaでできますか？
0113１３２人目の素数さん2005/10/30(日) 01:16:21
>>111
　>>99
　　RandomRelist[x_List] := Block[{n = x, p}, Do[p = Random[Integer, {1, i}]; n[[{p, i}]] = n[[{i, p}]], {i, Length[x]}]; n]
　>>100
　　HypothesisTests.mを読む
　>>102
　　RotateLabel
0114１３２人目の素数さん2005/10/31(月) 05:23:47
age
0115１３２人目の素数さん2005/11/01(火) 12:00:22
ver. 3 くそ
ver. 4 本来３だったもの
ver. 4.1 バグ取り
ver. 4.2 新たなバグ付き
ver. 5 どうでもいい
ver. 5.1 バグ取り
ver. 5.2 新たなバグ付き
ver. 6.0 さらなるバグ付き
0116１３２人目の素数さん2005/11/01(火) 20:49:47
でﾜﾚﾏﾃｨｶは？
0117１３２人目の素数さん2005/11/04(金) 22:35:23
放流するか
0118１３２人目の素数さん2005/11/04(金) 23:08:31
ny使い方わからん
0119１３２人目の素数さん2005/11/05(土) 18:30:42
>>117
Student版買ったのに・・・
登録とかしなくて使えんの？
0120１３２人目の素数さん2005/11/06(日) 01:53:51
nyはいかん．矢風億がいい．
0121１３２人目の素数さん2005/11/06(日) 06:20:31
>>120
億漁ったけどイマイチだな・・・
以前なら自サイトへのリンクとかあって色んなソフト販売してる奴らが多かったから
直接聞くことができたが、今は取り締まりが厳しいのかな。
ライムやっても罠入りばっかだしな・・・。
0122名無しさん＠Mathematica2005/11/06(日) 13:04:32
>>117
お使いのはVersion4や4.1とか4.2でしょうか
もしそうでしたら[email protected]へご連絡いただければ
便利なつーるをおわたしします
0123１３２人目の素数さん2005/11/07(月) 01:29:29
変数の値が、0.000000000000001
だった場合に、これを 0.000000000000001　
と出力させる方法はないですか？？
それか、1*10^(-14) と出力する方法はないでしょうか。

mathematica　の出力を 他のプログラムに食わせることを考えてます。
CForm が期待してるものに近いんだけど、そのまま使えないです。いまは、、CForm の出力に、
ラッパーかませてます。おしえてくださいませ。
0124１３２人目の素数さん2005/11/07(月) 07:26:44
5.1の日本語が使える鍵genはまだぁ〜
0125１３２人目の素数さん2005/11/07(月) 13:55:24
>>123
0.000000000000001 だったら

fooForm[n_] := NumberForm[n, ExponentFunction -> (Null &)]

1*10^(-14) なら

barForm[n_] := NumberForm[n,
NumberFormat -> (SequenceForm[#1, "*", #2, "^", #3] &),
ExponentFunction -> Identity
]
0126１３２人目の素数さん2005/11/08(火) 01:20:26
>>124
汗水垂らして働いて買えよ。
0127１３２人目の素数さん2005/11/09(水) 18:52:43
鼻水垂らして割れ物あさってまつ
0128１３２人目の素数さん2005/11/10(木) 19:35:06
NDsolveで得られた補完関数ってどうやって微分するのでしょうか？
0129１３２人目の素数さん2005/11/14(月) 19:11:11
D
0130１３２人目の素数さん2005/11/16(水) 22:11:26
http://www.pcsoftoff.com/2bu.htm
Mathematica5.1大安売り
0131１３２人目の素数さん2005/11/17(木) 00:08:24
>>130
IPアドレス 219.101.229.147
ホスト名 www.pcsoftoff.com
IPアドレス
　割当国　※ 日本 (JP)

国内IPか。へぇ。
0132１３２人目の素数さん2005/11/17(木) 02:34:42
>>130
ACCSに通報しますた。
0133１３２人目の素数さん2005/11/17(木) 13:59:11
日本語も使えるのかな？
多分使えないモノだなwww
0134１３２人目の素数さん2005/11/17(木) 16:04:50
>>130
高杉。
0135１３２人目の素数さん2005/11/18(金) 17:43:32
>130
　日本語も使えます。64ビットのOSなんか使用している人は、いるんですかね。
0136１３２人目の素数さん2005/11/20(日) 08:26:07
日本語使えるって、英語版に入力できるだけでしょ
メニューやヘルプは英語
日本語版の管理は厳重
Mathematica万歳
0137１３２人目の素数さん2005/11/20(日) 08:53:14
まてまてかっていくらくらいで買えるの？
0138１３２人目の素数さん2005/11/20(日) 20:23:54
>>137
ttp://www.amazon.co.jp/exec/obidos/ASIN/B000B6GGC4/
0139１３２人目の素数さん2005/11/24(木) 14:58:37
全然わからないのですいませんが質問させてくださいm(_ _)m

p1(r,z)=∫exp[(j*k*r'^2)/2*(1/z-1/D)]*J0((k*r*r')/z)*r'dr'

先生にこの式を鳥瞰図にしてくれって頼まれたのですが
もうずっとわからんままで、先生に聞いても教えてくれんのです(ﾉД`)

j=虚数、k=1500π、D=45,J0=0次ベッセル関数,r=0〜20,z=20〜60
積分は定積分で、範囲は0〜17です。
あ、ソフトはMathcadとMuPADとMathematica5.1とかいろいろためしてます
Mathematica5.1は英語であんまよくわからんのですけどね(´･ω･｀)
どなたかわかる方いたら助けてください（ノ＿・。）
01401392005/11/25(金) 18:14:25
すいませんなんとかできたっぽいです（；´Д｀）
0141１３２人目の素数さん2005/11/25(金) 18:56:49
>>140
BesselJ出すのに時間かかったんじゃね？
後学の為にMathematicaの式コピペして
0142１３２人目の素数さん2005/11/26(土) 07:39:16
age
0143１３２人目の素数さん2005/11/28(月) 22:02:53
>>130
高杉て話にナラナイ。そんな値段なら学割で正規品を買った方がマシかな。
0144１３２人目の素数さん2005/12/01(木) 21:36:18
世紀版の方が高いだろ
0145１３２人目の素数さん2005/12/02(金) 18:21:32
読み込みんだビットマップ画像のサイズを縮小することはできるでしょうか？
50×50ピクセル　->　25×25ピクセルのようにしたいのですが。。。
0146１３２人目の素数さん2005/12/04(日) 09:48:10
画像ソフトでやれ
0147１３２人目の素数さん2005/12/06(火) 23:00:26
文字式で
(sqrt{A*B} + sqrt{A*C}) / A
みたいなやつを、分母分子を√Aで割って
(sqrt{B} + sqrt{C}) / sqrt{A}
に変形して表示したい場合はどうすればいいんでしょう？

SimplifyもFullSimplifyも無視されてしまいます、、、
0148１３２人目の素数さん2005/12/07(水) 00:16:02
>>147
できたよ
\!\(f = \(\@\(A\ B\) + \@\(A\ C\)\)\/A\)
\!\(PowerExpand[f] /. {\@A -> A\/\@A} // FullSimplify\)
こういう時は新しいRule［→］を作ってやるとうまくいく
01491472005/12/07(水) 00:20:47
ごめん･･･PowerExpand[(sqrt{A*B} + sqrt{A*C}) / A] // Simplifyでも逝けるじゃん･･･吊ってくる　orz
01501472005/12/07(水) 02:27:21
>>148=149 (ですよね？)

御蔭様で、どちらの方法でも完璧にできましたです。
即レスありがとうございました！
0151１３２人目の素数さん2005/12/07(水) 19:58:19
age
0152１３２人目の素数さん2005/12/08(木) 20:55:35
われまてぃか
万歳
0153１３２人目の素数さん2005/12/09(金) 01:10:39
どなたか
http://science4.2ch.net/test/read.cgi/math/1133617862/441
の計算をよろしくお願いします
無知識の者が免許の学科試験に合格する確率です
0154１３２人目の素数さん2005/12/09(金) 12:16:04
ttp://www.uploda.org/file/uporg256633.jpg
これでおｋか？つーか海外にはガウス記号って言葉無いんだな('A`)ハジメテシッタ
ttp://blog.livedoor.jp/calc/archives/717594.html
0155１３２人目の素数さん2005/12/09(金) 13:44:06
>>154
mathematicaでのガウス関数と等価なのはFloorだよ（床関数）
それでやったらすんなりいくと思われ
再うｐキボンヌ
0156１３２人目の素数さん2005/12/10(土) 02:10:17
>>155
ドモ(・ω・)ノ
うちの研究室じゃガウス記号なんて絶対使わないから初めて知ったよｗ
ttp://www.uploda.org/file/uporg257191.jpg

＞無知識の者が免許の学科試験に合格する確率です
絶対受かんないなｗソース希望
0157１３２人目の素数さん2005/12/10(土) 18:53:48
ファイルにカンマ区切りで数字を書き出す際、
　 　 　 -6
0.1,1.3 10
みたいになるのをふせぐにはどうしたらいいですか？
0.1,1.3E-6
と扱いやすい形で保存したいです。
0158１３２人目の素数さん2005/12/10(土) 19:56:39
>>156
どうもありがとう(・ω・)ノ
確かに絶対受からないねｗ
-19乗オーダーなんて物理とかでしかみたことねーや（電気素量は確かこれぐらいのオーダーだったような）

ソースは441のすぐ上だよ
http://science4.2ch.net/test/read.cgi/math/1133617862/412n-
0159１３２人目の素数さん2005/12/13(火) 19:43:54
ListPlotでプロットした際、OlotJoinedをTrueにすると各点が自動的に結ばれますが、
任意の2点を結ぶことはできるのでしょうか？
例えば、
{{x1,y1},{x2,y2},{x3,y3},{x4,y4}}
をプロットしたときに
{x1,y1}と{x2,y2}を、また{x3,y3}と{x4,y4}を結んで
{x2,y2}と{x3,y3}の間は結ばれていない
・・・というのは無理でしょうか？
0160GiantLeaves ◆6fN.Sojv5w 2005/12/13(火) 20:49:59
talk:>>159 線分を描く関数があった気がするから、探してみよう。
01611592005/12/14(水) 12:44:58
>>160
ParametricPlotを使って書くことができました
ありがとうございました
01621592005/12/15(木) 16:55:55
またまた失礼します
x軸にはx1,x2に，y軸には自動で破線の目盛線を入れようと思い、ヘルプの通り

sty1=Dashing[{0.05,0.05}]
Plot[... , GridLines->{{{x1,sty1},{x2,sty1}},{Automatic,sty1}}]

としたのですが，y軸の目盛は破線になっていたのですが，
x軸の目盛はデフォルトの青線でも指定した破線でもなく
黒の実線で表示されてしまいます。

どうすればよいのでしょうか？
よろしくお願いします
0163１３２人目の素数さん2005/12/15(木) 22:04:04
がんがれ!
0164１３２人目の素数さん2005/12/28(水) 15:26:34
NDSolveで微分方程式を解くとき、独立変数の区間を指定して解きますが、
従属変数の区間を設定する（もしくはある値を超えたら計算終了とか）
ことは出来ますか？（無限ループになる可能性もあるわけですが...）
0165１３２人目の素数さん2005/12/31(土) 12:31:32
気合いでできる
Mathematicaに不可能はない!
0166１３２人目の素数さん2005/12/31(土) 14:32:03
【かっこ悪い】建部崩れ、見参！【情けないｗ】
http://science4.2ch.net/test/read.cgi/math/1135765594
【夢vs】結果を出せば職はある？ｗ【現実】
http://science4.2ch.net/test/read.cgi/math/1134888899
【事実】研究しても、ポスト無し！【愕然】
http://science4.2ch.net/test/read.cgi/math/1134089493
関連：【建部 】斎藤毅先生【Invent】
http://science4.2ch.net/test/read.cgi/math/1134743220
0167１３２人目の素数さん2006/01/02(月) 05:14:25
884
0168１３２人目の素数さん2006/01/13(金) 14:23:41
例えばx^2+y^2<=1の面積を出したいときy=f(x)の形にせずに出すことって
できますか？
x^3+y^3<=3xyとかx^2+y^2+sinx+siny<=1のようなy=f(x)にしずらいような
図形の面積とか出してみたいんですけど。
0169１３２人目の素数さん2006/01/13(金) 23:37:03
>>168
<< Calculus`Integration`;
NIntegrate[Boole[x^2+y^2<=1],{x,-∞,∞},{y,-∞,∞}]
NIntegrate[Boole[x^3+y^3<=3x y],{x,0,∞},{y,0,∞}]
こんなのとか。NIntegrateでなくIntegrateの場合、ちゃんと計算して
くれない時がある。それから、定義域はある程度自分で面倒を見ないと
いけない。
0170１３２人目の素数さん2006/01/19(木) 14:04:47
Integrate[Boole[x^2+y^2+z^2+Sin[4x]+Sin[4y]+Sin[4z]<=1],{x,-∞,∞},{y,-∞,∞},
{z,-∞,∞}]
これをやったら結果は0になった・・・。どうみても体積があるようにみえるんだが・・・。
0171１３２人目の素数さん2006/01/20(金) 17:55:08
ある関数を定義して，その実行中になんらかのエラーが発生した場合
処理を途中で止めたいのですがどうすればよいでしょうか？
エラーが発生した場合の戻り値かなにかがあればIfとBreakまたはAbortでできそうなのですが。。。
■ このスレッドは過去ログ倉庫に格納されています