〓 Mathematica 2 〓

1 名前：１３２人目の素数さん mailto:sage [04/03/06 09:51] 　　　　　　　　_....._{{ 〃 　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､ 　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､ .　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-< 　 　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）) 　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ. 　　　',　',　|Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　 ﾞ', 　　　 ヽ､', l:!Kノ}. 　 　 f:_.)iﾞi: ﾘ !　l　ル　 　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉ ﾞ ,ﾘ l ﾘ'´ 　 　　 .　　　　 ',|!!､ 　 　r‐┐ 　 ｀ ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣ 　 　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　 ＜　 Mathematica に関する話題はここに書いてね！ 　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　 　 　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿ 　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､ 　 　　/＼＼　　　　',　　 }　　 //｀ヽ 　★前スレ 〓Mathematica〓 science2.2ch.net/test/read.cgi/math/985023298/ 381 名前：１３２人目の素数さん [04/07/27 17:38] >>374 ContourPlotで行けるんでね？ 382 名前：１３２人目の素数さん mailto:sage [04/07/27 17:47] >>379 バーカ、市ね。くたばれ 383 名前：１３２人目の素数さん [04/07/28 03:25] TraditionalForm[ Cos[x/5-1] ] いじょ 384 名前：１３２人目の素数さん [04/07/28 03:27] \!\(TraditionalForm[Cos[x\^2 + x/5 - 1]]\) 385 名前：１３２人目の素数さん [04/07/29 23:30] ArcTan(163/217)でArcTanを求めるにはどのようにすればいいのでしょうか？ 最終的に" ° "で答えが欲しいのですが... 今日体験版を落として、触っているのですが分からなくて・・・ 386 名前：１３２人目の素数さん [04/07/29 23:35] 電卓で計算 387 名前：385 mailto:sage [04/07/29 23:41] 簡素化して書いているので実際はもう少し複雑で代入を繰り返さなければならないのでなんとか教えて頂けないでしょうか・・・ 388 名前：385 mailto:sage [04/07/29 23:56] すみません、間違ってました。 ArcTan(163/217)の答えを出したいのですが 答えも ArcTan(163/217)とでて" ° "で答えが欲しいのですが... 389 名前：１３２人目の素数さん [04/07/29 23:59] >>338 \!\(N[ArcTan[163\/217]\ 180\/π]\) 390 名前：385 mailto:sage [04/07/30 00:10] 本当に助かりました！どうもありがとうございましたm(_ _)m 391 名前：１３２人目の素数さん mailto:sage [04/07/31 16:43] おれはN[]よりも//Nの方を良く使うな 392 名前：１３２人目の素数さん [04/08/11 01:33] 461 393 名前：１３２人目の素数さん [04/08/12 01:54] mathematica student 4.1を使っていたのですが今年の3月で卒業になったので、 ５．０proにup gradeすることに決めました。 そこで聞きたいのですが 今現在パソコンを買い換えたのでmathematica　student版が インストールされていません こういう場合でもアップグレードできるのでしょうか？ 394 名前：１３２人目の素数さん [04/08/12 17:14] >>393 大丈夫。 とりあえずここら辺でも見ておけ。 ttp://www.varsitywave.co.jp/products/mathematica/html/math-upgrade.html 395 名前：１３２人目の素数さん [04/08/19 17:36] x軸として x={x1,x2,x3,...} y軸として y={y1,y2,y3,....} という感じにデータセットがあるんですが（要素数は同じ）、 これを test = { {x1,y1}, {x2,y2},},....} といった感じのx,yを対応させたリストに するにはどうしたらいいでしょうか。すべてただの実数です。 どなたかご教示願います・・・。 396 名前：１３２人目の素数さん [04/08/19 22:48] 254 397 名前：１３２人目の素数さん [04/08/20 12:10] >395 それはできますね マニュアルかヘルプを読みましょうよ 398 名前：１３２人目の素数さん mailto:sage [04/08/20 12:48] transpose 399 名前：１３２人目の素数さん mailto:sage [04/08/22 03:04] >>397,398 ありがとうございました。 簡単にできました。申し訳ないです。 400 名前：１３２人目の素数さん [04/08/28 07:32] 400 401 名前：１３２人目の素数さん [04/08/29 11:55] Ver5を持っている人へ ＧＵＩＫＩＴが出ています。 公式ホームページから ライセンスナンバーだけでダウンロードできます。 インストールする場所と ヘルプファイルの更新をお忘れなく。 402 名前：１３２人目の素数さん [04/08/31 19:12] NxNのでっかい対称行列の固有値を求めたいのですが、 Cのプログラムでmatrix.txtと 1.0 1.1 1.2 ... ←N個 1.1 2.0 1.3 ... ←N個 ... ↑N個 の形式でN行で出力させてから ReadList["matrix.txt",Number,RecordLists->True] Eigenvalues[%] とMathematicaに読み込ませて固有値を計算させようとしています。 Nが200位になるとファイルの読み込みで非常に時間がかかってしまい、 必要としているN=1000程度だと2日ほど放置しましたが読み込みが完了しません。 固有値の計算で時間がかかるのならわかる気がするんですが、 ファイルの読み込みの段階で時間がかかるのが納得いきません。 これは仕様で仕方ないのでしょうか？ それとも私の読み込ませ方に問題があるのでしょうか？ 計算機はXeonの2GDualメモリ1G（だったかな…） VineLinuxでMathematicaのVersionは4.2.1.0です 突然研究で固有値を求めることになったMathematica超初心者ですが、 ご教授いただければと思います。゜゜(´□｀｡)°゜。 403 名前：FeaturesOfTheGod ◆UdoWOLrsDM [04/08/31 19:15] Re:>402 私はよく分からないけど、読み込みを何回も行っているのではないだろうか？ もしそうでないとしても、固有値の計算は時間がかかる。(少なくとも私の知っているアルゴリズムでは。) 404 名前：１３２人目の素数さん mailto:sage [04/08/31 22:01] >>402 読み込みにではなく、表示に時間がかかっている。 どうしても1000*1000の行列を眺めたいのでなければ   ReadList["matrix.txt",Number,RecordLists->True];   Eigenvalues[%] とするといいだろう。 405 名前：１３２人目の素数さん [04/08/31 23:09] デカルト座標から極座標に変換するように、下記のようにしたんですが この角度のパラメータをラジアン表記でなく度数表記で一発出力させる ためには、どういじればいいのでしょうか。 << Calculus`VectorAnalysis` SetCoordinates[Cartesian[x, y, z]]; x = 1.; y = 1.; z = 1.; CoordinatesFromCartesian[{x, y, z}, Spherical] 406 名前：１３２人目の素数さん mailto:sage [04/08/31 23:35] >>405 CoordinatesFromCartesian[{x, y, z}, Spherical] / Degree HoldForm[Degree] 407 名前：１３２人目の素数さん [04/09/01 00:21] 405です >>406 ありがとうございます。実行してみました。 返数のうち最初の一つは長さなので、角度である後二つにDegreeを かけたいのですが、方法はあるのでしょうか。 すいません、初歩的なことで。 408 名前：１３２人目の素数さん mailto:sage [04/09/01 01:17] >>407 すまん。こっちこそ初歩的に間違えた。 CoordinatesFromCartesian[{x, y, z}, Spherical]  {1, HoldForm[Degree] / Degree, HoldForm[Degree] / Degree} # ラジアンに戻す時は # x = Coord...<snip> # ReleaseHold[x] 409 名前：１３２人目の素数さん [04/09/01 01:35] 405 >>408 ありがとうございます。できました。 勉強させていただきます。 410 名前：402 mailto:sage [04/09/01 10:52] >>404 ありがとうございました。それで解決いたしました。 そんなことだっとは盲点でした……。 次のような記事も発見しました。 私の場合ここまでする必要はありませんでしたが、参考までに。 『ReadList: Is it really slow?』 ttp://forums.wolfram.com/mathgroup/archive/2000/Jul/msg00264.html 411 名前：４１１ [04/09/02 19:38] Animate[]やFor[]で生成される複数の３Dグラフィックスを 同じ座標内に全て表示するようにしたいんですけど、どうすればいいんでしょうか？ 412 名前：１３２人目の素数さん mailto:sage [04/09/04 15:39] >>411 SelectionAnimate か、 表示されたものをマウスで畳んで、 右クリックしてAnimateみたいなメニューを選ぶ。 413 名前：１３２人目の素数さん [04/09/04 18:02] 馬鹿な質問ですが、 ギリシャ文字はマウスを使ってボックスから選択しないと使えないのでしょうか。 キーボードから入力するショートカットがあれば教えてください。 釣りではないのでよろしく。 414 名前：１３２人目の素数さん mailto:sage [04/09/04 21:48] >>413 ぜってー釣りだな 415 名前：１３２人目の素数さん mailto:sage [04/09/04 22:10] >>413 \[CapitalAlpha] \[Alpha] <Esc>Alpha<Esc> <Esc>alpha<Esc> 416 名前：FeaturesOfTheGod ◆UdoWOLrsDM [04/09/04 22:17] Re:>415 何故もっと簡単な方法を教えてやらぬのだ？(私もだ。) 417 名前：４１１ [04/09/05 10:24] >>412 それでShowをつかって複数のグラフィックスを同一座標内に表示したように グラフィックスを重ねて表示する事ができるんでしょうか？ もうちょっと具体的に教えてくれるとうれしいんですけど・・ 418 名前：１３２人目の素数さん [04/09/05 15:33] 絵を描きたいんですが、何かいい物ありませんか？ 例えば、ドラえもんとかドラゴンボールとか… 419 名前：FeaturesOfTheGod ◆UdoWOLrsDM [04/09/05 15:36] Re:>418 図形をたくさん組み合わせる。楕円、長方形、線分、その他いろいろある。 420 名前：１３２人目の素数さん mailto:sage [04/09/05 16:25] >>412 アニメーションがしたいわけではなくて、単に重ねたいのなら 最初からShowか元のPlotなどで重ねろ。 Animatieの出力を利用するのはNotebook操作系の関数を使えば できないことはないが、二階から目薬。 421 名前：４１１ [04/09/05 17:26] >>420 二回から目薬とは？？ どのようなNotebook操作系の関数を利用すればできますか？ 422 名前：１３２人目の素数さん [04/09/05 17:29] 有名な諺 423 名前：１３２人目の素数さん [04/09/10 15:33:35] 224 424 名前：４１１ [04/09/12 19:28:00] Mathematicaで4次元グラフィックスが比較的簡単に書けるようなパッケージありますか？ 425 名前：１３２人目の素数さん [04/09/14 20:41:52] {-X*a*(1+X^2)(d+d*X^2+e)}/(b-f+b*X^2)+X*(c-f)+X^3*c+ {(b-f+b*X^2)*e}/{X*(d+d*X^2+e)}==0,X 426 名前：１３２人目の素数さん [04/09/14 20:44:02] 誰か解いてください。お願いします。 427 名前：１３２人目の素数さん mailto:sage [04/09/14 22:22:47] Mathematicaで解いたら？ 428 名前：１３２人目の素数さん [04/09/14 22:58:42] 今、Mathematicaがないから解いてと頼んでいます。 429 名前：１３２人目の素数さん [04/09/19 17:58:16] 552 430 名前：１３２人目の素数さん mailto:sage [04/09/21 16:57:02] ノートブックを編集中のキーバインディングを変える方法 というか割り当ての方法をどなたかご存じないですか？ 例えば コントロール＋ｋ を行末までの文字の削除に当てるとか。 431 名前：１３２人目の素数さん [04/09/21 19:14:09] mathematicaはちょっとした計算に向くけど、 なんか大きなことやろうとすれば、matlab 使うなー。 mathematicaは言語の外見が汚なすぎ。 てかセンス悪いよな、Wolframは 432 名前：１３２人目の素数さん [04/09/21 19:15:47] mathematica での開発は効率が悪すぎ 433 名前：FeaturesOfTheGod ◆UdoWOLrsDM [04/09/21 20:33:05] Re:>432 そもそも、mathematicaは研究のためにある。 434 名前：１３２人目の素数さん mailto:sage [04/09/21 21:22:14] そもそもFeaturesOfTheGod ◆UdoWOLrsDM は馬鹿にされるためにある。 435 名前：FeaturesOfTheGod ◆UdoWOLrsDM [04/09/21 21:28:32] Re:>434 馬鹿にされるのはお前だけで十分。 436 名前：１３２人目の素数さん [04/09/25 04:05:10] >>431 >mathematicaは言語の外見が汚なすぎ。 >てかセンス悪いよな、Wolframは そこはかなり同意する 437 名前：１３２人目の素数さん mailto:sage [04/09/25 09:52:28] コマンドラインで手軽に計算できるというのと、制御構文使って ある程度のプログラムを組めるというのを両立させようとするから ああいうことになるんだよな。 438 名前：１３２人目の素数さん [04/09/26 16:45:39] >>437 そうですね。mathematicaは開発者が多いからグラフィックスやGUIだけが 他のソフトより良いだけであって、Wolfram本人もキモいし、言語もキモいし、 やたらプログラミングに時間がかかるだけですね。 もうmathematicaは逝ってよしだと思いました 439 名前：１３２人目の素数さん [04/09/26 16:51:10] mathematica使いはWolframのように若ハゲしそうですね。 がんばってください。キモっ! 440 名前：１３２人目の素数さん [04/09/26 16:53:03] ところで、「制御構文使って・・・」の言語はいくつも存在しますよ。 Rubyとかは特にね。 matlab言語も（・∀・）ｲｲ!!よ 441 名前：１３２人目の素数さん mailto:sage [04/09/27 14:21:50] MATHEMATICA撲滅しよーぜ 442 名前：１３２人目の素数さん [04/09/27 19:07:27] Mathematica氏ね。 443 名前：１３２人目の素数さん [04/10/03 18:13:37] 987 444 名前：１３２人目の素数さん [04/10/04 15:09:33] mathmaticaで例えばSinを変数名として使っちゃた場合とか どうやったら元の関数として機能してくれるの？ 445 名前：FeaturesOfTheGod ◆UdoWOLrsDM [04/10/04 19:37:46] Re:>444 Sinをどうやって変数名に使うのだ？protected属性を外した？ 446 名前：FeaturesOfTheGod ◆UdoWOLrsDM [04/10/04 19:38:30] Re:>444 protectedじゃなくて、Protectedね。とりあえず、Sin[Pi]などで正常に値が出るなら大丈夫だと思うけど。 447 名前：１３２人目の素数さん mailto:sage [04/10/04 21:15:25] >>444 評価の仕組み上無理だったと思います。 どうしてもというなら、Symbol["Sin"]でもどうぞ。 Func[Sin_] := {Sin, Symbol["Sin"][Sin]} Func[3.0] # /* Example.c */ # int x = 3 # f(int x) { return x; } # /* さて、f(3)の返り値は？ */ >>446 そういう場合にはちゃんとDefinitionやFullDefinitionを使いましょう。 # 前身時代からなんで知識があるわけでもないのに考えずに適当なこと # ばかりこたえるかな。 448 名前：１３２人目の素数さん [04/10/07 00:22:39] グーグルで検索していてここにたどり着きました。 >>405-408 のプログラムを利用させてもらっています。 これに改良を加え座標を方程式で出してその値を極座標に変換できるように しようと試みたのですがうまくいきません。 次のようにしてみました。 << Calculus`VectorAnalysis` Solve[{x + y + z == 1, 　 x - y + z == 3, x + y - z == 5}, {x, y, z}] CoordinatesFromCartesian[{x, y, z} /. %, Spherical] {1, HoldForm[Degree]/Degree, HoldForm[Degree]/Degree} 返数は CoordinatesFromCartesian[{{4, -1, -2}}, Spherical] （以下同じような感じで二つ返ってきてしまいます。） どうもただ式のx, y, zに方程式の解が代入されただけの感じのものです。 どうすれば成功するのでしょうか。アドバイスもらえないでしょうか。 449 名前：１３２人目の素数さん mailto:sage [04/10/07 04:21:53] どうも「評価」というのがわからない。 Table[If[PrimeQ[i], i, ""], {i, 2, 100}] を実行すると、 {2, 3, , 5, , 7, , , } が返るが、 list = {2}; Do[ If[PrimeQ[i], Append[list, i]], {i, 3, 100}]; list とやっても {2} となるだけ。 なんでなのでしょうか。 450 名前：１３２人目の素数さん mailto:sage [04/10/07 12:53:46] >>448   CoordinatesFromCartesian[{{4, -1, -2}}, Spherical] でなくて   CoordinatesFromCartesian[{4, -1, 2}, Spherical] {1,  でないといけない。Solveが複数答えを返すこともあるので ans = Solve[{x^2 + y + z == 1, x^2 - y + z == 3, x^2 + y - z == 5}, {x, y,        z}] Map[CoordinatesFromCartesian[{x, y, z} /. #, Spherical] {1,          HoldForm[Degree]/(Pi / 180), HoldForm[Degree]/(Pi / 180)} &, ans] など。 >>449 list = {2} Append[list, 3] list と入力すると何がおこるでしょう？ 451 名前：１３２人目の素数さん mailto:sage [04/10/07 14:08:09] >>450 あうう… 「評価」じゃなかったのですね。 思いこみであちこちにEvaluateつけて考えあぐねておりました Thanks a lot! 452 名前：１３２人目の素数さん mailto:sage [04/10/10 21:32:56] 見当違いのことをしているのかもしれませんが、教えてください。 Mathematicaでの数値計算の精度は機械精度以上には上げることが出来ないでしょうか。 $MachinePrecision Out[1]=15.9546 a1=SetPrecision[Sqrt[2-0.1],100]; NumberForm[a1,DigitBlock->10,NumberSeparator-> " " ] NumberForm= 1.3784048752 0902210683 0047960102 0924746990 2038574218 7500000000 ＼ 0000000000 0000000000 0000000000 000000000 a2=Sqrt[SetPrecision[2.-0.1,100]]; NumberForm[a2,DigitBlock->10,NumberSeparator-> " " ] //NumberForm= 1.3784048752 0902214457 7973787659 2686126529 8801258806 7606899563 ＼ 5042810459 9163905598 4762095066 7510687484 a1とa2を比較すると、機械精度が保証する桁以下で両者が異なっています。 どちらが正なのでしょうか。SetPrecisionの使い方がまずいのでしょうか。 453 名前：１３２人目の素数さん mailto:sage [04/10/10 22:39:04] >>452 Sqrt[2 - 0.1]の結果の数値を精度100桁にするというのと、 (2. - 0.1)の結果の数値を精度100桁としてからSqrtをとるのとは 全く別のことですが。 もちろんどちらもN[Sqrt[2 - 1/10], 100]とは異なりますね。 454 名前：１３２人目の素数さん mailto:sage [04/10/16 06:43:32] 709 455 名前：１３２人目の素数さん [04/10/16 12:12:51] 455 456 名前：１３２人目の素数さん mailto:sage [04/10/16 21:27:04] >>452 ひょっとして　０．１と１／１０が同じと思っていませんか？ 457 名前：１３２人目の素数さん mailto:sage [04/10/17 17:04:18] >>456 一緒と考えておりました... これで納得できました。ありがとうございました。 458 名前：１３２人目の素数さん mailto:sage [04/10/19 02:55:15] マスマティカはお腹が一杯になりまつか？ 459 名前：１３２人目の素数さん mailto:sage [04/10/22 14:47:57] ある分布に従う乱数を発生させて標本平均値を求めたいのですが、 発生させた乱数のリストから任意の数をランダム抽出するにはどうしたらいいでしょうか？ 460 名前：１３２人目の素数さん mailto:sage [04/10/22 18:01:07] >>459 乱数から「ランダム抽出」って意味が無いだろ。 Table[YourRandomFunction[], {i, 10}] とでもすればいいんじゃないか？ 461 名前：Not, used with Select [04/10/25 01:34:56] Select の条件部分に Not を入れてはいけないのでしょうか．うまく行かない理由が分かりません． たとえば Select[Range[10], PrimeQ] Select[Range[10], ! PrimeQ] にたいしてはアウトプットが{2,3,5,7}と{}となり， ２つ目が意図したものになりません． Select[{b, a, d, c, e}, (# == a) &] Select[{b, a, d, c, e}, (! (# == a)) &] にはアウトプットが{a}と{}となり， これも Not (!) を使った２つ目が意図とはちがいます． ちなみに Select[Range[10], (# == 4) &] Select[Range[10], (! (# == 4)) &] だと意図した通りに {4}と{1,2,3,5,6,7,8,9,10}となります．どうして？ ちなみに問題回避法としては　!MemberQ というのを使いました．もっとマシな方法あれば知りたし． 462 名前：１３２人目の素数さん mailto:sage [04/10/25 01:57:45] >>461 ! PrimeQはNot[PrimeQ]ですね。 Not[PrimeQ][x]とNot[PrimeQ[x]]は違いますね。 Not[PrimeQ[#]]& か Composition[Not, PrimeQ]を使ってください。 !MemberQなんてものをいったいどう使ったのか興味津津です。 463 名前：Not, used with Select (unresolved) [04/10/26 23:51:59] >>462 >!�PrimeQはNot[PrimeQ]ですね。 0-1逆の特性関数を作るわけではないようですね． 普通の関数が否定できないのは分かるが，述語関数は否定できてもまずくなさそうに思えますが． >Not[PrimeQ][x]とNot[PrimeQ[x]]は違いますね。 前者は意味をなさないということなんですかね． [ところで同じ意味をなさない物でも， ! PrimeQ[3]だと何も出ないが(! PrimeQ[3])だと False と出るのは， 前者は括弧の位置が曖昧だが，後者は(! PrimeQ)[3]ではないと見なされる ということなのだろうか．一番外側に括弧がつくだけで差が出るのは妙だが．] >Not[PrimeQ[#]]&�か�Composition[Not,�PrimeQ]を使ってください。 ありがとうございました．続く． 464 名前：Not, used with Select (unresolved) [04/10/26 23:52:33] ところが本来の目標の方なんですが， Select[{b, a, d, c, e}, ! (# == a) &] Select[{b, a, d, c, e}, Not[# == a] &] ではどうしても {} となります．分からない〜！！！ 等しくなる可能性があるから a=b, a=c などは False と言えないということでしょうか． 値はどうであれ(というか値は取らない)， 「シンボルとして a とは異なるものを取り出す」 ことを表現できればいいのですが． どこかに{a,b,c,d,e}は(異なる要素からなる)集合だといっておけばいいのだろうか． そんなコマンドあるのかな． >!MemberQなんてものをいったいどう使ったのか興味津津です。 Select[{b, a, d, c, e}, ! MemberQ[{c}, #] &]でも Select[{b, a, d, c, e}, Not[MemberQ[{c}, #]] &]でも {b, a, d, e}を出しました． 465 名前：１３２人目の素数さん mailto:sage [04/10/27 00:11:51] >>463 > ! PrimeQ[3]だと何も出ないが(! PrimeQ[3])だと False と出るのは， !が「行の先頭」に来た時はshellでのコマンド実行を意味する。 > 等しくなる可能性があるから a=b, a=c などは False と言えないということでしょうか． そのとおり。===を使えばいい。 > Select[{b, a, d, c, e}, ! MemberQ[{c}, #] &]でも それはあたりまえだね。 偶然、Not[MemberQ[ ]]& という形で使っていたわけだ。 466 名前：sage [04/10/27 12:04:06] From Not, used with Select (unresolved) >>465 行けそうです．ありがとうございました． ===は初耳でした(Mathematica 使い始めてまだ2週間くらいなもので)ので， その周辺のコマンドを勉強しておきます． 467 名前：１３２人目の素数さん [04/11/02 00:22:41] 323 468 名前：１３２人目の素数さん [04/11/04 10:06:12] z=-x^2+y+1 |x|,|y|<=3 において z>=0 の部分の地形を等高線で表せ。 という課題が出まして、 ContourPlot[-x^2 + y + 1, {x, -3, 3}, {y, -3, 3}, Contours -> 10, PlotRange -> {0, 4} といった具合にやらせてみたんですが、これだと真っ黒で見た目が非常に悪いんです。 色を使って高さを表すにはどうしたらよろしいでしょうか？ どなたかアドバイスお願いします。 469 名前：１３２人目の素数さん mailto:sage [04/11/04 19:48:54] > 468 color[i_] := Which[Mod[i, 1] < 0.25, CMYKColor[0, 1 - 4 Mod[i, 1], 1, 0], Mod[i, 1] < 0.5, CMYKColor[4 (Mod[i, 1] - 0.25), 0, 1, 0], Mod[i, 1] < 0.75, CMYKColor[1, 4 (Mod[i, 1] - 0.5), 1 - 4 (Mod[i, 1] - 0.5), 0], Mod[i, 1] ≥ 0.75, CMYKColor[1, 1, 0, 4 (Mod[i, 1] - 0.75)]] ContourPlot[-x^2 + y + 1, {x, -3, 3}, {y, -3, 3}, Contours ->10, ColorFunction -> (color[1 - 0.66*(#1 + 0.5)]&), PlotRange -> {0, 4}] で、どう？ 470 名前：１３２人目の素数さん [04/11/04 20:10:57] し 　　っ 　　　　か 　　　　　　し 　　　　　　　　、 　　　　　　　　　　ほ 　　　　　　　　　　　　ん 　　　　　　　　　　　　　　と 　　　　　　　　　　　　　　　　に 　　　　　　　　　　　　　　　　　　キ 　　　　　　　　　　　　　　　　　　　　タ 　　　　　　　　　　　　　　　　　　　　　　ネ 　　　　　　　　　　　　　　　　　　　　　　　　ー 　　　　　　　　　　　　　　　　　　　　　　　　　　言 　　　　　　　　　　　　　　　　　　　　　　　　　　　　語 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　だ 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　な 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　こ 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　り 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　ゃ 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　w 471 名前：ChaosicSoul ◆/yaJbLAHGw [04/11/04 20:27:15] Mathematicaの関数名は、ごく一部を除いて省略が無い。 どうにかならぬか？ 472 名前：パンチョ伊藤 [04/11/04 20:37:50] お前久しぶりやな >>471 > Mathematicaの関数名は、ごく一部を除いて省略が無い。 > どうにかならぬか？ マジレスしとくか。 InputAliases オプションとかで短くしとけ。 escape にキーバインドされてるコマンドを \ とかに新規でキーバインド登録しとくと 入力が早くなる。 もし、本気で楽にテクニカルコンピューティングをしたいのなら matlab に乗りかえることもよかろうな。 あれはいい。仕事が mathematica の(1/5) t で終わるからな。 mathematicaではマジでプログラムする気にならん。 故に自分の関数も増えていかん；だから効率も悪い。お前もわかるだろ？ 473 名前：１３２人目の素数さん [04/11/04 21:55:01] どこまでもウザイ奴だな＞King 474 名前：１３２人目の素数さん [04/11/04 22:18:14] お前もウザイな>>473 475 名前：１３２人目の素数さん [04/11/04 22:33:55] NIntegrate[{NIntegrate[ｆ[x], {x,0,∞}]}g(x),(x,0,∞)] (f(x)をxで積分したものにg(x)をかけて、さらにそれをxで積分したもの) こういった二重積分はmathematicaでは出来ないのでしょうか？ 狙った値がでなくて非常に困っております・・ アドバイスよろしくお願いします。 476 名前：１３２人目の素数さん mailto:sage [04/11/04 22:35:35] >>468 >>469のほうがきれいではあるが、簡単なものも。 ContourPlot[-x^2 + y + 1, {x, -3, 3}, {y, -3, 3}, Contours -> 10,    PlotRange -> {0, 4}, ColorFunction -> (Hue[(0.1 + #)/1.1] &)] >>471 補間のキーボードショートカットを使え。 おまえ前も同じこと書いてただろ。学習しないやつだな… 今は使ってないとか言って逃げたんだったか。 荒しはMathematicaが使えなくてトラウマになってるのかね。 477 名前：１３２人目の素数さん mailto:sage [04/11/04 22:39:47] >>475 > (f(x)をxで積分したものにg(x)をかけて、さらにそれをxで積分したもの) それ二重積分じゃなくて、別々に積分してただ掛ければいいでしょ。 > こういった > 狙った値 何が聞きたいのかわからない。質問は具体的に。 478 名前：475 mailto:sage [04/11/05 00:49:26] orz 問題設定は、>>475の数値積分を行うと確率0.95という値が出てくるというものです。 そこで実際に計算を行ったところ目的の0.95が出てこなくてなにかおかしいのかと思いました・・ 479 名前：１３２人目の素数さん [04/11/05 09:30:32] >>469 さん >>476 さん 468です。ありがとうございました。 実は,この問題さらに続きがありまして。 z=-x^2+y+1 とする。 (1)|x|,|y|<=3 において z>=0 の部分の地形を等高線で表せ。 (2)地図上に二点 A(-1,0) B(2,3) を最短距離で結ぶ道がある。 この道を実際にたどるときの道のりはどれだけか。 また,この道の最高地点の高さ(海抜)を示せ。 (3)地点 (1/2,2/3) における接平面が,x軸,y軸,z軸と形成する三角錐の体積 はいくらか。 といった問題なのですが。 正直取っ掛かりすら見つけられないほど難儀しております。 どなたかアドバイスお願いします。 教えて君ですみません。 480 名前：マンコ [04/11/05 09:40:04] つか誰だよおまえ。コテハンも付けれんで調子にのんなよ 481 名前：マンコ [04/11/05 09:41:21] 「エッ、俺？」って思った 479。 心配すんな。>>480 はおまえに向けたわけじゃない >>476 だよ。

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef