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

〓 Mathematica 3 〓

レス数が950を超えています。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/
0896１３２人目の素数さん2007/09/29(土) 23:15:01
こんにちは、

y=625.733*((2/5)*(1-x)*a2^2-(4/105)*(1+2*x)*a2^3);
があります。
このとき、a2は
a2->-((7*(-1+x))/(1+2*x))-s1;
です。

ｓ１は
nが１のときは、ｓ１＝０
ｎが２以上のときは、s1=(92*n)/(y*10^6)
となります。

ｎを増やして、ｙが5.9に一番、近づいたときの
ｎを計算したいのですが、どのように計算すれば
良いでしょうか？
下記は、ｎを増やして、数値計算で求めたものです。
直接計算する方法を教えて下さい。
x=0.767476;
For[n=1,n<2*10^4,n++,
If[n==1,s1=0];
If[n>1,s1=(92*n)/(y1*10^6)];
y=625.733*((2/5)*(1-x)*a2^2-(4/105)*(1+2*x)*a2^3);
y1=y/.a2->-((7*(-1+x))/(1+2*x))-s1;
If[y1-5.9<0,Print["n=",n]];
If[y1-5.9<0,Print["y1=",y1]];
];

0897１３２人目の素数さん2007/09/30(日) 21:09:07
＞＞896
Maple でやってもらうと
問題の式を、x, y, n の関数として
f:= (x,y,n)-> y-(625.733*((2/5)*(1-x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^2-(4/105)*(1+2*x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^3));
で定義して、f( 0.767476,5.9,n)=0 の解 n を求めると n=13820.40245、59148.69709 となりました。

ということです。Mathematicaでは

x=0.767476;
z=Solve[y-(625.733*((2/5)*(1-x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^2-(4/105)*(1+2*x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^3))-5.9Š0,n];
Simplify[z]

では計算できませんが、どこが悪いのでしょうか？


0898１３２人目の素数さん2007/10/01(月) 20:56:01
>897

失礼しました。下記で求まりました。

x=0.767476;
z=Solve[y-(625.733*((2/5)*(1-x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^2-(4/105)*(1+2*x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^3))==0,n];
y=5.9
Simplify[z]
0899１３２人目の素数さん2007/10/02(火) 21:02:09
>898
やはり、解が一致しませんでした。
y=5.9
のときは、近い値が得られるのですが、
ｙ＝0.1
では答えが異なります。
如何でしょうか？
0900１３２人目の素数さん2007/10/03(水) 09:57:34
よく分からないけど
先に y=5.9 ってやってから Solve した方がよさそうですよ
>>898のだと丸め誤差の影響か虚数部がひどいことになりました

y=0.1のときは
{{n -> -347.016}, {n -> 651.842}, {n -> 742.06}}
ではだめなの？
0901１３２人目の素数さん2007/10/03(水) 20:46:10
>>900
お返事ありがとうございます。
>ではだめなの？
すいません。駄目なのです。
質問の本質は変更しませんが、
複雑な式は、分かりにくいだけなので止めて、
式を下記に変更したいと思います。
下記が、“０”に一番近づくのは
“７”付近なのですが、
このように数値計算をして求めるのではなく
式（たぶん微分方程式）を
解いて、下記と同様な解を
得るには、どのような式（微分方程式）
を作ればよいでしょうか？

a2=.;
For[n=1,n<2*10^1,n++,
If[n==1,s1=0];
If[n>1,s1=n/y];
y=3/2*a2^2+a2^3;
a2=2-s1;
f=5;
If[y-f<0,Print["n=",n]];
If[y-f<0,Print["y=",N[y]]];
If[y-f<0,Break[]];
];

もし、正しい式が得られましたら、ｆの値を
変更しても、正しい解が得られるはずです。
0902１３２人目の素数さん2007/10/03(水) 22:46:35
質問の本質はなんなのさ。
0903１３２人目の素数さん2007/10/04(木) 00:27:12
>902
お返事ありがとうございます。
＞質問の本質はなんなのさ。
誤解を招く書き方をして申し訳ございません。
本来の目的は、896に書いたプログラム（数値計算）
から導かれる答えを、数値計算ではなく、式を作って
解いて、同じ答えを導きだしたかったのです。しかし
896に書いたプログラムは、やたら複雑でわかりにくいので
901のとおりプログラムを単純にしました。

私の希望はあくまでも、数値計算で求められる解を
数式を解いて、得ることです。
901で答えをご教示頂きましたら、それをそのまま
896に私が置き換えるつもりです。

0904１３２人目の素数さん2007/10/04(木) 00:59:25
教えてください(>_<)
整数a、bを係数とする2次方程式x^2＋ax＋b＝0が有理数の解p/q(p、qは互いに素な整数)をもつならば、その解は必ず整数であることを証明せよ。
0905１３２人目の素数さん2007/10/04(木) 08:44:53
>>901
Mathematicaの問題じゃなくなったのでよそでやってくれ。
数学の初心者向けの板で馬鹿にされてくるといいよ。
0906１３２人目の素数さん2007/10/04(木) 20:46:16
>>905
お返事ありがとうございます。

>Mathematicaの問題じゃなくなったのでよそでやってくれ。
数％は、mathematicaの問題でもあるのでは？
>数学の初心者向けの板で馬鹿にされてくるといいよ。
まあ、そう言わずに、解っているのなら、せめて
答えだけでも教えて下さいよ！！直ぐに退場します。
この問題は、初心者じゃわからないです、、、
でも、かっこいいこと言って、本当は、答えが？では

0907１３２人目の素数さん2007/10/05(金) 09:57:50
ぷ。
0908１３２人目の素数さん2007/10/05(金) 11:09:51
答えは必ずあって、方程式には必ず解析解がある。
つまらなそうな世界だな。
0909１３２人目の素数さん2007/10/05(金) 13:00:10
Premier入ってホームユーザのアップデートきたんだけど、
パスワードの取得がやっかい。
Webに詳細な個人情報とパスワードが必要な理由ってのを打ち込んだらきたパスワードが４日間だけ有効って代物。

まともなパスワードほしければこれを見ろってんでみたらPDFファイル。
それ印刷して手書きで個人情報書き込んでパスワードがほしい理由を英語で説明して
アメリカへファックスし、approvalが降りるの待てときたもんだ。

何十万も払った上、事務手続きまで要求されるとは。
0910１３２人目の素数さん2007/10/05(金) 14:47:13
そんな面倒なこともあるの?
俺もPremierだけど、「あんたはPremierだからこっからダウンロードしてくれ。
ちなみにパスワードはこれ」とか言うメールが来て、別に面倒はなかったけど。
0911１３２人目の素数さん2007/10/05(金) 15:22:38
前にダウンロードで失敗したことがある。途中で切れると二回はできない。
お前はもうダウンロードしたからだめと表示される。途中で切れててもね。

連絡して再認可が下るのを待てばいいんだけど、また失敗すると面倒なんで
今回はあえてCD送付を選んだわけ。そしたらこちらは事務仕事が強制だったw


ところで起動して眺めたら、デバッガーがついたみたいだけどこれって前からあったっけ？
もし使いやすいデバッガーなら嬉しい驚きなんで、とりあえず簡単なプログラム
書いて試したけど、やはりマニュアル読まないとだめみたい。
もう使ってみた人いたら感想教えて。
0912１３２人目の素数さん2007/10/05(金) 20:31:51
>>908
>つまらなそうな世界だな。
皆様、超多忙にも関わらず、結局結論は出ず、ただ時間の無駄使いでした。
結論：お騒がせしました事をお詫び致します。
0913１３２人目の素数さん2007/10/05(金) 21:05:45
誤爆？
0914１３２人目の素数さん2007/10/09(火) 00:01:52
>>895
ありがとう。
0915あ2007/10/14(日) 00:52:10
P[x]という式を長軸a、短軸bの楕円部分で積分するにはどうしたらよいか教えてくれませんか？
0916１３２人目の素数さん2007/10/14(日) 06:36:15
>>915
例
Integrate[x^2 y^2 Boole[x^2/a^2 + y^2/b^2 <= 1],{x,-Infinity,Infinity},{y,-Infinity,Infinity},Assumptions->{a>0,b>0}]
0917あ2007/10/15(月) 01:04:03
>>916
ありがとうございます！実行したらうまくいきました！感謝です！！
｡･ﾟ･(ﾉД`)･ﾟ･｡
0918１３２人目の素数さん2007/10/15(月) 01:21:30
連投すんません(´Д｀ヽ)(/´Д｀)/
mathematicaで例えばT[x]とR[x]を畳み込み積分させるにはどうすればいいですか？
ググってもわかりまへんでした(#+_+)
0919１３２人目の素数さん2007/10/16(火) 22:15:13
ver.5とver.6って共存できる？
0920１３２人目の素数さん2007/10/16(火) 23:53:21
>919
問題なくできる。まあ同じ設定ファイルを読もうとするから書き方に気をつける
必要あるけど
0921１３２人目の素数さん2007/10/17(水) 16:03:19
>>920
ありがとう！
ところで、
> まあ同じ設定ファイルを読もうとするから書き方に気をつける
というのは、どういうことでしょうか？
6と5で書き方が異なるということですか？
0922１３２人目の素数さん2007/10/18(木) 01:12:35
>921
添付のライブラリの構成が大幅に変わってるので、5で入れてたNeedsとかをその
まま6に読ませるとWarningでまくりってこと。
0923１３２人目の素数さん2007/10/18(木) 14:55:37
mathematica6の日本語の体験版がないのはなぜ?
0924１３２人目の素数さん2007/10/18(木) 15:24:46
>>922
ご丁寧に有難うございました。気をつけます！
0925１３２人目の素数さん2007/10/19(金) 01:46:17
UnitStepを使って半径rの円内部にだけP[x]を掛けたいんだけど
θの範囲(0〜π)の指定ってどうしたら出来ますか？？

[現在の式]
T[x_] :=P[x]*(UnitStep[x+r*Cos[θ]]-UnitStep[x-r*Cos[θ]])
0926１３２人目の素数さん2007/10/20(土) 00:03:24
どなたか mathematica 6 のノートブックのファイル形式を
5.2 の形式に変換する方法をご存知の方おられませんか。

一度 6 で開くと 5.2 で開けなくなりました。
もちろん、関数まで自動で書き直してもらう事は期待してません。
（RandomIinteger[] から Random[Integer, x] に戻すとか）

私のやってる計算では、5.2 の方が３倍速いので、6 で書いた
ものを 5.2 のノートブックにしたいのであります。
0927１３２人目の素数さん2007/10/20(土) 00:50:37
>>926
やっぱりむやみやたらに６にするのは危険だな。
もう少し待っていよう・・・
0928１３２人目の素数さん2007/10/20(土) 01:26:08
拡張子maっていつのバージョンのファイルですあ？
09299262007/10/20(土) 01:46:14
>>926
>>927
すいません、 3倍じゃなくて1.3倍です。ひどいタイポでした・・・
しかし 6 はメモリの消費量もまたすごいです。

.m ファイルにするのはちょっとなあ。
0930１３２人目の素数さん2007/10/20(土) 18:23:25
見栄えを浴するためにメモリを食うのかな？　Ｖｉｓｔａ傾向というか。
0931１３２人目の素数さん2007/10/20(土) 19:33:22
英語版のmathematica使ってるんやけど
グラフの軸に
X軸距離
って名前つけて、その図をワードにコピペしたらX以外文字化けしました(ToT)
コレどーにかなりませんか？
0932１３２人目の素数さん2007/10/20(土) 19:39:14
>>931です
すんません、自己解決しました
Copy As⇒Bitmap
でコピペしたら文字化けせんかったです
0933１３２人目の素数さん2007/10/20(土) 21:56:35
数学オリンピックのリーダー格だった開成から
５大会連続で代表が出てないのは異常だな。
選考方法に問題がありそうだな。
まあ詳しく言わないけど。
0934１３２人目の素数さん2007/10/20(土) 22:36:15
Mathematica4です。ListPlotでプロットする点を○や△などの記号に変えることはできますか。あと、できれば各記号ごとに中間値で直線を引きたいです。
無理なのでしたら、ほかにどのような方法（座標指定して直接記号入力など）がいいと思いますか。
0935１３２人目の素数さん2007/10/20(土) 23:37:11
>934
マニュアルを読むのが良いと思うよ。MultipleListPlotとかPlotJoined->Trueとか
書いてあるから。
0936１３２人目の素数さん2007/10/20(土) 23:59:40
Mathematicaを初めて購入しようと検討中ですが
Win用とMac用、どちらが良いでしょうか。
因みにパソは両方あります。
0937１３２人目の素数さん2007/10/23(火) 00:23:28
メモリをたくさん使える方。
0938１３２人目の素数さん2007/10/24(水) 16:04:51
ParametricPlot3D[{(2 + 0.5Cos[u])Cos[v], (2 + 0.5Cos[u])Sin[v], 0.5Sin[u]}, {u, 0, 2Pi}, {v, 0, 2Pi}]
を使用し３Dの図形を書いたのですが、それをどうすれば回転させることができますでしょうか？よろしくお願いします。
0939１３２人目の素数さん2007/10/27(土) 14:49:46
マウスで図形をぐりぐり動かすと言う意味？
Ver6 なら特に設定なしでできるし
Ver5,4あたりなら
<<RealTime3D`
ParametricPlot3D[{(2 + 0.5Cos[u])Cos[v], (2 + 0.5Cos[u])Sin[v], 0.5Sin[u]}, {u, 0, 2Pi}, {v, 0, 2Pi}]
でできたような気がした
0940１３２人目の素数さん2007/11/01(木) 01:08:26
Mathematica5.2

・計算終了予想時間を計算中に出力することはできますか？
・カーネルの中身を次回起動時にも参照できるように残すことはできますか？
0941１３２人目の素数さん2007/11/04(日) 05:43:47
Ver6ってやっぱお布施だったのか？あんまレビューみたいなのないし
0942１３２人目の素数さん2007/11/10(土) 11:38:29
6っていったい何が一番変わったのですか？
0943１３２人目の素数さん2007/11/10(土) 13:22:31
メジャーバージョンが奇数から偶数に変わった
0944１３２人目の素数さん2007/11/11(日) 00:02:19
グラフはきれいになった。
インタラクティブ性がよくなったそうだけど、使いそうもないのでわからん。
デバッガーが入ったようで喜んだんだけど、使い方がわからない。ヘルブで検索してもWebが呼ばれるだけ、それらしいところ探してたら一巡しちまって結局分からず。
昔のノートブック読み込んでセーブするとver6の形式になりver5ではまともに読めなくなって、これは困った。情報センターのシステムはver5だったりするから。ver5形式でセーブってオプションは見つからなかった。
0945１３２人目の素数さん2007/11/11(日) 23:26:25
Iｆ文で、
If[条件,　式1,　式2]
とした場合、式1や式2を複数の式にするにはどうしたらいいんでしょうか？

0946１３２人目の素数さん2007/11/11(日) 23:59:00
それだけでは何がしたいのかわからんｗ
0947１３２人目の素数さん2007/11/12(月) 00:01:48
>>945
; で継げる. つまり

If[条件, 式 1a; 式1b; 式 1c, 式 2a; 式2b; 式 2c]

TRUE なら 式 1a; 式1b; 式 1c
FALSE なら 式 2a; 式2b; 式 2c
0948１３２人目の素数さん2007/11/12(月) 15:03:14
>>947
私も；を使ってみたのですが、そのようにすると、
式 1a、式1b、式 2a、式2bが評価されない答えが返ってきてしまいます。
aからc全て評価されるようにするにはどうしたらよいのでしょうか？
0949１３２人目の素数さん2007/11/12(月) 23:01:36
ひょうかはされているのじゃ
0950１３２人目の素数さん2007/11/13(火) 01:56:41
studentバージョン
ってどんな制限付くのですか？
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
にしろ
レス数が950を超えています。1000を超えると書き込みができなくなります。