〓Mathematica〓

1 名前：１３２人目の素数さん [2001/03/20(火) 02:34] ないので立てました。 Mathematicaについて語るスレです。 大学で初めて出会った時は感動したなぁ。(NeXT版…) Plot3Dでよく遊んだものです。 四次方程式の一般解を出力して、壁に貼ったりしてました（ｗ 445 名前：１３２人目の素数さん mailto:??? [03/06/02 18:40] 再インストールしる! 446 名前：１３２人目の素数さん [03/06/02 20:41] Mathematicaって使いにくくない？ライブラリとか凄いとは思うけど 447 名前：mathmania ◆uvIGneQQBs [03/06/02 20:46] コマンドが長いのは確かに使いづらい。 もっとも、私は最近Mathematicaをやっていない。 C言語(ForTranも少し)を使ってる。 理由はというと、、、Mathematicaが使えなくなったから。 まぁ、使いやすさの面でいえば、大きなプログラムには、C言語の方がいいだろう。 448 名前：444 mailto:sage [03/06/02 23:50] >>445 再インストールしても同じエラーが出るんです。 設定ファイルが残ってるんじゃないかと疑ってるのですが その場所が分かりません。 449 名前：１３２人目の素数さん mailto:sage [03/06/03 22:42] >>446 原始的なプログラミング言語と不統一な仕様の mapleよりは使いやすいと思う。 言語としてつかえてなんぼだから 素人にはお薦めできない、って感じはあるかもね。 >>447 補完機能をつかうべし。 450 名前：１３２人目の素数さん [03/06/07 17:45] 初心者にお勧めのMathemticaの 本は何でしょうか？ 451 名前：１３２人目の素数さん mailto:??? [03/06/07 21:09] >448 OSを再インストールしる! それでだいじょうぶだ! 安心しろ >450 白石さんのがいいぞ! 安心して買え 452 名前：１３２人目の素数さん [03/06/07 21:25] あなたが探してるのってこれだよね？二日間無料だしね♪ alink3.uic.to/user/angeler.html 453 名前：１３２人目の素数さん [03/06/07 22:15] mathmaniaスレかと思った 　　　　　　　　　　　　と言ってみる 454 名前：１３２人目の素数さん mailto:sage [03/06/08 00:02] >>450     Mathematicaプログラミング J.ゲイロード･N.カーミン･R.ウエリン を読んで、その後で     プログラミングMATHEMATICA. Ｒ・メーダー 著、時田 節 を読むとよくわかると思われ。 特にメーダーの本は(全くの初心者以外には)一読の価値あり。 非常に下手な翻訳だけど。 455 名前：１３２人目の素数さん [03/06/08 00:09] >>248 Windows Me はやめた方がいいですよ。 2KかXPにしましょう。 （マジレス） 456 名前：１３２人目の素数さん [03/06/08 02:09] Mathematicaって有効数字が６桁の表示で、例えば 1.23456789とか1234.56789って直接入力すると 1.23456や1234.56と表示されます。 ＃　表示だけがそうなるだけで内部的には入力した数のままです これを以降の入力からすべて指定した有効数字で表示させたいのですが どうやったら良いのでしょうか？NumberFormではそれで指定した場合だけ 表示されるのでこれを以降の入力すべてで有効にしたいのでよろしく お願いいたします。 457 名前：１３２人目の素数さん [03/06/08 10:38] >>451 >>454 どうもありがとう！ 458 名前：１３２人目の素数さん [03/06/08 10:49] >>456 ちょっと邪魔くさいけど、関数を定義すればどう？たとえば 　F[x_]:=N[x, 30]; 　F[E^3] とするとe^3が30桁表示されるはず。 459 名前：mathmania ◆uvIGneQQBs [03/06/08 13:02] 456の方法でも有効数字は確かに表示される。 だが、SetPrecision[E^3,Infinity]か、SetAccuracy[E^3,Infinity]はどうだろう？ まぁ、分数表示になるが、その数に2^52(?)を掛ければ小数点表示がわかるだろう。 あぁ、それと、私はMathematica復活したわけではないので。 460 名前：444 mailto:sage [03/06/08 21:28] >>451 結局OS再インストールしましたよ。。。 したくなかったけど。 >>455 大学のPCだから自由にできません。 経費で2000買ってくれなかった。 461 名前：１３２人目の素数さん [03/06/09 17:46] >>458-459 Mathematicaでの以降の入力すべてで指定した有効数字で表示させたいということ でしょうから、その方法ではその関数内の引数にだけ有効であって意味がありません。 何かシステムの設定でできそうな気もしますがそれらしいものが見当たりませんね。 どうやるんだろうか？ 462 名前：１３２人目の素数さん [03/06/09 18:00] >>461 SetPrecisionは計算精度を変更させるだけであってその表示は依然として 有効数字６桁です。 463 名前：１３２人目の素数さん [03/06/10 19:55] >>456 編集→環境設定→書式設定→式の書式→表示設定→PrintPrecision 464 名前：１３２人目の素数さん [03/06/10 23:57] >>463 書式設定でしたか。見落としてました。大変助かりました。ありがとうございます。 465 名前：mathmania [03/06/17 23:17] 駄スレ漏れスレage 466 名前：１３２人目の素数さん [03/06/18 00:04] あなたが探してるのってこれだよね？この中にあったよ♪ endou.kir.jp/betu/linkvp/linkvp.html 467 名前：１３２人目の素数さん [03/06/21 17:16] 4次元データを視覚化したいのですが、mathematicaに限らず 何かよい方法(ツール)はありませんか? 等高線図を複数枚、平行になるよう並べて ナナメ上から見下ろすようにすればよいと思うのですが、 よいアプリケーションを探しています。 他の方法でも、可能ならば何でもいいです。 468 名前：１３２人目の素数さん mailto:sage [03/06/22 17:01] >>467 いってることが良くわかりません。等高線図を斜めから見下ろす？ mathematicaでの実現ならば、三次元にスライスしたときの データをどうやって表示するかという話になるのだろうけど、 それはデータの形式に依存するからなんともいえないし。 469 名前：467 [03/06/23 13:26] >>468 間違いがありました。すみません。 まず画面に出したいのは4次元データ中の3次元分なので、 「3次元データを表示したい」と訂正します。 データは格子状データです。 古い人間なのでfortran形式で言うと do k=1,nz do j=1,ny do i=1,nx write(6,*) x(i),y(j),z(k),data(i,j,k) end do end do end do の形式の等間隔のデータと思って下さい。 やりたいのは、これまた解かりにくい表現かもしれませんが gnuplotで set contour base unset surface set view 60,10 くらいの時にプロットした等高線図を複数枚、一つの画像にプロットしたいのです。 470 名前：１３２人目の素数さん mailto:sage [03/06/23 20:03] >>469 gnuplot云々はわかりませんが、 ListContourPlot3Dかな? 471 名前：１３２人目の素数さん mailto:sage [03/06/24 22:42] Ver.5 ですね。 それに $1880 5.0 Eng vs \448800 日本語版 4.2only US ってこれでも Premier Service 付きだって。。。。 潮時かね 472 名前：467 [03/06/26 18:39] 忙しかったので、見るの遅くなって済みません。 >>470 >ListContourPlot3Dかな? 手持ちのマニュアルには ListContourPlotとListPlot3Dしか見付かりません。。 (ver5からの新機能ですか?) それら二つから推測すると、 ListContourPlot3D[{{{},{}},{{},{}}}] のような3重ネストの配列を入れてプロットするのでしょうか? 値段高い… 473 名前：１３２人目の素数さん mailto:sage [03/06/26 23:19] >>472 4.2にはあります。 Onlinehelpで出てくるはずです。 または、 <<Graphics` ?ListContourPlot3D 474 名前：467 [03/06/27 14:28] みなさま、御親切に 本当にありがとうございます！ ちょっといろいろ試してみます。 475 名前：１３２人目の素数さん [03/06/27 15:22] 足してnになる1以上の整数の組み合わせを並べる時、並べ方が最大になるような 整数の組み合わせをMathematica（あるいはC言語でも）で求める方法を教えてください。 わかりにくいので例を挙げますと　n=5　のとき 足して5になる1以上の整数の組み合わせは 5, 41, 32, 311, 221, 2111, 11111　であり 並べ方が何通りあるかはそれぞれ　1, 2, 2, 3, 3, 4, 1 で最大値は　2111　の場合で　4　に、 n=6の場合は　321　と　2211　の場合で最大値は　6　になります。 一般に、数字の個数（311なら3個、11111なら5個）がm個、1がk_1個、2がk_2個、…nがk_n個の 場合、並べ方はm!/(k_1!・k_2!・・・k_n!)になります。 476 名前：mathmania ◆uvIGneQQBs [03/06/27 15:36] Re:>475 最大値をどうやって計算するかは分かるだろう。 あとは、k_1,…,k_nの与え方だが、 これはk_n,k_{n-1},…,k_2,k_1の順に再帰的に与えればいいのではないか？ n=5のときを例に挙げると、 まず、k_5=1ととり、1*k_1+2*k_2+3*k_3+4*k_4=0になるk_1,k_2,k_3,k_4の組合わせは0,0,0,0しかない。 次にk_5=0ととるとき、1*k_1+2*k_2+3*k_3+4*k_4=5となるk_1,k_2,k_3,k_4を選ぶことを考える。(a) (a)の中でk_4=1ととり、1*k_1+2*k_2+3*k_3=1となるようにk_1,k_2,k_3を選ぶようにする。(k_1=0,k_2=k_3=0しかない。) (a)の中で、今度はk_4=0ととり、1*k_1+2*k_2+3*k_3=5となるようにk_1,k_2,k_3を選ぶ。(b) (b)の中で、k_3=1ととると、(k_2,k_1)=(1,0)と(k_2,k_1)=(0,2)となる。 (b)の中で、k_3=0ととると、k_2=2,1,0の3通りにわかれる。 k_2=2ととるとき、k_1=1となり、k_2=1ととるとき、k_1=3となり、k_2=0ととると、k_1=5となる。 この文では分かりにくいかもしれないが、とにかくやってみてくれ。 477 名前：474 mailto:sage [03/06/27 16:10] >>476 その考え方はわかりました。 初心者のものでプログラムの書き方を勉強してきます。 478 名前：１３２人目の素数さん mailto:sage [03/06/27 16:11] 477=475です 479 名前：１３２人目の素数さん mailto:sage [03/06/27 17:06] >>475 どうぞ。 << DiscreteMath`Combinatorica` << Statistics`DataManipulation` Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@     Partitions[5] 480 名前：１３２人目の素数さん mailto:??? [03/06/27 20:12] Mathematica Ver.5の情報はどこにあるんですか? 何が変わるのでしようか おしえてね! 481 名前：１３２人目の素数さん [03/06/27 20:24] 過去にレスがついていなかったのですが、 ますます出来るMATHEMATICA-標準偏- www.catkk.co.jp/mathcd/index.html どうでしょうか？ 482 名前：１３２人目の素数さん mailto:sage [03/06/27 20:25] >>480 www.wolfram.com 483 名前：475 mailto:sage [03/06/27 23:37] >>479 よくわからないのですが、Partitions[5]の　5　を　n　に変えれば 足してnになる1以上の整数の組み合わせの並べ方が最大に なる整数の組み合わせが求められるんでしょうか？ 484 名前：483 mailto:sage [03/06/27 23:39] すいません。”　5　を　ｎ　に”　は　”　5　を任意の整数に”という意味です 485 名前：１３２人目の素数さん mailto:sage [03/06/27 23:43] >>483 f[n_Integer] := Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@       Partitions[n] ってこと? もしかして、数学っぽいnについての式を欲しがってるんですか? それならそれは数学の問題。 プログラムで求めることはできません。 Cでなんて論外。 486 名前：479 == 485 mailto:sage [03/06/27 23:45] >>484 See f in >>485. 487 名前：475 mailto:sage [03/06/27 23:54] >>486 ありがとうございます。手元にMathematicaがないので 少し勉強して、学校ので試してみたいと思います。 数学的に（nを使って）解くことが出来たらそれが一番いいのですが、 それはここで聞くことではないので他のスレで聞いてみます。 488 名前：１３２人目の素数さん [03/07/04 13:36] 再セットアップしたらアクティベーションできなくなった 489 名前：１３２人目の素数さん mailto:??? [03/07/04 19:39] www.wolfram.com/ うわーん! Ver.4.2の情報しかないじゃないか。 しかも日本語だー 490 名前：１３２人目の素数さん mailto:sage [03/07/05 01:36] >489 どこかに English ってあるからそこから英語ページへ。 トップにあるよ。 まあ大規模行列と数値計算の速度改良が主ですね。 試行錯誤でアルゴリズム作るときは、全作業時間において コンパイル系とためはるつもりなのかな なお日本では日本語版のみしか売れない。 誰だよな日本語ヘルプキット強制バンドル計画したやつは。 491 名前：１３２人目の素数さん mailto:??? [03/07/05 21:11] ttp://www.wolfram.com/products/mathematica/newin5/ に五版の概要がありますね。 英語版でいいから新しいやつを動かしてみたいねえ。 ほんと。 492 名前：１３２人目の素数さん [03/07/05 23:02] 私が作りました♪ nuts.free-city.net/index.html 493 名前：488 mailto:sage [03/07/07 15:52] 誰か助けてくれ 494 名前：１３２人目の素数さん mailto:sage [03/07/07 19:47] >>493 サポートに連絡しろ。 495 名前：１３２人目の素数さん mailto:sage [03/07/08 01:13] >>494 わかったやってみる �dくす 496 名前：475 mailto:sage [03/07/08 16:22] >>479 Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@ Partitions[5]　　を実行したところ Part::"partd": "パート指定5[[1]]の長さはオブジェクトの深さを超えています．" 　 1/5[[1]]!　 となったのですが、直すべきところがわからないので教えてください。 497 名前：１３２人目の素数さん mailto:sage [03/07/08 21:15] >>496 ちゃんと上の二行を入力してください。 # そのエラーはPartitions[5]が評価されていないため。 498 名前：１３２人目の素数さん [03/07/08 22:23] 最大尤度法で推定をしていて、結構複雑なモデルなんですが、推定したい変数が５個 あるんです。それで、尤度関数を作って、偏微分してFindRootつかって求めようした んですが、出来ないんです。 そこで、正規分布使って、普通に正規分布を推定してみようって思って、同じような 方法を使ったんですが、これでもうまくいきません。エラーとしては {{初期値の尤度}},{{初期値の尤度}}は、｛a,b}={1,2}では長さ２のリストになって いません というような感じです。最大尤度法をやってるプログラムのサンプルどっかないかなあ って探したんですが、見つからないし・・・・。泣きそうなんですけど、どなたかお助け ください。 499 名前：１３２人目の素数さん mailto:sage [03/07/08 22:43] >>498 その説明で他人がなにかわかると思うか？ >  んですが、出来ないんです。 ふーん出来ないんだ。(どこがどう出来ないのか全くわからん) > {{初期値の尤度}},{{初期値の尤度}}は、｛a,b}={1,2}では長さ２のリストになって 文法間違いだな。 (syntax errorなんだからソース見ないとどうしようもない) >  って探したんですが、見つからないし・・・・。泣きそうなんですけど、どなたかお助け まぁ冷静になれ。 500 名前：１３２人目の素数さん mailto:sage [03/07/08 23:08] >>497 わかりました。 501 名前：498 [03/07/08 23:11] そうですか・・・・すいません。どっかサンプルプログラム探すしかない ですかね・・・ 正規分布を推定する場合、乱数を出して、その後、尤度関数を出せたとして K1[a_,b_]=D[L[a,b],a] K2[a_,b_]=D[L[a,b],b] として、 FindRoot[{K1=0,K2=0},{a,-1},{b,2}] としてるんですけど。。この部分では、間違いがないとは思ってるんですが 尤度関数が間違ってるのかな・・・・。はあー泣きそう。。。 502 名前：１３２人目の素数さん [03/07/08 23:12] ☆頑張ってまーす！！☆女の子が作ったサイトです☆ 　　　　　　　☆見て見て！！ yahooo.s2.x-beat.com/linkvp/linkvp.html 503 名前：１３２人目の素数さん mailto:sage [03/07/08 23:20] >>501 Mathematicaの基礎からやりなおしってかんじ。 とりあえずダメなところは指摘しとくが、 マニュアルか教科書(上でも挙がってる)読み直せ。 > K1[a_,b_]=D[L[a,b],a] = と := の違いが良くわかってるのなら何もいわないが。 >  FindRoot[{K1=0,K2=0},{a,-1},{b,2}] K1 = 0 は K1 という変数を 0にする。 書くなら {K1[a, b] == 0, K2[a,b] == 0} 。 504 名前：498 [03/07/08 23:20] あー勿論 FindRoot[{K1[a,b]=0,K2[a,b]=0},{a,-1},{b,2}] としてます。。。。。 サンプルプログラム載ってるサイトってないですかねえ・・・・。。。 505 名前：503 mailto:sage [03/07/08 23:22] >>504 その2行目をもういちど読み直せ。 506 名前：498 [03/07/08 23:23] :=と=の違いってえーと、すみません、良く分ってないです。。。 確か:=は、なんか固定するみたいな意味だった気が・・・・ 507 名前：498 [03/07/08 23:25] >>505 ==ですね・・・すみません、それは、やってます・・・・。 508 名前：503 mailto:sage [03/07/08 23:28] >>506 で、聞いてる方としてはあいかわらずどこがうまく行かないのか わからないんだが。 で、どううまくいかいの？ エラーメッセージは？ 509 名前：498 [03/07/08 23:37] エラーメッセージは、FindRootのところなんですよー。 >>498に書いたようなエラーメッセージが出ます。 FindRootの文法は、合ってると思うんですが。。。 510 名前：503 mailto:sage [03/07/08 23:45] >>509 498を読んだところでは第一パラグラフと第二パラグラフは 別の手段を試してるようにしか読めない(FindRootで失敗したから別の手段)。 どの関数でエラーが出力されたかもわからない。 エラーがあなたの言語に翻訳されているがそれが適切かどうかわからない。 すでに書いたがsyntax errorはソースをみないと何もわからない。 技術的なことを質問する方法をしらない、 シンタックスエラーが自分で解決できない、 レベルだったらソース張るべきだよ。 変にプライドがあるのかそういう「初心者」は多いけどね。 情報を出さずに、自分のやってることは間違ってない!、みたいな。 たんに誰かが作ったライブラリが欲しいのなら、 まず、googleなりwolframのサイトなりを検索し、 もっと適切な質問のしかたをすべき。 おれは今日はもう寝ます。 511 名前：503 mailto:sage [03/07/09 00:02] >>509 寝ようと思ってたが気づいたので。 多分 L が リストを返してるからだ。 512 名前：１３２人目の素数さん [03/07/09 14:24] ndist2=NormalDistribution[0,4]; p1=Table[ Random[ndist2],{100}]; L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}] K1[a_,b_]=D[L[a,b],a] K2[a_,b_]=D[L[a,b],b] FindRoot[{R1[a,b]==0,R2[a,b]==0},{a,{1,2}},{b,{-2,4}}] ソースは、これです。 >>511 あーそれは、どうすれば、解消できますかね。リストで返さずに関数として認識 してもらわないとつらいっすね・・・。 513 名前：498 [03/07/09 15:12] あっ最後の式は、 FindRoot[{K1[a,b]==0,K2[a,b]==0},{a,{1,2}},{b,{-2,4}}] です。 514 名前：１３２人目の素数さん [03/07/09 16:05] 今課題をやっています。 誰か偏微分係数の出し方を教えてください。 ちなみにMathematika 4.1を使用しています。 ヘルプの使い方がわからなくて・・・ 515 名前：514 [03/07/09 16:19] 解決すますた。 失礼。 516 名前：503 mailto:sage [03/07/09 17:37] >>498 ソースはコピペしてそのまま張れ。 タイプミスもあるしな。 πってPiだよな？ <<Statistics`NormalDistribution` してるんだよな？ 原因はp1のまえにある'{'だ。'{'はリスト構築子。   L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}] ↓におきかえろ。   L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-((p1[[i]] - b)^2/(2 *a^2)))],{i,1,100}] 517 名前：１３２人目の素数さん [03/07/09 21:13] このスレッドの人から見てＳ言語はどうですか？ 社会学系だとmatlabやmathematicaより使われていますが… 518 名前：498 [03/07/10 16:44] >>516 ありがとうございました。それで正規分布の推定は、できました。ただ、本当に 求めたいのは、 L[a1_, a2_, b1_, b2_, b3_] = Sum[y1[[i]]*Log[F[a1 - (b1* x1[[i]] + b2*x2[[i]] + b3* x3[[i]])]] + y2[[i]]* Log[F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])] - F[ a1 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])]] + y3[[i]]*Log[ 1 - F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])]], {i, 1, 100}]; という結構複雑なもので、これも｛を（にしてみたんですが、やはりリストで返してる みたいです。Sumの扱い方が悪いんでしょうか・・・・ てか、見にくくてすみません。 519 名前：１３２人目の素数さん mailto:sage [03/07/10 19:28] >>518 a1,a2,b1,b2,b3,F,y1,... のどれかにリストが入ってんだろ… 同じことを何度も書くが、部分だけ張るなよ… プログラミング自体初めてなのかも知れないが、レベル低すぎる。 一冊Mathematicaの本読め。 520 名前：498 [03/07/10 21:10] すみません、一遍に貼れないんで２回に分けますけど、これが全部です。。。 <<Statistics`ContinuousDistributions` ndist=NormalDistribution[0,1]; x1=Table[ Random[ndist],{100}]; x2=Table[ Random[ndist],{100}]; x3=Table[ Random[ndist],{100}]; y=Table[2*x1[[i]]+3*x2[[i]]-x3[[i]],{i,100}]; y1=Table[If[y[[i]]<=-5,1,0],{i,100}]; y2=Table[If[-5<y[[i]]<=5,1,0],{i,100}]; y3=Table[If[5<y[[i]],1,0],{i,100}]; 521 名前：498 [03/07/10 21:11] F[x_] = Integrate[Exp[z]/{1 + Exp[z]}^2, {z, -Infinity, x}]; L[a1_,a2_,b1_,b2_,b3_]=Sum[y1[[i]]*Log[F[a1-( b1*x1[[i]]+b2*x2[[i]]+b3* x3[[i]])]]+ y2[[i]]*Log[F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]- F[a1-(b1*x1[[i]]+b2*x2[[i]]+ b3*x3[[i]])]]+y3[[i]]*Log[1- F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]],{i,1,100}]; K1[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a1]; K2[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a2]; K3[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b1]; K4[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b2]; K5[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b3]; FindRoot[ {K1[a1,a2,b1,b2,b3]\[Equal]0, K2[a1,a2,b1,b2,b3]\[Equal]0, K3[a1,a2,b1,b2,b3]\[Equal]0, K4[a1,a2,b1,b2,b3]\[Equal]0, K5[a1,a2,b1,b2,b3]\[Equal]0}, {a1,0},{a2,1},{b1,1},{b2,3},{b3,1} ] 522 名前：503 mailto:sage [03/07/10 21:20] >>521 F[x_]の定義を良く読め。 それからMathematicaでプログラミングするときは 出力を確認しながらやれ。 F[x_] の行をセミコロン無しで打つか、 F[x]とか試しに打ってみながらやってれば間違えないだろ？ あとは自分でやれ。 523 名前：１３２人目の素数さん mailto:sage [03/07/11 14:51] Solve[]や、NSolve[]で、求めた複数の解のうち、 ひとつだけを使いたいのですが、 たとえばt = %[[-1]], t = Last[%]などとしても、 {x ->0.234hogehoge}と　t　に入ってしまいます。 この x->を取り除くにはどうすればいいのでしょうか？ 524 名前：523 mailto:sage [03/07/11 15:00] x /. % でできました。トンクス 525 名前：１３２人目の素数さん mailto:sage [03/07/12 00:43] >521 えー、 Sum[<<>>,{i,1,100}] ってあるけど 数値やるなら NSum じゃ駄目なん？ どうせFindRoot なんだし 526 名前：１３２人目の素数さん [03/07/14 23:27] あのー四捨五入するような関数ってないですか？切り捨てでもいいんですけど。 マニュアルで調べてみたんですが見つからなくて・・・・ 527 名前：１３２人目の素数さん [03/07/14 23:28] Roundじゃないの？ 528 名前：１３２人目の素数さん [03/07/14 23:36] 0.5足して切り捨てろ。 529 名前：山崎 渉 mailto:（＾＾） [03/07/15 12:38] 　__∧＿∧_ 　|（　　＾＾ ）|　＜寝るぽ（＾＾） 　|＼⌒⌒⌒＼ 　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉 　　 ~￣￣￣￣ 530 名前：１３２人目の素数さん [03/07/16 23:06] shop.ponycanyon.co.jp/kids/p14796.html 531 名前：１３２人目の素数さん [03/07/17 07:26] age 532 名前：１３２人目の素数さん [03/07/17 14:11] >>530 これなら3,800円か。安いな。 533 名前：１３２人目の素数さん [03/07/17 16:26] 質問なんですが、例えば x=k/10 y=k+1 For[k=1,k<11,k++, 　　r=x+y] みたいに繰り返し計算をさせる場合に、このrは、毎回更新されてしまうんで その結果を格納してリストの形にしたいとするとどーすればいいですか? 勿論、これぐらい単純だとTable使えばいいんですが、やりたいのはもっと ややこしいんで、結果をリストにする格納する方法を教えてください。 よろしくお願いします。 534 名前：１３２人目の素数さん mailto:sage [03/07/17 19:13] >>533 単純であろうがなかろうがTableでいいはずだ。 どうしてもというなら l ={} For[k=1,k<11,k++,r=x+y;l = Prepend[r,l]] 等。 しかしmathematicaではFor文を使う必要はない(その上 遅いので使うべきでもない)ので再考することを勧める。 535 名前：１３２人目の素数さん mailto:sage [03/07/18 14:39] とぇ 536 名前：１３２人目の素数さん [03/07/22 12:25] FindMinimumを使って最小値とそのときの変数の値を求めると、 {234,{a -> 3,b -> 2}} みたいな感じだと思うんですが、これをいっぱい出して、a,bの平均値を出すには どうすればいいでしょうか？ 537 名前：１３２人目の素数さん mailto:sage [03/07/23 01:12] >>546 漠然としすぎ。 こういうことか？ << Statistics`DescriptiveStatistics` data = {x, y} /. (#[[2]] &) /@         Flatten[Table[             FindMinimum[Sin[x + y], {x, i}, {y, j}], {i, 0, 3}, {j, 0, 3}],           1]; Mean /@ Transpose @ data 538 名前：１３２人目の素数さん mailto:sage [03/07/23 01:12] >>537 s/546/537 539 名前：537=538 mailto:sage [03/07/23 01:13] 536駄・・・まちがえすぎだ。 540 名前：536 [03/07/24 13:07] すみません、537のプログラム見て解読してたんですけど、ちょっと良くわかんなくて。 たとえば、 FindMinimum[f[a,b],{a,1},{b,2}] とした場合、{最小値,{a -> 0.89898},{b->0.1}] みたいにした場合、a->の部分を消すのはどうすればいいですか？ 541 名前：536 [03/07/24 14:36] なんか読みにくいですね・・えーと置換なんかでうまくa->の部分を消せない かなあと・・・無理なんですかねえ。。。 542 名前：１３２人目の素数さん [03/07/24 14:37] 開店セール見に来てーーーーーーーーーー！DVDを見ようー！GO！ ▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼ 開店セール1枚　500円より→→→www.dvd-exp.com　 アダルトDVDが激安販売→→→www.dvd-exp.com　 アダルトDVDの価格はなーんと1枚900円→→→www.dvd-exp.com　 最短でご注文いただいた翌日にお届け！→→→www.dvd-exp.com　 セルビデオ店の人気商品が24時間お買い物が出来ます。 見やすい画面・注文簡単・時間をかけてゆっくり選んでください！ 　　　　　　　　　　www.dvd-exp.com←へGO！ 　　　　　　　　　　　DVDエクスプレス 543 名前：１３２人目の素数さん mailto:sage [03/07/24 15:18] >>541 簡単に消せるように"->"つきで出力されている。 s = {a -> 1, b -> 2} {a,b} /. s 544 名前：１３２人目の素数さん [03/07/24 22:08] s はリストみたいになっているから s[[1]] や s[[2]] とかで取り出しても良いんじゃない？ 545 名前：１３２人目の素数さん mailto:sage [03/07/25 04:10] >>544 543で解決してます。”a->”みたいなのを取り除きたいってのが541の質問だろ。 ってか、房は黙ってろ（ｗ

---

---

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

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