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

〓 Mathematica 3 〓

レス数が1000を超えています。これ以上書き込みはできません。
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/
0951１３２人目の素数さん2007/11/14(水) 14:46:17
その教育機関の学生でないとライセンス違反というだけ。
プログラムに制限はないよ。
0952１３２人目の素数さん2007/11/14(水) 22:19:59
>>948
If[条件，a；b；c；｛a，b，c｝，d；e；f；｛d，e，f｝]

TRUE なら{a，b，c}
FALSE なら{d，e，f}
というふうにリスト形式で出力される。

>>950
Studentバージョンの場合、OS再インストールするだけで同じPCでも再認証が必要。
学校の住所、学校の電話番号、学校名、ユーザーの氏名、再認証理由、等をWEB上で書かなければならない。
再認証を何回も繰り返していると、学生証明をFAXで求められる。この場合WEB認証は一時的に不可能になる。

うざいことこの上ないので、バイナリエディタで（ｒｙ
0953１３２人目の素数さん2007/11/15(木) 00:04:14
デスクトップとは別に
notePCで使う場合
また金払わないとダメなんじゃん....
知らんかった
0954１３２人目の素数さん2007/11/15(木) 00:05:42
ここまで厳格なアカデミックライセンスも凄いね
0955１３２人目の素数さん2007/11/15(木) 00:30:59
マイクロソフト様でも
デスクトップとノートの２台にインストールすることを認めてくれているのに
0956１３２人目の素数さん2007/11/16(金) 14:46:19
ver6.0（Win）にGTensorM入れてOptionsの内容もちゃんと調整したんだが
たとえばqload[SCHW]でシュバルツシルト計量読ませたら

　xSuperscript[1] = r, xSuperscript[2] = th, xSuperscript[3] = ph, xSuperscript[4] = t

みたいに出て、本来表示されるはずのシュバルツシルト解の数式が

　ds Superscript[2] = ds[SCHW]

って出てきてろくに使い物にならないんだけど、何かもうちょっと設定弄らなきゃダメかな？
卒研で使うことになって、研究室側が買うので入れろってことで初めて入れて今色々弄ってるんだけど。
研究室のPC（ver5）の方で動いてるのは確認したんだが…もしかして6で使えなくなった？
0957１３２人目の素数さん2007/11/16(金) 15:03:33
間違えた…GRTensorMね。
0958１３２人目の素数さん2007/11/17(土) 01:06:33
どうせ仕様に無い機能使ってプログラムが書いてあるんだろ。
仕様に無い機能は、版が変わって振る舞いが変わっても文句はいえない。
0959１３２人目の素数さん2007/11/17(土) 23:42:32
やっぱver6 はだめか。
0960１３２人目の素数さん2007/11/18(日) 01:38:00
>>958
あー、やっぱそんなとこか…。
0961１３２人目の素数さん2007/11/18(日) 11:37:35
根拠0の憶測を信じるなよ
0962１３２人目の素数さん2007/11/22(木) 09:24:10
y(x,t)=A sin(kx-ωt)のアニメーションはどのようにやればいいのでしょうか？
基本的な質問でスマン
09631stVirtue ◆.NHnubyYck 2007/11/22(木) 12:27:21
Reply:>>962 Javaでなんとかしろ。
0964１３２人目の素数さん2007/11/28(水) 08:35:27
y[x_,t_]:=a Sin[k x-\[Omega] t];
Animate[Plot[y[x,t]/.{a->1,k->1,\[Omega]->2\[Pi]},{x,0,10}],{t,0,10}]

あとは適当にいじれ、ただし6.0じゃないと動かんぞ
0965１３２人目の素数さん2007/11/28(水) 23:53:53
横軸は同じで、縦軸のスケールが違う２つのグラフ(一方は10000程度、もう一つは20程度)を重ねたいのですが、
どのようにすればいいのでしょうか？
Showだとスケールの同じものしか出来ませんでした
09661stVirtue ◆.NHnubyYck 2007/11/29(木) 00:13:08
よくみたら Mathematica のスレッドだった。
0967１３２人目の素数さん2007/12/02(日) 06:30:02
成分数任意のベクトル関数演算を定義することは可能ですか？
0968１３２人目の素数さん2007/12/02(日) 12:31:09
すみません、補足します。

f[指定された多変数]=（任意の多変数からなるスカラー関数）

というタイプの数式を仮定して、入力式に対して「指定された多変数」主体の
Gradや内積をとり、また「任意の多変数」に実際の数値も入力して結果を
出力するようなサブルーチンを作りたいのですが、この場合入力式とは別に

・数式に用いた変数のリスト
・演算に使う変数のリスト
・数式に入力する定数のリスト

