〓 Mathematica 3 〓
レス数が900を超えています。1000を超えると表示できなくなるよ。
0001132人目の素数さん
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/
0821132人目の素数さん
2007/08/14(火) 10:41:56ありがとう。
0822132人目の素数さん
2007/08/17(金) 18:11:45(互換性を保ってるのがスゴイ)
0823132人目の素数さん
2007/08/17(金) 19:01:11確かに凄いが、値段がそれ以上に凄い!
個人が買うにはあまりにもきびしい、、、。
0824132人目の素数さん
2007/08/18(土) 10:28:51バージョンアップだけでも大変。
0825132人目の素数さん
2007/08/18(土) 10:47:516は使ってないけど、5.2ですでにそうだった。インストール過程で
mathIdとかいうPC固有の番号みたいのだして、それをwolfram researchへ
おくるとそれに合わせた文字列を送ってきて、それを入れると動く。
0826132人目の素数さん
2007/08/18(土) 10:50:16苛々する
あまのじゃく
くそったれ(`o´)
●●●らめ!!
あ〜ちょっとすっきり(--;)失礼しました
0827132人目の素数さん
2007/08/18(土) 12:30:28それはStudent版ではなくて一般の個人が買う40万円のパッケージでも同じ?
0828132人目の素数さん
2007/08/19(日) 11:59:57少なくともStudent版からPremierにアップグレードしたパッケージでは同じだった
0829132人目の素数さん
2007/08/20(月) 13:31:440830132人目の素数さん
2007/08/21(火) 19:59:000831132人目の素数さん
2007/08/21(火) 20:26:02似たようなフリーソフトもあるし。
0832132人目の素数さん
2007/08/23(木) 16:40:25例えばSinh[5]とか入れても入力と同じ答えが
出力されます。
どうすれば具体的な値を計算することが出来るのでしょうか?
0833132人目の素数さん
2007/08/23(木) 17:29:180834132人目の素数さん
2007/08/23(木) 18:59:240835132人目の素数さん
2007/08/23(木) 19:28:510836132人目の素数さん
2007/08/23(木) 20:18:320837132人目の素数さん
2007/08/24(金) 00:36:06個人で買うほどの>>835がきっと面白いネタを投下してくれる。
0838132人目の素数さん
2007/08/24(金) 12:23:44と書くと余計にネタっぽいが、小田急の新宿西口店に売ってるんだよ。
最近は休日に動画作って遊んでる。
0839132人目の素数さん
2007/08/24(金) 12:34:150840132人目の素数さん
2007/08/24(金) 12:35:000841132人目の素数さん
2007/08/24(金) 15:48:26>>835がネタだなんて言ってないよ。個人の趣味で買うくらいだから
無理矢理使わされてる学生と違ってちゃんと使ってるんだろうし、
きっと面白い話のネタを投下してくれるだろうと。
0842842
2007/08/24(金) 19:22:250843132人目の素数さん
2007/08/24(金) 21:26:49が変化してなくて残念。
StringForm["``", 1`4] とか入力すると
1.`3.9999999999999996
とか珍妙な結果になるんだよなー。まあMathematica的にはもっともな理由がある
のかも知れんが、ここはやっぱり「1.000」と表示して欲しいところ。
Wolframの中の人はStringForm使わんのかな。レポート書きでは実数をフォーマッ
トしてプリントしたいことが頻繁にあるんだが。
0844132人目の素数さん
2007/08/24(金) 21:29:190845132人目の素数さん
2007/08/24(金) 22:54:34835じゃないけど、じゃあ夏休みの宿題でやったネタを投下しよう。fooが何をす
る関数かわかるかなー。(実行すると一発でわかるけど)
foo[n_] :=
If[n <= 1, {{1}},
Module[{m = foo[n - 1],
h = ({#} &) /@ Range[(n - 1)^2 + 1, n^2 - n],
w = Range[n^2 - n + 1, n^2]},
If[EvenQ[n], Prepend[Join[m, Reverse[h], 2], Reverse[w]],
Append[Join[h, m, 2], w]]]];
0846132人目の素数さん
2007/08/24(金) 23:42:36君がやりたいのはこうだろ。
StringForm["``", OutputForm[1`4]]
0847132人目の素数さん
2007/08/25(土) 00:27:00いや、引数にHogeFormをかましてやれば良いのは分かってるけど、書式付き出力
書くたびにいちいちNumberFormとかScientificFormとかEngineeringForm
とか引数毎に10数文字もタイプするのがクソやってらんねーと言う話です。
凝ったオブジェクトならともかく実数くらいデフォルトでそれなりに整形してく
れてもバチはあたんないと思うんだけどなー。
0848132人目の素数さん
2007/08/25(土) 00:31:170849132人目の素数さん
2007/08/25(土) 04:45:19君がどの書式が欲しいのかは君しか知らないと思うのだが。
0850132人目の素数さん
2007/08/25(土) 07:58:20せめて笑顔をふりまいてあげてください><
0851132人目の素数さん
2007/08/25(土) 08:17:560852132人目の素数さん
2007/08/25(土) 12:35:540853132人目の素数さん
2007/08/25(土) 14:26:140854132人目の素数さん
2007/08/25(土) 14:35:170855132人目の素数さん
2007/08/25(土) 14:50:080856132人目の素数さん
2007/08/25(土) 18:17:44Ctrl - だと右斜め下になってしまいます
解決策はありますか?
0857132人目の素数さん
2007/08/25(土) 19:41:54あなたの入力だと、数列などのindexと認識されます。Σの式とは完全に意味が異なります。
インストールしてすぐのdefault状態で出てくる
BasicInput.nb(標準入力)でΣを選び、入力したい項目に式を書いて下さい。
0858132人目の素数さん
2007/08/25(土) 20:43:120859132人目の素数さん
2007/08/25(土) 22:54:43MuPAD などがあります。
0860132人目の素数さん
2007/08/25(土) 23:10:21MuPADは売りもの。無料版はもはや存在しない。
0861132人目の素数さん
2007/08/25(土) 23:12:27回答ありがとうございます。
コマンドは同じなのでしょうか?
0862132人目の素数さん
2007/08/25(土) 23:51:00たたかう どうぐ
じゅもん にげる
0863132人目の素数さん
2007/08/26(日) 00:43:580864132人目の素数さん
2007/08/26(日) 01:06:09無くなったのか…
大学生時代にお世話になったよ・゚・(ノД`)・゚・。
0865132人目の素数さん
2007/08/26(日) 01:17:31ありがとうございます。パレットから選ぶしかないんですね・・・
ΣやΠ記号の下の文字を含めて
ショートカットだけで書くことはできないでしょうか・・・
さっき
苦労の末
ベクトルa→ を a ctrl 7 Esc rule Esc で書けることを知りました。
0866858
2007/08/26(日) 01:18:290867132人目の素数さん
2007/08/26(日) 01:53:33Esc sum Esc Ctrl+Plus n=1 Ctrl+% 100 Ctrl+Space n
な感じでできると思う。ver.6でしか確認してないけど。intやPiはsumを置き換えてね。
Ctrl+%とCtrl+Spaceは、汎用的なショートカットなので覚えておくと吉。
0868132人目の素数さん
2007/08/26(日) 01:58:11いまならMaximaかな。
yacasはさすがにだめか。オレはまだぎりぎり戦力として使ってるけど。
可換代数とか個別分野ではそれぞれ分野ごとによくできたソフトがあります。
と書いてるけど、可換代数と関連する代数幾何以外のソフトがよくできてるか知りません。
0869132人目の素数さん
2007/08/26(日) 02:16:05うおぉぉ〜すげ〜
ありがとうございました((ヾ(。・ω・)ノ☆゚+.ァリガトゥ♪
沢山のΣとΠが出てくる数式で・・・
パレットにもううんざりしてたところです
0870869
2007/08/26(日) 02:20:51Word2007で数式を書こう!という本が売られていました。
凄い凄いって著者は言うんだけど、パレットのオンパレード
・・・おまけにもちろん式展開などの計算は自分でして
1つ1つ打ち込みしないといけないから・・意味ないですよね・・・
0871132人目の素数さん
2007/08/26(日) 02:22:33計算結果込みで、数式がほとんどの文章だとMathematicaが向いてる。
もっとこだわる人はTeXだろうな。
0872132人目の素数さん
2007/08/26(日) 03:37:190873132人目の素数さん
2007/08/26(日) 03:47:07WordもTexもOOoMathも計算はできないなw
0875869
2007/08/26(日) 23:45:29f/@{0,1,2,3,4,5}]としたら、上手く出力されたんですが
いかんせん{f0,f1,f2,f3,f4,f5}というように、続けて横にズラズラ〜っと
でてきたので見難いです。
これを
f1
f2
f3
f4
f5
というようにタテに表示するにはどうしたらいいですか?
0876132人目の素数さん
2007/08/26(日) 23:49:240877132人目の素数さん
2007/08/26(日) 23:50:28MathematicaではComputer Modernのフォントって無いんでしょうか?
Ctrl + Shift + T でも 綺麗じゃないよね・・・
特に dx の d がいただけない・・・
0878132人目の素数さん
2007/08/26(日) 23:52:39それだと( )が出るんじゃないか?
0879132人目の素数さん
2007/08/27(月) 00:05:27Column
いろいろあるよ
0880132人目の素数さん
2007/09/01(土) 03:43:53ListLogPlot: y軸の目盛を対数にする
x軸が時間でy軸が対数の散布図を描きたい時はいったいどうすれば
いいんだ… orz。
WolframはDateListPlotの内部で使っているTicks関数を公開してくれたら良いと
思った。つーか俺が知らないだけ?
0881132人目の素数さん
2007/09/01(土) 13:07:480882132人目の素数さん
2007/09/01(土) 15:55:21式の本数と変数の数がnである連立方程式はどうすれば解くことができるでしょうか?
具体的にはx[0]〜x[n-1]を変数とした以下の連立方程式です。kは定数です。
x[0]-3 x[1]-2 x[2]-x[3]==4/k
x[0]+2 x[1]-3 x[2]-2 x[3]-x[4]==5/k
i>=3 && i<=n-3のとき
x[i-3]+2 x[i-2]+3 x[i-1]-3 x[i]-2 x[i+1]-x[i+2]==6/k・・・・(1)
x[n-5]+2 x[n-4]+3 x[n-3]-2 x[n-2]-x[n-1]==5/k
x[n-4]+2 x[n-3]+3 x[n-2]-x[n-1]==4/k
x[n-3]+2 x[n-2]+3 x[n-1]==3/k
です。x[i]を求めることはできるでしょうか。
上記の式をeqnlistに格納して(i>=3 && i<=n-3を除く)
RSolve[eqnlist,x[i],i]でいいような気がしますが、(1)式の差分方程式と他の境界条件の式が矛盾するためダメではないかと思っています。
もちろん、nにこちらで何かしらの定数を与えれば、Solveで解くことができますが、nに何も代入せずnの時の一般解を求めたいです。
よろしくお願いします。
0883132人目の素数さん
2007/09/01(土) 22:25:24自分で形を見るだけならそれで良いけど、人に見せる用のグラフなのでy軸には
あの対数表記のtickが欲しいのだ。今回は時間もないので結局Excelで描いた。
負けた気分(何にだ)。
0884132人目の素数さん
2007/09/02(日) 12:28:08ListLogPlotからY軸の目盛りを取り出す
yTicks = (Ticks /. AbsoluteOptions[ListLogPlot[data], Ticks])[[2]];
DateListPlotで対数加工したデータをプロットし、さっきの目盛を合成
DateListPlot[data /. {x_, y_} -> {x, Log[y]},
FrameTicks -> {Automatic, yTicks}]
これでめでたくプロットできる。
逆をやろうとしてDateListPlotからFrameTicksを取り出そうとしたら
AbsoluteOptionsがエラーを吐いて取り出せないのは、中で何か変なこと
してるんだろうか…。
0885132人目の素数さん
2007/09/09(日) 22:02:185.xまでは結果を丸めていたようだけど
6では打ち切りになっている
0886132人目の素数さん
2007/09/22(土) 19:41:370887132人目の素数さん
2007/09/23(日) 11:46:360888132人目の素数さん
2007/09/23(日) 18:56:570889132人目の素数さん
2007/09/24(月) 00:49:09PoissonDistribution[5]
という表記であっていますよね?
なぜか何もアウトプットされないのですが・・
0890132人目の素数さん
2007/09/24(月) 09:33:430891132人目の素数さん
2007/09/24(月) 17:50:16<< Statistics`DiscreteDistributions`は書いたか?
誰か教えて欲しい
混合分布に従う乱数を作ろうとしているんだが、Randomでは
「Real,Integer,あるいはComplexでなければなりません」と出て作れない。
それぞれの乱数を作って合体させたもので最尤推定量を求めると、
混合分布としての式とは数値がかなりずれている。
何かいい手があったら教えて欲しい。V5.2で作成中。
0892132人目の素数さん
2007/09/25(火) 23:57:52Ver.2 , Ver.4 , Ver. 6
特にVer.3はクソだった。Ver.1は機能不足。Ver.5はVer.4とそんなにかわらん
0893132人目の素数さん
2007/09/27(木) 23:55:38インタラクティブうんねんとは関係しない重い計算ごりごりやらせてるだけなんで。
0894132人目の素数さん
2007/09/28(金) 05:05:48数学系の人は前者を使う人が多そうですが
0895132人目の素数さん
2007/09/28(金) 08:33:54Matlabは行列計算システム。
前者は主に記号処理に重点を置いていて、
後者は主に数値計算に重点を置いている。
どっちでもできることもたくさんあるけども
うまく使いわけないと目的によっては非常に効率が悪くなる。
0896132人目の素数さん
2007/09/29(土) 23:15:01y=625.733*((2/5)*(1-x)*a2^2-(4/105)*(1+2*x)*a2^3);
があります。
このとき、a2は
a2->-((7*(-1+x))/(1+2*x))-s1;
です。
s1は
nが1のときは、s1=0
nが2以上のときは、s1=(92*n)/(y*10^6)
となります。
nを増やして、yが5.9に一番、近づいたときの
nを計算したいのですが、どのように計算すれば
良いでしょうか?
下記は、nを増やして、数値計算で求めたものです。
直接計算する方法を教えて下さい。
x=0.767476;
For[n=1,n<2*10^4,n++,
If[n==1,s1=0];
If[n>1,s1=(92*n)/(y1*10^6)];
y=625.733*((2/5)*(1-x)*a2^2-(4/105)*(1+2*x)*a2^3);
y1=y/.a2->-((7*(-1+x))/(1+2*x))-s1;
If[y1-5.9<0,Print["n=",n]];
If[y1-5.9<0,Print["y1=",y1]];
];
0897132人目の素数さん
2007/09/30(日) 21:09:07Maple でやってもらうと
問題の式を、x, y, n の関数として
f:= (x,y,n)-> y-(625.733*((2/5)*(1-x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^2-(4/105)*(1+2*x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^3));
で定義して、f( 0.767476,5.9,n)=0 の解 n を求めると n=13820.40245、59148.69709 となりました。
ということです。Mathematicaでは
x=0.767476;
z=Solve[y-(625.733*((2/5)*(1-x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^2-(4/105)*(1+2*x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^3))-5.9Š0,n];
Simplify[z]
では計算できませんが、どこが悪いのでしょうか?
0898132人目の素数さん
2007/10/01(月) 20:56:01失礼しました。下記で求まりました。
x=0.767476;
z=Solve[y-(625.733*((2/5)*(1-x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^2-(4/105)*(1+2*x)*(-((7*(-1+x))/(1+2*x))-(92*n)/(y*10^6))^3))==0,n];
y=5.9
Simplify[z]
0899132人目の素数さん
2007/10/02(火) 21:02:09やはり、解が一致しませんでした。
y=5.9
のときは、近い値が得られるのですが、
y=0.1
では答えが異なります。
如何でしょうか?
0900132人目の素数さん
2007/10/03(水) 09:57:34先に y=5.9 ってやってから Solve した方がよさそうですよ
>>898のだと丸め誤差の影響か虚数部がひどいことになりました
y=0.1のときは
{{n -> -347.016}, {n -> 651.842}, {n -> 742.06}}
ではだめなの?
0901132人目の素数さん
2007/10/03(水) 20:46:10お返事ありがとうございます。
>ではだめなの?
すいません。駄目なのです。
質問の本質は変更しませんが、
複雑な式は、分かりにくいだけなので止めて、
式を下記に変更したいと思います。
下記が、“0”に一番近づくのは
“7”付近なのですが、
このように数値計算をして求めるのではなく
式(たぶん微分方程式)を
解いて、下記と同様な解を
得るには、どのような式(微分方程式)
を作ればよいでしょうか?
a2=.;
For[n=1,n<2*10^1,n++,
If[n==1,s1=0];
If[n>1,s1=n/y];
y=3/2*a2^2+a2^3;
a2=2-s1;
f=5;
If[y-f<0,Print["n=",n]];
If[y-f<0,Print["y=",N[y]]];
If[y-f<0,Break[]];
];
もし、正しい式が得られましたら、fの値を
変更しても、正しい解が得られるはずです。
0902132人目の素数さん
2007/10/03(水) 22:46:350903132人目の素数さん
2007/10/04(木) 00:27:12お返事ありがとうございます。
>質問の本質はなんなのさ。
誤解を招く書き方をして申し訳ございません。
本来の目的は、896に書いたプログラム(数値計算)
から導かれる答えを、数値計算ではなく、式を作って
解いて、同じ答えを導きだしたかったのです。しかし
896に書いたプログラムは、やたら複雑でわかりにくいので
901のとおりプログラムを単純にしました。
私の希望はあくまでも、数値計算で求められる解を
数式を解いて、得ることです。
901で答えをご教示頂きましたら、それをそのまま
896に私が置き換えるつもりです。
0904132人目の素数さん
2007/10/04(木) 00:59:25整数a、bを係数とする2次方程式x^2+ax+b=0が有理数の解p/q(p、qは互いに素な整数)をもつならば、その解は必ず整数であることを証明せよ。
0905132人目の素数さん
2007/10/04(木) 08:44:53Mathematicaの問題じゃなくなったのでよそでやってくれ。
数学の初心者向けの板で馬鹿にされてくるといいよ。
0906132人目の素数さん
2007/10/04(木) 20:46:16お返事ありがとうございます。
>Mathematicaの問題じゃなくなったのでよそでやってくれ。
数%は、mathematicaの問題でもあるのでは?
>数学の初心者向けの板で馬鹿にされてくるといいよ。
まあ、そう言わずに、解っているのなら、せめて
答えだけでも教えて下さいよ!!直ぐに退場します。
この問題は、初心者じゃわからないです、、、
でも、かっこいいこと言って、本当は、答えが?では
0907132人目の素数さん
2007/10/05(金) 09:57:500908132人目の素数さん
2007/10/05(金) 11:09:51つまらなそうな世界だな。
0909132人目の素数さん
2007/10/05(金) 13:00:10パスワードの取得がやっかい。
Webに詳細な個人情報とパスワードが必要な理由ってのを打ち込んだらきたパスワードが4日間だけ有効って代物。
まともなパスワードほしければこれを見ろってんでみたらPDFファイル。
それ印刷して手書きで個人情報書き込んでパスワードがほしい理由を英語で説明して
アメリカへファックスし、approvalが降りるの待てときたもんだ。
何十万も払った上、事務手続きまで要求されるとは。
0910132人目の素数さん
2007/10/05(金) 14:47:13俺もPremierだけど、「あんたはPremierだからこっからダウンロードしてくれ。
ちなみにパスワードはこれ」とか言うメールが来て、別に面倒はなかったけど。
0911132人目の素数さん
2007/10/05(金) 15:22:38お前はもうダウンロードしたからだめと表示される。途中で切れててもね。
連絡して再認可が下るのを待てばいいんだけど、また失敗すると面倒なんで
今回はあえてCD送付を選んだわけ。そしたらこちらは事務仕事が強制だったw
ところで起動して眺めたら、デバッガーがついたみたいだけどこれって前からあったっけ?
もし使いやすいデバッガーなら嬉しい驚きなんで、とりあえず簡単なプログラム
書いて試したけど、やはりマニュアル読まないとだめみたい。
もう使ってみた人いたら感想教えて。
0912132人目の素数さん
2007/10/05(金) 20:31:51>つまらなそうな世界だな。
皆様、超多忙にも関わらず、結局結論は出ず、ただ時間の無駄使いでした。
結論:お騒がせしました事をお詫び致します。
0913132人目の素数さん
2007/10/05(金) 21:05:450914132人目の素数さん
2007/10/09(火) 00:01:52ありがとう。
0915あ
2007/10/14(日) 00:52:100916132人目の素数さん
2007/10/14(日) 06:36:15例
Integrate[x^2 y^2 Boole[x^2/a^2 + y^2/b^2 <= 1],{x,-Infinity,Infinity},{y,-Infinity,Infinity},Assumptions->{a>0,b>0}]
0918132人目の素数さん
2007/10/15(月) 01:21:30mathematicaで例えばT[x]とR[x]を畳み込み積分させるにはどうすればいいですか?
ググってもわかりまへんでした(#+_+)
0919132人目の素数さん
2007/10/16(火) 22:15:130920132人目の素数さん
2007/10/16(火) 23:53:21問題なくできる。まあ同じ設定ファイルを読もうとするから書き方に気をつける
必要あるけど
レス数が900を超えています。1000を超えると表示できなくなるよ。