〓 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/
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これでどう？
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:22u2=[ [ { 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:210975１３２人目の素数さん
2007/12/08(土) 01:31:22参考に。
http://support.wolfram.com/mathematica/graphics/2d/twoaxisgraph.ja.html
0976１３２人目の素数さん
2007/12/08(土) 04:03:07aa={-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符号の切り替わる場所のリスト、でいいんだよね？
短くなったけどどうだろうこれ
cc = Join @@ Position[Sign /@ Times @@@ Transpose[{Drop[aa, -1], Drop[aa, 1]}], -1]
0978１３２人目の素数さん
2007/12/08(土) 09:46:20x^(2/3) (x + (x^2) x^(7/3))
に
x=k y^3 (aは任意の実数)
を代入して全ての括弧を取り払いxの関数に展開したいのですが、どうすればいいでしょうか？
0979１３２人目の素数さん
2007/12/08(土) 09:49:370980１３２人目の素数さん
2007/12/08(土) 10:17:46まず、
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ありがとうございます！
まだまだ素人なのでちょっと見た感じでは文法が理解できないですが；
週明けに早速試してみようと思います。
0982１３２人目の素数さん
2007/12/08(土) 15:40:000が入っている場合があるのか、あったらどうするのかがちょっと気になる
0983１３２人目の素数さん
2007/12/08(土) 16:10:36だめみたいです。
k^5 y^15 + k y^3 (k y^3)^(2/3)
↑
これが残ってしまって。
次の目的として最高次数を取得したいので、正確に評価できるよう
この部分も開いてしまいたいのですが…。
0984１３２人目の素数さん
2007/12/08(土) 16:37:200985１３２人目の素数さん
2007/12/08(土) 18:14:48えっと、>>976ではわかりやすいよう簡単に書きましたが、
実際のaaにあたる部分は特殊関数を用いた関数をサンプリングしたものなので、
ぴったり0になることは考えていないです。
しかもその関数はただの特殊関数ではなく、色々と複雑な処理を行っている
（さらにパラメータふたつを変えながらTableで回す）ので、
根の近似値を初期値としてFindRootで解くということが（現状の私の技術では）できないのです。
この方法で解いた値を初期値として解く事は勿論可能でしょうが、現状以上に膨大な時間が
かかることになる＆そこまでするほど厳密な値を必要としていないので。
0986980
2007/12/08(土) 20:22:12Collect[expr, y,Simplify[#,Assumptions -> y>0] &]
だと
k^(5/3) y^5 + k^5 y^15
となるね。
Assumptionsでyが実数としてもy<0でもだめみたいだね。
どういうことなのかいまひとつよくわからない、、、。
0987980
2007/12/08(土) 20:30:04Simplify[expr, Assumptions -> y > 0]
だけでそうなる。Collectは必要ない、、。
0988１３２人目の素数さん
2007/12/08(土) 20:42:03ありがとうございます。
変数は当面は正で問題ないのでそれでやってみます。
バリバリ数式変形してから結果を入れたいときには鬱陶しいけど
数値入力する時のためにうかつに開けないようにしてるのかな…？
0989１３２人目の素数さん
2007/12/08(土) 23:17:27遅レスだが
http://support.wolfram.co.jp/applicationpacks/general/apppacksinv6.html
から引用「バージョン6.0より前のアプリケーションパッケージには
Mathematica 6.0との互換性はないとお考えください．お試しいただくことは
できますが，多くの機能が以前のようには動作しなくなっています．」
wolfram researchは何を考えてるのかね。
0990１３２人目の素数さん
2007/12/08(土) 23:24:530991１３２人目の素数さん
2007/12/09(日) 00:12:49えらく古そうだけど何の話？
0992１３２人目の素数さん
2007/12/09(日) 09:07:500993１３２人目の素数さん
2007/12/09(日) 16:37:210994１３２人目の素数さん
2007/12/09(日) 20:42:300995１３２人目の素数さん
2007/12/09(日) 21:00:13Mathematica Ver.6以下はクソ 4
にしろ
0996１３２人目の素数さん
2007/12/10(月) 16:10:00http://science6.2ch.net/test/read.cgi/math/1197270001/
0997１３２人目の素数さん
2007/12/10(月) 16:37:200998１３２人目の素数さん
2007/12/11(火) 05:36:230999１３２人目の素数さん
2007/12/11(火) 05:46:381000１３２人目の素数さん
2007/12/11(火) 05:47:131001１００１
Over 1000Threadもう書けないので、新しいスレッドを立ててくださいです。。。
レス数が1000を超えています。これ以上書き込みはできません。