を同時に入力させて、リストの長さ、要素の取得などを応用した反復演算などで
地道に組んでいくしかないんでしょうか？
0969１３２人目の素数さん2007/12/04(火) 05:24:16
基本的な質問ですみません。

data = {{0, 0}, {1, 10}, {2, 20}, {3, 30}, {4, 100}, {5, 200}};

gr1 = ListPlot[data,
PlotStyle -> {RGBColor[1, 0, 0]}]

上のgr1のグラフと、

Do[Plot[a x ( x - 1), {x, 0, 10},
PlotStyle -> {RGBColor[0, 0, 1]}], {a, 1, 5, 1}]

このようなアニメーションのグラフを重ねて出すには
どうすれば良いのでしょうか？＾＾；
よろしくお願いします。
0970１３２人目の素数さん2007/12/04(火) 14:48:02
このような３元連立方程式を解きたいのですが、どのようにすればいいでしょうか。
知恵をお貸しください。

(a - 234.25 + b*c)*(1/1.01713 - 1/1.00036) -
c*(115/1.01713 - 156/1.00036) == 500,
(a - 234.25 + b*c)*(1/1.00973 - 1/1.00077) -
c*(119/1.00973 - 154/1.00077) == 500,
(a - 234.25 + b*c)*(1/1.01027 - 1/1.00005) -
c*(115/1.01027 - 157/1.00005) == 500
0971１３２人目の素数さん2007/12/06(木) 21:17:44
>>969
これでどう？

Animate[Show[
Plot[a x (x - 1), {x, 0, 10}, PlotRange -> {{0, 10}, {0, 200}},
PlotStyle -> {RGBColor[0, 0, 1]}], gr1], {a, 1, 5, 1}]
0972１３２人目の素数さん2007/12/06(木) 21:43:44
一本買ったら、
自宅用ライセンスというのを
もう一本もらった

自宅でも仕事しろ　ってか
0973１３２人目の素数さん2007/12/07(金) 12:55:22
u1=[ [ { arctan( ( 1-k ) * cot( u1*v*w ) -b*tan(u1*v*w) ) } /( 2*a) ] + ( q*π ) ] / [ { arctan( (1-k) * cot(v*w) -b*tan(v*w) ) } / (2*a) ]
u2=[ [ { arctan( ( 1-k ) * cot( u2*v*w ) -b*tan(u2*v*w) ) } /( 2*a) ] + ( r*π ) ] / [ { arctan( (1-k) * cot(v*w) -b*tan(v*w) ) } / (2*a) ]

上の2式において

a=sqrt(k)*(y/x-x/y)
b=k*(y/x)^2-(x/y)^2

mathematica流の書き方じゃなくてすみません。
これら2つの非線型方程式からxとwを数値解として求めるプログラムをつくりたいのですが・・・

x,w 以外の値はすべて既知で、
u1,u2 は 1<u1<u2<=3 である任意の実数
q,rは 0または1
それ以外は任意の正の実数
分母は必ず正の値(になるはず)
最適解を求める際に、分子のarctan(・・・)の項が負の値になってしまう時にはqまたはrを「1」にしてπを加える
それ以外の場合にはq,rは「0」

と、こんな感じなのですが,どなたかサクっと作ってくださいませんでしょうか・・・
0974１３２人目の素数さん2007/12/07(金) 13:26:21
おまえ文系だろ
0975１３２人目の素数さん2007/12/08(土) 01:31:22
>>965
参考に。
http://support.wolfram.com/mathematica/graphics/2d/twoaxisgraph.ja.html
0976１３２人目の素数さん2007/12/08(土) 04:03:07
質問です。

aa={-3,-1,1,2,3,5,2,1,-2,-4,-4,-3,1,5,1,-1,-2}

のようなリストがあるとします。
このaaから

cc={2,8,12,15}

のようなリストを作りたいのですが、これをスマートに行う計算方法を探しています。
現在は、

bb=Table[If[aa[[n]]*aa[[n+1]]<0,n,neednot],{n,1,Length[aa]-1}];
cc=Select[bb,NumberQ];

という感じで処理して計算できてはいるのですが、ちょっと無理やり感があるので、
もっと直接的に行えればよいのですが…
0977１３２人目の素数さん2007/12/08(土) 04:51:48
>>976
符号の切り替わる場所のリスト、でいいんだよね？
短くなったけどどうだろうこれ

cc = Join @@ Position[Sign /@ Times @@@ Transpose[{Drop[aa, -1], Drop[aa, 1]}], -1]
0978１３２人目の素数さん2007/12/08(土) 09:46:20
たとえば
x^(2/3) (x + (x^2) x^(7/3))
に
x=k y^3　(aは任意の実数)
を代入して全ての括弧を取り払いxの関数に展開したいのですが、どうすればいいでしょうか？
0979１３２人目の素数さん2007/12/08(土) 09:49:37
yの関数です、すみません。
0980１３２人目の素数さん2007/12/08(土) 10:17:46
>>978
まず、
expr=x^(2/3) (x + (x^2) x^(7/3)) /. x -> k y^3

