〓 Mathematica 2 〓
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
NGNG, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .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/
0369１３２人目の素数さん
NGNG>>329で動かないということは、Mathematicaのバージョンが違うのだろうか。
特異点もどきの場所が正確に特定できるなら>>368の一行目が利用できるが、
そうでないなら、オプションにMethod -> MonteCarloを指定すると
解決するかもしれない。
0370１３２人目の素数さん
NGNGf[x_] := If[x == 0, 0.0, Sin[x]/x];
NIntegrate[f[x], {x, -1, 1}]
とすれば大丈夫
0371１３２人目の素数さん
NGNGx -> 0 で Sin[x]/x -> 1 だから、
g[x_]:=If[x==0,1.0,Sin[x]/x]
としてみると
f[x_]:=If[x==0,0.0,Sin[x]/x]
の積分と僅かに違います
正しいグループは
N[Integrate[Sin[x]/x,{x,-1,1}]]
NIntegrate[g[x],{x,-1,1}]
別なグループは
NIntegrate[Sin[x]/x,{-1,0,1}]
NIntegrate[f[x],{x,-1,1}]
でグループ間の差は 2.2*10^(-16) でした。
面白いですね。
0372１３２人目の素数さん
NGNG面白いか？
アルゴリズムの都合でちょっと誤差がでただけだろ
(どっちも通常の意味の積分値は同じ)。
誤差が少ないところは流石なのかもしれないが。
0373１３２人目の素数さん
NGNG0374１３２人目の素数さん
NGNG0375１３２人目の素数さん
NGNGf[x_]:=If[x==0,1.0,Sin[x]/x] の間違いでした。すんません。
0376すいません教えてください
NGNG[0,1]でe＾xの積分の近似値を求めその値を実際の積分の値と比較せよ。
ってゆうもんだいなんですけどどうやったらいいんですか？誰か教えてくださいませんか？？おねがいします
0377１３２人目の素数さん
NGNG0378すいません教えてください
NGNG0379KingOfKingMathematician ◆H06dYXOYLA
NGNGなくして困っている
0380１３２人目の素数さん
NGNG0381１３２人目の素数さん
NGNGContourPlotで行けるんでね？
0382１３２人目の素数さん
NGNGバーカ、市ね。くたばれ
0383１３２人目の素数さん
NGNGCos[x/5-1]
]
いじょ
0384１３２人目の素数さん
NGNG0385１３２人目の素数さん
NGNG最終的に" ° "で答えが欲しいのですが...
今日体験版を落として、触っているのですが分からなくて・・・
0386１３２人目の素数さん
NGNG0387385
NGNG0388385
NGNGArcTan(163/217)の答えを出したいのですが
答えも ArcTan(163/217)とでて" ° "で答えが欲しいのですが...
0389１３２人目の素数さん
NGNG\!\(N[ArcTan[163\/217]\ 180\/π]\)
0390385
NGNG0391１３２人目の素数さん
NGNG0392１３２人目の素数さん
NGNG0393１３２人目の素数さん
NGNG５．０proにup gradeすることに決めました。
そこで聞きたいのですが
今現在パソコンを買い換えたのでmathematica student版が
インストールされていません
こういう場合でもアップグレードできるのでしょうか？
0394１３２人目の素数さん
NGNG大丈夫。
とりあえずここら辺でも見ておけ。
ttp://www.varsitywave.co.jp/products/mathematica/html/math-upgrade.html
0395１３２人目の素数さん
NGNGy軸として y={y1,y2,y3,....}
という感じにデータセットがあるんですが（要素数は同じ）、
これを
test = { {x1,y1}, {x2,y2},},....} といった感じのx,yを対応させたリストに
するにはどうしたらいいでしょうか。すべてただの実数です。
どなたかご教示願います・・・。
0396１３２人目の素数さん
NGNG0397１３２人目の素数さん
NGNGそれはできますね
マニュアルかヘルプを読みましょうよ
0398１３２人目の素数さん
NGNG0399１３２人目の素数さん
NGNGありがとうございました。
簡単にできました。申し訳ないです。
0400１３２人目の素数さん
NGNG0401１３２人目の素数さん
NGNGＧＵＩＫＩＴが出ています。
公式ホームページから
ライセンスナンバーだけでダウンロードできます。
インストールする場所と
ヘルプファイルの更新をお忘れなく。
0402１３２人目の素数さん
NGNGCのプログラムで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超初心者ですが、
ご教授いただければと思います。゜゜(´□｀｡)°゜。
0403FeaturesOfTheGod ◆UdoWOLrsDM
NGNGもしそうでないとしても、固有値の計算は時間がかかる。(少なくとも私の知っているアルゴリズムでは。)
0404１３２人目の素数さん
NGNG読み込みにではなく、表示に時間がかかっている。
どうしても1000*1000の行列を眺めたいのでなければ
ReadList["matrix.txt",Number,RecordLists->True];
Eigenvalues[%]
とするといいだろう。
0405１３２人目の素数さん
NGNGこの角度のパラメータをラジアン表記でなく度数表記で一発出力させる
ためには、どういじればいいのでしょうか。
<< Calculus`VectorAnalysis`
SetCoordinates[Cartesian[x, y, z]];
x = 1.;
y = 1.;
z = 1.;
CoordinatesFromCartesian[{x, y, z}, Spherical]
0406１３２人目の素数さん
NGNGCoordinatesFromCartesian[{x, y, z}, Spherical] / Degree HoldForm[Degree]
0407１３２人目の素数さん
NGNG>>406
ありがとうございます。実行してみました。
返数のうち最初の一つは長さなので、角度である後二つにDegreeを
かけたいのですが、方法はあるのでしょうか。
すいません、初歩的なことで。
0408１３２人目の素数さん
NGNGすまん。こっちこそ初歩的に間違えた。
CoordinatesFromCartesian[{x, y, z}, Spherical] {1, HoldForm[Degree] / Degree, HoldForm[Degree] / Degree}
# ラジアンに戻す時は
# x = Coord...<snip>
# ReleaseHold[x]
0409１３２人目の素数さん
NGNG>>408
ありがとうございます。できました。
勉強させていただきます。
0410402
NGNGありがとうございました。それで解決いたしました。
そんなことだっとは盲点でした……。
次のような記事も発見しました。
私の場合ここまでする必要はありませんでしたが、参考までに。
『ReadList: Is it really slow?』
ttp://forums.wolfram.com/mathgroup/archive/2000/Jul/msg00264.html
0411４１１
NGNG同じ座標内に全て表示するようにしたいんですけど、どうすればいいんでしょうか？
0412１３２人目の素数さん
NGNGSelectionAnimate か、
表示されたものをマウスで畳んで、
右クリックしてAnimateみたいなメニューを選ぶ。
0413１３２人目の素数さん
NGNGギリシャ文字はマウスを使ってボックスから選択しないと使えないのでしょうか。
キーボードから入力するショートカットがあれば教えてください。
釣りではないのでよろしく。
0414１３２人目の素数さん
NGNGぜってー釣りだな
0415１３２人目の素数さん
NGNG\[CapitalAlpha]
\[Alpha]
<Esc>Alpha<Esc>
<Esc>alpha<Esc>
0416FeaturesOfTheGod ◆UdoWOLrsDM
NGNG0417４１１
NGNGそれでShowをつかって複数のグラフィックスを同一座標内に表示したように
グラフィックスを重ねて表示する事ができるんでしょうか？
もうちょっと具体的に教えてくれるとうれしいんですけど・・
0418１３２人目の素数さん
NGNG例えば、ドラえもんとかドラゴンボールとか…
0419FeaturesOfTheGod ◆UdoWOLrsDM
NGNG0420１３２人目の素数さん
NGNGアニメーションがしたいわけではなくて、単に重ねたいのなら
最初からShowか元のPlotなどで重ねろ。
Animatieの出力を利用するのはNotebook操作系の関数を使えば
できないことはないが、二階から目薬。
0421４１１
NGNG二回から目薬とは？？
どのようなNotebook操作系の関数を利用すればできますか？
0422１３２人目の素数さん
NGNG0423１３２人目の素数さん
NGNG0424４１１
NGNG0425１３２人目の素数さん
NGNG{(b-f+b*X^2)*e}/{X*(d+d*X^2+e)}==0,X
0426１３２人目の素数さん
NGNG0427１３２人目の素数さん
NGNG0428１３２人目の素数さん
NGNG0429１３２人目の素数さん
NGNG0430１３２人目の素数さん
NGNGというか割り当ての方法をどなたかご存じないですか？
例えば
コントロール＋ｋ
を行末までの文字の削除に当てるとか。
0431１３２人目の素数さん
NGNGなんか大きなことやろうとすれば、matlab 使うなー。
mathematicaは言語の外見が汚なすぎ。
てかセンス悪いよな、Wolframは
0432１３２人目の素数さん
NGNG0433FeaturesOfTheGod ◆UdoWOLrsDM
NGNGそもそも、mathematicaは研究のためにある。
0434１３２人目の素数さん
NGNG0435FeaturesOfTheGod ◆UdoWOLrsDM
NGNG0436１３２人目の素数さん
NGNG>mathematicaは言語の外見が汚なすぎ。
>てかセンス悪いよな、Wolframは
そこはかなり同意する
0437１３２人目の素数さん
NGNGある程度のプログラムを組めるというのを両立させようとするから
ああいうことになるんだよな。
0438１３２人目の素数さん
NGNGそうですね。mathematicaは開発者が多いからグラフィックスやGUIだけが
他のソフトより良いだけであって、Wolfram本人もキモいし、言語もキモいし、
やたらプログラミングに時間がかかるだけですね。
もうmathematicaは逝ってよしだと思いました
0439１３２人目の素数さん
NGNGがんばってください。キモっ!
0440１３２人目の素数さん
NGNGRubyとかは特にね。
matlab言語も（・∀・）ｲｲ!!よ
0441１３２人目の素数さん
NGNGMATHEMATICA撲滅しよーぜ
0442１３２人目の素数さん
NGNGMathematica氏ね。
0443１３２人目の素数さん
NGNG0444１３２人目の素数さん
NGNGどうやったら元の関数として機能してくれるの？
0445FeaturesOfTheGod ◆UdoWOLrsDM
NGNG0446FeaturesOfTheGod ◆UdoWOLrsDM
NGNG0447１３２人目の素数さん
NGNG評価の仕組み上無理だったと思います。
どうしてもというなら、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を使いましょう。
# 前身時代からなんで知識があるわけでもないのに考えずに適当なこと
# ばかりこたえるかな。
0448１３２人目の素数さん
NGNG>>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に方程式の解が代入されただけの感じのものです。
どうすれば成功するのでしょうか。アドバイスもらえないでしょうか。
0449１３２人目の素数さん
NGNGTable[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}
となるだけ。
なんでなのでしょうか。
0450１３２人目の素数さん
NGNGCoordinatesFromCartesian[{{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
と入力すると何がおこるでしょう？
0451１３２人目の素数さん
NGNGあうう…
「評価」じゃなかったのですね。
思いこみであちこちにEvaluateつけて考えあぐねておりました
Thanks a lot!
0452１３２人目の素数さん
NGNGMathematicaでの数値計算の精度は機械精度以上には上げることが出来ないでしょうか。
$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の使い方がまずいのでしょうか。
0453１３２人目の素数さん
NGNGSqrt[2 - 0.1]の結果の数値を精度100桁にするというのと、
(2. - 0.1)の結果の数値を精度100桁としてからSqrtをとるのとは
全く別のことですが。
もちろんどちらもN[Sqrt[2 - 1/10], 100]とは異なりますね。
0454１３２人目の素数さん
NGNG0455１３２人目の素数さん
NGNG0456１３２人目の素数さん
NGNGひょっとして ０．１と１／１０が同じと思っていませんか？
0457１３２人目の素数さん
NGNG一緒と考えておりました...
これで納得できました。ありがとうございました。
0458１３２人目の素数さん
NGNG0459１３２人目の素数さん
NGNG発生させた乱数のリストから任意の数をランダム抽出するにはどうしたらいいでしょうか？
0460１３２人目の素数さん
NGNG乱数から「ランダム抽出」って意味が無いだろ。
Table[YourRandomFunction[], {i, 10}]
とでもすればいいんじゃないか？
0461Not, used with Select
NGNGたとえば
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 というのを使いました．もっとマシな方法あれば知りたし．
0462１３２人目の素数さん
NGNG! PrimeQはNot[PrimeQ]ですね。
Not[PrimeQ][x]とNot[PrimeQ[x]]は違いますね。
Not[PrimeQ[#]]& か Composition[Not, PrimeQ]を使ってください。
!MemberQなんてものをいったいどう使ったのか興味津津です。
0463Not, used with Select (unresolved)
NGNG>!PrimeQはNot[PrimeQ]ですね。
0-1逆の特性関数を作るわけではないようですね．
普通の関数が否定できないのは分かるが，述語関数は否定できてもまずくなさそうに思えますが．
>Not[PrimeQ][x]とNot[PrimeQ[x]]は違いますね。
前者は意味をなさないということなんですかね．
[ところで同じ意味をなさない物でも，
! PrimeQ[3]だと何も出ないが(! PrimeQ[3])だと False と出るのは，
前者は括弧の位置が曖昧だが，後者は(! PrimeQ)[3]ではないと見なされる
ということなのだろうか．一番外側に括弧がつくだけで差が出るのは妙だが．]
>Not[PrimeQ[#]]&かComposition[Not,PrimeQ]を使ってください。
ありがとうございました．続く．
0464Not, used with Select (unresolved)
NGNGSelect[{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}を出しました．
0465１３２人目の素数さん
NGNG> ! PrimeQ[3]だと何も出ないが(! PrimeQ[3])だと False と出るのは，
!が「行の先頭」に来た時はshellでのコマンド実行を意味する。
> 等しくなる可能性があるから a=b, a=c などは False と言えないということでしょうか．
そのとおり。===を使えばいい。
> Select[{b, a, d, c, e}, ! MemberQ[{c}, #] &]でも
それはあたりまえだね。
偶然、Not[MemberQ[ ]]& という形で使っていたわけだ。
0466sage
NGNG>>465
行けそうです．ありがとうございました．
===は初耳でした(Mathematica 使い始めてまだ2週間くらいなもので)ので，
その周辺のコマンドを勉強しておきます．
0467１３２人目の素数さん
NGNG0468１３２人目の素数さん
NGNG|x|,|y|<=3 において z>=0 の部分の地形を等高線で表せ。
という課題が出まして、
ContourPlot[-x^2 + y + 1, {x, -3, 3}, {y, -3, 3}, Contours -> 10,
PlotRange -> {0, 4}
といった具合にやらせてみたんですが、これだと真っ黒で見た目が非常に悪いんです。
色を使って高さを表すにはどうしたらよろしいでしょうか？
どなたかアドバイスお願いします。
■ このスレッドは過去ログ倉庫に格納されています