Mathematicaプログラミング 質問箱 その1
■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん
NGNG本家
http://www.wolfram.com/index.ja.html
http://www.jip.co.jp/si/soft/mathematica/
Mathematicaユーザー会（有料）
http://www.jp-mathusers.org/
0278デフォルトの名無しさん
2009/12/26(土) 17:30:36u[x_] := x^(2/5)
Plot[u[x], {x, 0, 5}, AxesLabel -> {"x", "u(x)"}]
上記でかいたグラフに点(a,u(a))における接線を書き加えるModuleを作りたいんです
tangent[a_] := Module[u'[x]*(x - a) + u[a]]
としてtangent[1.5]などと入れるとそのaの値における接線が書き加えられたグラフを描くようにしたいのですが
うまくいきません
0279デフォルトの名無しさん
2009/12/28(月) 15:27:47Module は使わなかったのですが、
u[x_] := x^(2/5)
g1 = Plot[u[x], {x, 0, 5}, AxesLabel -> {"x", "u(x)"}]
tangent[a_] := Show[g1, Plot[u'[a]*(x - a) + u[a], {x, 0, 5}]]
でやりたいことが実現されるでしょうか?
0280デフォルトの名無しさん
2010/01/13(水) 12:37:31u[x_] := x^(2/5)
tangent[a_] := Module[a],
Plot[u[x], u'[a]*(x - a) + u[a], {x, 0, 5}, AxesLabel -> {"x", "u(x)"}]
こんな感じにしたいのですができませんでした・・・
0281デフォルトの名無しさん
2010/01/16(土) 19:07:44こんな感じでいけます。
u[x_] := x^(2/5)
tangent[a_] :=
Plot[{u[x], u'[a] (x - a) + u[a]}, {x, 0, 5}, AxesLabel -> {"x", "u[x]"}]
tangent[1]
0282デフォルトの名無しさん
2010/01/21(木) 18:20:044次元データのcontourを書きたいのです．
ControurPlot3Dを使えば良さそうなのは分かるのですが，使い方がわかりません．
陰関数f(x,y,z)=0の球面上にg(x,y,z)の大きさでcontourを描きたいんです．
http://tech.ckme.co.jp/img/sphere8cdata_matlab.png
こんな感じです
よろしくお願いします！
0283デフォルトの名無しさん
2010/02/12(金) 11:47:20とりあえず以下のプログラムになりました。
In[1] := Epidemic[n_, p_] :=
Module[{choices = {{1, 0}, {0, 1}, {-1, 0}, {0, -1}},
reject, pickAndChoose, select, newpers},
reject = {};
pickAndChoose :=
(select = #[[2, Random[Integer, {1, Length[#[[2]]]}]]];
If[Random[] <= p,
newPers =
Complement[Union[Map[Function[y, y + select],
choices], #[[2]]],
{select}, #[[1]], reject];
{Join[#[[1]], {select}], newPers},
reject = Join[reject, {select}];
{#[[1]], Comple, ent[#[[2]], {select}]}]) &;
FixedPoint[pickAndChoose, {{{0, 0}}, choices}, n,
SameTest -> (#2[[2]] == {} ||
Length[#2[[1]]] == n &)][[1]]
]
上を実行すると
SetDelayed::write : In[1]のタグInはProtectedです．
と、エラーが出るんですけど、どうしたら解決しますか？
初心者でヘルプを見てもよくわからなくて・・・わかる方いましたら教えてください。
0284デフォルトの名無しさん
2010/03/09(火) 18:00:420285デフォルトの名無しさん
2010/03/29(月) 20:23:05なるほど。ありがとう
0286デフォルトの名無しさん
2010/06/09(水) 14:23:01ただし、東西南北をそれぞれ、｛１、０｝｛ー１、０｝｛０、ー１｝｛０、１｝で
表すことにする。
（ヒント：りすとdirectotions={{1,0},{-1,0},{0,-1},{0,1}}を作り、1、2、3、4の範囲の
乱数kを発生させ、directotionsのk番目の要素をとるとランダムな方角が得られる）
すいませんが練習だと思って解いてくれませんか？課題が出せないんでお願い
0287デフォルトの名無しさん
2010/06/10(木) 08:30:170288デフォルトの名無しさん
2010/06/14(月) 17:32:48例えば、f4[1234]=10.
わかる人います？教えてくれませんか？
0289デフォルトの名無しさん
2010/06/14(月) 20:51:23ヒント：1234 をなんとかして {1,2,3,4} の形にして足す。
0290デフォルトの名無しさん
2010/06/16(水) 13:53:39正整数nのケタ数を与える関数f2を書け。たとえば、f2[40172]＝５
0291デフォルトの名無しさん
2010/06/16(水) 14:57:57問題間違えました。
０より大きく１より小さい既約分数で分母がn以下のものを
小さい方から並べてできるリストを出力する関数f5[n]を書け。たとえば、
f5[n]={1/4,1/3,1/2,1/3,1/4}
（このような数列をnに対するFarey数列といい、いろいろおもしろい性質が知られている。
たとえば、隣り合う項をa/b,c/dとすると、必ずab-cd=-1となっている。）
すいませんがお願いします。
0292デフォルトの名無しさん
2010/06/16(水) 15:22:37あまり美しくないが。
f5[n_] := DeleteDuplicates[Sort[Flatten[Table[Table[Which[0 <= a/b <= 1, a/b, True, {}], {a, 0, n}], {b, 1, n}]]]]
0293デフォルトの名無しさん
2010/06/16(水) 15:43:480294デフォルトの名無しさん
2010/06/16(水) 16:06:39たとえば、f[1234]＝4321
r(n)=n/4となるような最小の正の整数を求めよ。
お願いします。意味わかりません。
0295デフォルトの名無しさん
2010/06/16(水) 17:21:11f[1234]＝4321
を作るとどうなる？
0296デフォルトの名無しさん
2010/06/16(水) 18:04:330297デフォルトの名無しさん
2010/06/16(水) 18:34:320298デフォルトの名無しさん
2010/06/18(金) 01:07:450299デフォルトの名無しさん
2010/06/18(金) 02:08:260300デフォルトの名無しさん
2010/06/18(金) 09:03:39もっと頭を柔らかく、Mathematica 的に考える。
Mathematica が得意なのは何か？
ヒント：>>288
0301デフォルトの名無しさん
2010/06/23(水) 13:58:45{3,1,1,5,2,1,2}み適応し、{3,1,5,2}になることを確かめよ。
（Unionを使うと{1,2,3,5}になるので注意）
なんか方法ある？
0302デフォルトの名無しさん
2010/06/24(木) 19:41:09「その数は2通りの仕方で2つの正の立方数（立方数とはn^3）の和で表される
自然数のうち最小のものです」さてその数は何か？
答えはたしか９１になればいいと思います。
僕にとってムズイので教えてください。
91 = 6＾3 + (-5)＾3 = 4＾3 + 3＾3
0303デフォルトの名無しさん
2010/06/25(金) 13:18:07「その数は2通りの仕方で2つの正の立方数（立方数とはn^3）の和で表される
自然数のうち最小のものです」さてその数は何か？
答えはたしか1729になればいいと思います。
僕にとってムズイので教えてください。
0304デフォルトの名無しさん
2010/06/30(水) 12:58:266より大きな偶数nに対し、n=p+q(p<=q)となるような素数の｛p,q｝を
与える関数goldbach[n]を書け。
goldbach[1000]を実行せよ。
goldbach[n_] := Table[{i, n - i}, {i, 3, n/2}]
までやったんだけど、goldbach[10]をやると、
{{3, 7}, {4, 6}, {5, 5}}になるんだ。どうすりゃ｛4，6｝消える？
PrimeQはどこで入れるの？最初から違う
0305デフォルトの名無しさん
2010/06/30(水) 14:51:49例えば１２３２１は対称数である。
正整数nに対し、それが対称数ならnを出力し、
そうでなければ、nとその逆転数を加え、それが対称数かどうか調べる。
この操作を対称数が得られるまで繰り返し、
得られた対称数を出力するような関数f[n]を書け。
たとえば、
78→78+87=165→165+561=726→726+627=1353→1353+3531=4884
なのでf[78]=4884である。
わかんないです。助けてください。
0306デフォルトの名無しさん
2010/07/02(金) 13:58:34そして、その数値と求積の公式や積分を利用して得られる理論値との違いを調べてください。
すなわち答えることは以下の4点です。
（１）どうゆう図形について考えましたか？
（２）その図形の面積をモンテカルロ法で求めるとどうゆう結果がでましたか？
（３）その図形の面積を数学で求めるとどうなりますか？
（４）あなたがモンテカルロ法で計算した結果は理論値に対して何％のずれがありますか？
0307デフォルトの名無しさん
2010/07/06(火) 00:48:57トランプのカードを
{H,1},{H,2}.....,{H,13}（ハート）
{S,1},{S,2}.....,{S,13}（スペード）
{D,1},{D,2}.....,{D,13}（ダイヤ）
{C,1},{C,2}.....,{C,13}（クラブ）
で表す。
（１）52枚のカードからランダムに5枚選び出す関数を書け。
（２）5枚の手札の（ポーカーにおける）点を下の表に従って与える関数を書け。
（３）（１）を1000回繰り返したときの平均点数を求めよ。
（表）
ワンペア 5点 ツーペア 10点 スリーカード 15点
ストレート 20点 フラッシュ 25点 フルハウス 30点
フォーカード 40点 ストレート・フラッシュ 50点
ロイヤル・ストレート・フラッシュ 100点
0308デフォルトの名無しさん
2010/07/07(水) 09:58:37solve(A.x=b,x)
0309デフォルトの名無しさん
2010/07/07(水) 12:27:580310デフォルトの名無しさん
2010/07/07(水) 18:54:410311デフォルトの名無しさん
2010/07/07(水) 19:26:37俺へたっぴだから。
0312デフォルトの名無しさん
2010/07/27(火) 00:21:29Hoge[Erf]ってすると、∫e^t2dtって感じの出力が得られるような関数です。
0313デフォルトの名無しさん
2010/07/27(火) 03:03:380314デフォルトの名無しさん
2010/07/27(火) 03:22:22なるほど。ありがとうございます。
0315デフォルトの名無しさん
2010/09/29(水) 23:06:460316デフォルトの名無しさん
2010/12/22(水) 21:05:11.nbファイルをもらったんだけど，
mathematicaでは簡単な計算をしたことはありますが，
プログラムをしたことがありません．
0317デフォルトの名無しさん
2010/12/22(水) 21:07:210318デフォルトの名無しさん
2010/12/22(水) 22:03:500319デフォルトの名無しさん
2011/03/09(水) 09:56:08.23によると
Names["List*"]
ってやるとListで始まるシンボルが列挙されるらしいぞ。
http://reference.wolfram.com/mathematica/ref/Information.html
によると
シンボルが分かったら
? シンボル
でいいらしいぞ。
0320デフォルトの名無しさん
2011/04/25(月) 01:19:47.590321デフォルトの名無しさん
2011/04/25(月) 08:22:31.65「Mathematica Cookbook 日本語版」出版記念講演
http://www.hulinks.co.jp/software/mathematica/seminar_cookbook.html
4月27日(水)14:00〜17:00
参加費無料
0322天使 ◆uL5esZLBSE
2011/07/03(日) 05:25:19.66変数に＄ついてる言語触ってるって事だよね
いちいちSHIFT＋４キーおして ＄ 打ちまくってる感触はどう？
はい次（笑）
0323天使 ◆uL5esZLBSE
2011/07/03(日) 18:09:10.72＜＜＜＜＜＜＜＜＜＜＜ こんなのがあるよ。 ＞＞＞＞＞＞＞＞＞＞＞(キリッ!キリッッッ！！ッ！！！
∨∨
＜＜＜＜＜＜＜＜＜＜＜＜ 「Mathematica Cookbook 日本語版」出版記念講演 ＞＞＞＞＞＞＞＞＞＞＞＞(キリッッッッ!キリッッッッ！！！！キリッッッッッ！
∨
＜＜＜＜＜＜＜ 4月27日(水)14:00〜17:00 ＞＞＞＞＞＞＞（←キリッ!!!きリッッ！！！キリ!!キリ!!ッッ！！
ゴミじゃねーか
0324天使 ◆uL5esZLBSE
2011/07/04(月) 19:05:29.37ここで、Ruby以外を選ぶ奴ってマジでなんなんだろうな
ゴミだよなぁほんと
0325デフォルトの名無しさん
2011/07/05(火) 02:02:10.53ペーペー野郎を助けて頂けませんか？
Mathematicaって、けっこうバグとかあったりするものです？
非線形方程式、非線形常微分方程式とか解かせているんですが、
確認のために、解を元の方程式に代入しても方程式が成り立た
なくなってしまうんです・・・
何か特別なことをしないといけないのですかねぇ・・・。
0326デフォルトの名無しさん
2011/07/05(火) 02:58:43.71〓 Mathematica 5 〓
http://kamome.2ch.net/test/read.cgi/math/1285859504/
0327デフォルトの名無しさん
2011/08/11(木) 15:38:22.06みたいなプロットできる？
0328デフォルトの名無しさん
2011/09/06(火) 13:52:32.94ListPointPlot3Dが近いかも
■ このスレッドは過去ログ倉庫に格納されています