expr2=Collect[expr,y]

でだめなの？
もし、kもyも実数でいいなら
Collect[expr, y,
Simplify[#,
Assumptions -> k \[Element] Reals && y \[Element] Reals] &]
とすれば簡単になる。aじゃなくてkだよね？
0981１３２人目の素数さん2007/12/08(土) 15:24:57
>>977
ありがとうございます！
まだまだ素人なのでちょっと見た感じでは文法が理解できないですが；
週明けに早速試してみようと思います。
0982１３２人目の素数さん2007/12/08(土) 15:40:00
>>976-977
0が入っている場合があるのか、あったらどうするのかがちょっと気になる
0983１３２人目の素数さん2007/12/08(土) 16:10:36
>>980
だめみたいです。
k^5 y^15 + k y^3 (k y^3)^(2/3)
　　　　　　　　　　　　　↑
　　　　　これが残ってしまって。

次の目的として最高次数を取得したいので、正確に評価できるよう
この部分も開いてしまいたいのですが…。
0984１３２人目の素数さん2007/12/08(土) 16:37:20
二年百四十六日。
0985１３２人目の素数さん2007/12/08(土) 18:14:48
>>982
えっと、>>976ではわかりやすいよう簡単に書きましたが、
実際のaaにあたる部分は特殊関数を用いた関数をサンプリングしたものなので、
ぴったり0になることは考えていないです。
しかもその関数はただの特殊関数ではなく、色々と複雑な処理を行っている
（さらにパラメータふたつを変えながらTableで回す）ので、
根の近似値を初期値としてFindRootで解くということが（現状の私の技術では）できないのです。
この方法で解いた値を初期値として解く事は勿論可能でしょうが、現状以上に膨大な時間が
かかることになる＆そこまでするほど厳密な値を必要としていないので。
09869802007/12/08(土) 20:22:12
>>983
Collect[expr, y,Simplify[#,Assumptions -> y>0] &]
だと
k^(5/3) y^5 + k^5 y^15
となるね。
Assumptionsでyが実数としてもy<0でもだめみたいだね。
どういうことなのかいまひとつよくわからない、、、。
09879802007/12/08(土) 20:30:04
追伸。
Simplify[expr, Assumptions -> y > 0]
だけでそうなる。Collectは必要ない、、。
0988１３２人目の素数さん2007/12/08(土) 20:42:03
>>987
ありがとうございます。
変数は当面は正で問題ないのでそれでやってみます。

バリバリ数式変形してから結果を入れたいときには鬱陶しいけど
数値入力する時のためにうかつに開けないようにしてるのかな…？
0989１３２人目の素数さん2007/12/08(土) 23:17:27
>>956
遅レスだが
http://support.wolfram.co.jp/applicationpacks/general/apppacksinv6.html
から引用「バージョン6.0より前のアプリケーションパッケージには
Mathematica 6.0との互換性はないとお考えください．お試しいただくことは
できますが，多くの機能が以前のようには動作しなくなっています．」

wolfram researchは何を考えてるのかね。
0990１３２人目の素数さん2007/12/08(土) 23:24:53
MS-DOSとWindowsのことを考えています
0991１３２人目の素数さん2007/12/09(日) 00:12:49
>>990
えらく古そうだけど何の話？
0992１３２人目の素数さん2007/12/09(日) 09:07:50
次スレおね
0993１３２人目の素数さん2007/12/09(日) 16:37:21
二年百四十七日。
0994１３２人目の素数さん2007/12/09(日) 20:42:30
スレタイが「Mathematica 4」だとバージョン4限定のスレに見えてしまい気になるます
0995１３２人目の素数さん2007/12/09(日) 21:00:13
それじゃ
Mathematica Ver.6以下はクソ　4
にしろ
0996１３２人目の素数さん2007/12/10(月) 16:10:00
〓 Mathematica 四 〓
http://science6.2ch.net/test/read.cgi/math/1197270001/
0997１３２人目の素数さん2007/12/10(月) 16:37:20
二年百四十八日。
0998１３２人目の素数さん2007/12/11(火) 05:36:23
.
0999１３２人目の素数さん2007/12/11(火) 05:46:38
999
1000１３２人目の素数さん2007/12/11(火) 05:47:13
1000
1001１００１Over 1000Thread
このスレッドは１０００を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。
レス数が1000を超えています。これ以上書き込みはできません。