トップページ⇒math
1001コメント292KB

〓 Mathematica 3 〓

■ このスレッドは過去ログ倉庫に格納されています
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/
0183132人目の素数さん2006/01/31(火) 01:12:30
>>182
こんな感じ?
MapThread[Max, {{{1, 2, 1}, {2, 3, 1}}, {{3, 0, 1}, {2, 1, 5}}}, 2]
0184132人目の素数さん2006/01/31(火) 03:11:14
>>175
お宅はヴァカ?
鍵genが出回っていること知らないの?
この鍵genで生成されるpassだと英語しか使えないけど、普通は不自由しないよ。
01851822006/01/31(火) 23:24:20
>> 183
ありがとうございます。行列をこねくりまさなきゃいけないかと覚悟してましたが、助かりました。
0186132人目の素数さん2006/02/01(水) 22:22:44
>>184
ありがとうございました。
通報させていただきました。
0187132人目の素数さん2006/02/02(木) 15:21:42
一辺がaの正四面体OABCがあります。
OA,OB,OC上に点P,Q,Rをとったところ、
四面体OPQRの体積が四面体OABCの1/3になりました。
このような三角形PQRの周および内部が通過する領域の体積を求めてください。

今↑の問題を解くためにMathematicaを使おうと考えています(普通に解けるぞ!という方いらしましたらご解答お願いしますm(_ _)m)。

上の方の相似比1/3(体積比1/27)の領域は通らないことはわかっています。
それで、下の方の「ぽっこりした領域」も通らないことがわかっています。
で、この「ぽっこりした領域」の体積を求めることができれば、
(四面体の体積)−(相似比1/3の領域)−(ぽっこりした領域)=求める領域
なので、この「ぽっこりした領域」の体積を求めようとしています。

これが相当むずかしくて、なかなか求められませんので、なにか知見を得るために
Mathematicaにこのぽっこり領域を描画させていろいろ切断してみたいなぁと考えています。

どのようにしたら良いか分かるかたがいらっしゃいましたら、教えていただけると幸いです。
0188132人目の素数さん2006/02/04(土) 14:09:22
5.2になりました
0189132人目の素数さん2006/02/07(火) 09:45:02
ArcTan(x)=0,初期値X=1から始まる方程式を「再帰的に」解くプログラムを書け。
という問題です。使用するソフトはマスマティカです。
級数展開などいろいろすると
ArcTan(x)=∫1/(1+t^2)dt
=(Σ (2k - 2)! / (2k - 1)!){x^(2k - 1) / (1 + x^2)^k}
となる所まで辿り着きました。
階乗が出ているからここから再帰的に解くプログラムを何とか作れるような予感はするのですが、
Σ (2k - 2)! / (2k - 1)!に関して、
f[k]をf[k-1]の式で表す、つまり漸化式のような形にして解けばいいのでしょうか。
再帰方程式自体はコマンドRsolveで解けるらしいのですが

0190132人目の素数さん2006/02/09(木) 23:56:25
アップグレードの料金が高すぎると思う。
MATHEMATICAもMAPLEもなんでこんなに高い値段で売るのだろうか?
そのような高い値段のソフトを学校はサイトライセンスだと安く買える
とかいう理由なのか導入して、それで教育したり教材を作るが、
それは結局価格の高いソフトの広告宣伝塔として学校がソフト会社の
戦略の一翼を担っていることに他ならない。そういうところでほぼ
人生の一時期にでも強制して使用させることで、依存性と習慣性を
高めて、一生ライセンス料とアップグレード料金を吸い上げ続ける
為のイニシエーションの機関となっているのも同然だ。
マイクロソフトのOSやワープロ、表計算、プレゼンテーションツール
に関しても同様だ。学校はそのような外国製のソフトで高額のものに
学生をはめ込み、海外の企業に利益を吸い上げさせる構造を支援している。
 かつてイギリスがインドを植民地として、繊維製品を作ることを禁じ、
イギリス製の綿製品を買わせたり、紅茶に入れる砂糖に重税を課して
いたように。官公庁や特に学校が海外製品でブラックボックスのものを
使うことをほぼ強要するのは、自国民を植民地化させ、自国のソフト産業
を壊滅させるために作業しているようなものである。
0191132人目の素数さん2006/02/11(土) 15:30:28
フリーのソフトつくって下さい。
0192132人目の素数さん2006/02/12(日) 01:49:43
つmupad,maxima,risa/asir
0193132人目の素数さん2006/02/14(火) 20:26:09
質問です。以下のような簡単な定積分の計算をさせたところ、
ものすごい長い、わけわかんない答えが出ました。
何ででしょうか?
Integrate[(beta*t + c)^(1/h), {t, i, j}]
0194132人目の素数さん2006/02/14(火) 22:50:44
>>193
Assumptionsを使って条件を限定すれば、それなりに短い答えが
返される。
0195132人目の素数さん2006/02/15(水) 23:25:35
質問です。MathematicaをインストールしていないPCに
Mathematicaからの文字をワードなどにコピペしたものを表示させたいのですが
あいにくフォントがなく文字化けします。
以前はマサチューセッツ工科大学のどこかにフォントリンクが張ってあったようですが
見つけれず… フォントのありかをご存知の方いますか?
0196◇ゾンビ@Linuxザウルス : ◆mzpjBMNq3U 2006/02/16(木) 11:20:18
>>195
WRIのサイトからMathReaderをダウンロードして
インストールしてみてはいかがでしょう
画面表示用のFontが入っていたと思います

ただ、MicrosoftWord上で表示できるにしても
うまく文字間隔や1などを再現できるかは
わかりません
0197132人目の素数さん2006/02/16(木) 18:13:07
wordに貼り付けるときに図で貼ってみれば?
0198132人目の素数さん2006/02/16(木) 23:53:54
LinuxのPCにインストールした。そのマシン名を仮にPC1としよう。
自宅のLinuxマシン(そのマシン名を仮にPC2としよう)から
SSHでPC1にログインしてPC2上のX画面でそのGUI版のMathematicaをPC1
で起動してやると、PC1上のXの画面上に、Mathematica用のフォントが
ありませんとかいうようなメッセージが出て、数式の表示もおかしいし
いろいろとまともに表示されない。これはどうしたら回避できるのだろうか?
Linuxは両方のマシンで同じFedora Core for AMD64 だ。
0199132人目の素数さん2006/02/17(金) 10:17:49
FAQだな.PC2にフォントをインストールしろ.
0200132人目の素数さん2006/02/18(土) 08:44:56
age
0201132人目の素数さん2006/02/18(土) 09:10:28
>>198
フォントサーバでいけるのでは?
0202132人目の素数さん2006/02/22(水) 06:48:56
>>190
フリーですが、何か?

>>198
両方にインスト。 コレ最強。 鍵gen使えば何台でも桶。
0203132人目の素数さん2006/02/22(水) 16:25:57
インテルマック版来たね。
人柱よろ。
0204132人目の素数さん2006/02/23(木) 10:28:26
しこ
0205132人目の素数さん2006/02/23(木) 16:57:47
mathematica で、XYデータセット(非対称)
における分散(半値全幅でもいい)とか求めたい
のですが、descriptive statisticsをつかえば
良いのはわかります。
具体にどうしたらいいのでしょうか?
0206132人目の素数さん2006/02/27(月) 10:53:29
亀の動いたあとが、樹のフラクタル図形になるような、亀の始めの動きtと
置き換えルールrule6を定義し、作図せよ。・・・さっぱりわかりません。
ヒント:枝の部分をxとおき、この部分を置き換えるルールを考えよ。だそうですが。
あと,
move[z]:=(p={{0,0},{0,1}};p[[1]])
move[f]:=(p[[1]]=p[[1]]+p[[2]];p[[1]])
kaiten[th_]:=p[[2]]={{Cos[th],Sin[th]},{-Sin[th],Cos[th]}}.p[[2]]
move[r]:=(kaiten[Pi/3];p[[1]])
move[l]:=(kaiten[-Pi/3];p[[1]])
が、前の問題で定義されています。
0207132人目の素数さん2006/03/02(木) 19:11:52
320
0208132人目の素数さん2006/03/04(土) 15:26:25
質問なのですが、

y=ax+.....

のグラフをx軸を縦に、y軸を横に表示したいのですが
どうすればいいでしょうか?
0209132人目の素数さん2006/03/05(日) 13:04:55
>>208
頭を90度かたむける
0210132人目の素数さん2006/03/16(木) 22:50:48
実験データをy=a/x+bの式にフィッティングしたいのですが、
データの中に{x.y}={0、100}を含んでいるので
無限式が出ると言われて上手くいきません。
どのようにすればいいのでしょうか。
0211132人目の素数さん2006/03/16(木) 23:20:23
>>210
0のでるデータを外す。
普通、実数は0で割れないでしょ。
02122102006/03/17(金) 20:11:16
>>211
ありがとうございます。やはりそうするしかないのでしょうか。
データ数が少ない、かつx=0付近が重要なのでそれは避けたいのですが。
0213132人目の素数さん2006/03/23(木) 09:38:11
>>210
自分で(フィッテングの)評価関数書いて、評価関数の最適化だけやらせるとか。
自分で書けば、「x=0付近が重要」と云うのも、それなりに盛り込めるでしょ?
0214132人目の素数さん2006/03/26(日) 15:03:14
0215132人目の素数さん2006/04/03(月) 21:43:37
そろそろ、新しいMathematicaでるのかな?
http://www.wolfram.co.jp/news/events/JMC2006/
0216132人目の素数さん2006/04/05(水) 17:20:18
ある数列が与えられていて、それがベータ関数の値であるときに、
元のベータ関数のパラメータはどうやって求めればいいんでしょうか。
V5.2を使っています。
0217132人目の素数さん2006/04/14(金) 21:03:53
wktk
0218132人目の素数さん2006/04/15(土) 12:01:53
age
0219132人目の素数さん2006/04/16(日) 01:39:42
0220132人目の素数さん2006/04/17(月) 01:42:18
Ver6マダー
0221132人目の素数さん2006/04/23(日) 22:59:53
マーダダヨー
0222132人目の素数さん2006/04/29(土) 22:05:27
マスマスティかで不等式を2DPlotで描くにはどうすればいいの?
領域の色づけがヘルプを見ても分からない。
0223132人目の素数さん2006/05/01(月) 10:41:08
Mathematicaでアイパターン描こうとしているのですがなかなかうまくいきません
やり方としてはNRZのランダムパルスを発生させそれを離散フーリエ変換してナイキストフィルタのリストと乗算
それを逆フーリエ変換して同一平面上にプロットすれば良いと思うのですがほかにやり方あったら教えてください

>>222
filledplot
0224132人目の素数さん2006/05/10(水) 12:56:10
√(sin 2x+2cos 5x)のグラフ書くために初めてmathematica使ったんですけど、
グラフがx軸に接しないのはなぜでしょうか?
0225132人目の素数さん2006/05/10(水) 22:20:26
>>224
接するグラフではないのだろう
0226132人目の素数さん2006/05/11(木) 16:36:46
FindRootの結果をグラフにしたいのですが、なんかスマートな方法ないですか?

普通にFindRootすると戻り値は
{hoge -> 51.6438}
みたいな文字列になってしまってるので困ってます。
0227132人目の素数さん2006/05/11(木) 21:16:10
>>226
ReplaceAllを使う
x /. FindRoot[x^2 == 2, {x, 1}]
みたいな感じで
0228132人目の素数さん2006/05/13(土) 12:15:14
>>224
Solve[(Sin[2x] + 2 Cos[5x])^0.5 == 0, x]
の結果、実数の範囲では
{x -> -1.5707963267948966`} {x -> 1.5707963267948966`}

Plotでは限界がある計算をさせてるみたい。
0229132人目の素数さん2006/05/13(土) 21:21:13
次のMathematicaの情報知ってる人いますか?
4.2から5のときのReduceはすごかったけどそれに匹敵するようなグレードアップ
があるといいなあ。
0230132人目の素数さん2006/05/14(日) 23:33:37
二つのベクトルから行列を作る方法を教えてください。
v . w^t のような感じで、v.Transpose[w] だとできません。
0231132人目の素数さん2006/05/15(月) 01:30:13
>>230
Outer[Times, w, v]
0232132人目の素数さん2006/05/18(木) 17:30:12
x^2(x-3)(2x+3)と2点で接する半径1.5の接点を求めたいのですがどうすればよろしいですか?
0233132人目の素数さん2006/05/20(土) 00:31:01
このソフトの最新バージョンナンバーはいくつですか?
0234132人目の素数さん2006/05/20(土) 00:34:29
5.2
0235132人目の素数さん2006/05/22(月) 16:22:35
ある曲線に、2点で接する円を書きたい時はどうすればいいんですか?
0236132人目の素数さん2006/05/24(水) 19:52:20
うちの教授が

「マテマチカ」

と呼ぶアプリはこれのことでしょうか?
0237132人目の素数さん2006/05/25(木) 16:19:16
ループを回して計算しているその途中経過をグラフィカルに表示するにはどうすればいいですか
計算が進むにつれ折れ線グラフがにゅーんと延びていくとかできればベストです。
0238132人目の素数さん2006/05/25(木) 19:41:40
>>236
その発音も間違いじゃない。
0239132人目の素数さん2006/05/26(金) 08:58:14
本屋さんの名前と混同してるんだろう
0240132人目の素数さん2006/05/27(土) 06:35:34
Mathematicaは線形連立方程式を数値的にではなく
代数的に解くことはできますか?
つまり
x=a+b
y=a-b
を解いて
a=(x+y)/2
b=(x-y)/2
のような答を、a,bに具体的な値を入れなくても
返してくれるのでしょうか?
0241132人目の素数さん2006/05/27(土) 06:40:20
できる
0242132人目の素数さん2006/05/27(土) 12:10:53
そろそろ、ver6.0が出るのかな?
0243132人目の素数さん2006/05/29(月) 23:55:30
n==0のとき、 f[x_,n_]:=1
それ以外のとき、f[x_,n_]:=x^n
となるようにf[x_,n_]を定義したいのですが、
どうしたらいいのでしょうか?
02442432006/05/30(火) 00:03:13
ちょっと調べたら見つかった
f[x_, n_] := If[n == 0, 1, x^n]
スレ汚しスマソ
0245132人目の素数さん2006/05/31(水) 14:40:52
出力がリストになるような関数をコンパイルするにはどうすればよいのでしょうか。
下記のようにうまくいきません


以下失敗例

In[1]:= test[X_, Y_] := {X + Y, X - Y}
In[2]:= comp = Compile[{X, Y}, test[X, Y]]
Out[2]= CompiledFunction[{X, Y}, ttest[X, Y], "-CompiledCode-"]
In[3]:= comp[1, 2]
CompiledFunction::cfse: コンパイルされた式(3.,-1.)はmachine-size real numberであるべきです.
CompiledFunction::cfex: "式2の外部評価においてエラーが発生しました.コンパイルなしで評価を続行します."
Out[3]={3, -1}
0246132人目の素数さん2006/06/06(火) 21:04:38
>>242
もう機能満載で大きな変化はないんとちゃう? 
0247132人目の素数さん2006/06/07(水) 21:13:56
Mapleみたいに、外観が変わるんじゃない?
0248132人目の素数さん2006/06/10(土) 08:40:56
mapleはだんだん教育ツール化している気がする。
mathematicaはどうだろう?
0249132人目の素数さん2006/06/11(日) 15:07:04
計算ツールになってる。
0250132人目の素数さん2006/06/15(木) 01:51:12
チュートリアルの例文そのまま
result = NDSolve[{
y''[t] + 0.1 y'[t] == -Sin[y[t]], y[0] == 1, y'[0] == 0}, y, {t, 0, 50}]
をコピペしただけなのに
NDSolve::deqnの最初の引数にはFalseではなく,方程式か方程式のリストが必要です.\"\>"}]\)
とかエラーが出やがるんですがどうしてでしょうか?
0251羽村2006/06/15(木) 16:20:27
>>250
Clear[y];
して再チャレンジ
0252UNI2006/06/16(金) 23:17:15
Mathematicaでアポロニウスの円を書きたいんですけど、どうすればよいでしょうか?
0253132人目の素数さん2006/06/17(土) 11:53:30
>>252

Print["アポロニウスの円"]
0254132人目の素数さん2006/06/17(土) 13:04:14
>>252
定義通りに式を書いて、それをPlotすればいいのでは?
0255132人目の素数さん2006/06/17(土) 20:57:27
噂に聞くMathematicaそのうち買おうと思ってますが、ver.6そろそろ出るそうでないですか。
メジャー ver upによってエンバグされたりすると悲しいので、今のうちに買おうかなと思ってるんですが、
何分はじめてなもんで、この判断が妥当かどうか微妙です。

皆様の過去のver up体験談(4 ==> 5で酷い目にあった or 全く問題なかった)などありましたら
聞かせて頂けないでしょうか
0256132人目の素数さん2006/06/17(土) 21:18:10
>>255
質問とはそれるけれど、
学生版以外は、Premier Serviceが無料で付くので、
ver.6に無料で変えれるよ。
0257132人目の素数さん2006/06/17(土) 21:52:50
>>255

ver4.0==>ver5.0(Win)で酷い目にあったこと
・立ち上げたときにピロロロンと鳴らなくなった
・境界条件に文字を含んだ簡単な微分方程式が解けなくなった
などなど
0258132人目の素数さん2006/06/17(土) 22:31:22
>>257
>立ち上げたときにピロロロンと鳴らなくなった
これは、環境設定のところで音を鳴らすことができる。

ver5.0はバグが多かったなぁ。
02592552006/06/17(土) 23:52:01
お返事ありがとうございました。

>>256
Academic(~20万)で買おうと思っとるんですが、どうもPremier Serviceは別売(~10万)のようです。
(今なら2万円程度おトクのようですが。。。)

>>257,258
むむ、やはり不具合などありましたか。

アドバイスありがとうございます。(比較的)枯れてる現バージョン購入の方向でいってみます
0260UNI2006/06/19(月) 22:46:12
お願いします!!Mathematicaを使って書けっていう問題が出ました・・・
問1 2点からの距離の比が一定な点の軌跡(アポロニウスの円)のグラフをかけ。
問2 1点と1直線から等距離にある点の軌跡(放物線)のグラフを書け。

コマンド教えてください!!
お願いします☆
0261GiantLeaves ◆6fN.Sojv5w 2006/06/20(火) 08:13:45
talk:>>260 軌跡がどうなるかあらかじめ計算して入力すればいいのではないのか?
0262UNI2006/06/20(火) 15:40:06
Mathematica上で、計算していきたいんです!!
お願いします☆
0263132人目の素数さん2006/06/20(火) 20:08:33
\!\(Appolonius[m_, n_] := Module[{a, b, r}, \[IndentingNewLine]a =
1 + n\^2\/\(m\^2 - n\^2\); \[IndentingNewLine]b =
0; \[IndentingNewLine]r =
1 + n\^2\/\(m\^2 - n\^2\) - m\/\(m + n\); \[IndentingNewLine]Show[
Graphics[{Point[{0, 0}], Point[{1, 0}], Circle[{a, b}, r], Point[{a,
b}]}, AspectRatio -> Automatic]]]\)
0264UNI2006/06/25(日) 11:54:58
Mathematicaで、x軸に垂直な直線を書きたいんですけど、どうすればいいですか?!
0265132人目の素数さん2006/06/25(日) 20:57:35
お前、それは中学校の数学で習う範囲だろ('A`)
小学生だってx=1とかx=0は知ってる人は知ってる
0266132人目の素数さん2006/06/25(日) 23:48:54
<<264

x軸に垂直な直線、例えば、x=1 を書きたいとき、

<<Graphics`

ImplicitPlot[x==1,{x,-5,5},{y,-5,5},AspectRatio -> Automatic]
0267132人目の素数さん2006/06/26(月) 17:29:00
20行1列のリストをX軸に、
もう一つの20行1列のリストをY軸に、といったプロットって出来ますか?
簡単に出来そうに見えて出来ない・・・
0268132人目の素数さん2006/06/27(火) 02:18:25
>>267
x={{x[1]},...,{x[20]}};
y={{y[1]},...,{y[20]}};
ListPlot[Flatten[Transpose[x_List]],Flatten[Transpose[y_List]]]
でいいんじゃない?
02692672006/06/27(火) 12:48:02
>>268さん
レスありがとうございます。
やってみたのですが、うまくいきませんでした・・・。
前の書き込みでは説明不足だったのですが、
x={{x[1]},...,{x[20]}};
y={{y[1]},...,{y[20]}};
のところは、私がやろうとしていたのは、
Excelのデータから抜き出したデータの羅列、例えば
x={{1},{2},{3}};
y={{3},{4},{5}};
のような感じで、番号を振ることができないのでうまくいかないのかなと思います。
ヘルプから、これらのようなxとyを、
n={{1,3}{2,4}{3,5}};
といった行列に直せば良いことは分かっていたのですが、この変換がクセモノで。。

勉強不足と言われればそれまでですが、
妙案がございましたら、レスしていただけると幸いです。
02702682006/06/27(火) 13:00:27
>>269
そういうnを得たいんだったら、
n=Transpose[{Flatten[x],Flatten[y]}]
でどう?
02712672006/06/27(火) 13:28:27
>>268さん
素早いレスありがとうございます!
希望の通りの行列ができました!!
本当に、ありがとうございます!!

もっと精進していこうと思います。
0272132人目の素数さん2006/06/28(水) 03:08:19
アルキメデスのらせんを書く方法を教えて下さい。
もうさっぱりです。
0273132人目の素数さん2006/06/28(水) 08:23:47
極座標で
r=nθ
を表現すればいい。
0274132人目の素数さん2006/06/28(水) 21:09:38
すいません、極座標の表現の仕方がわかりません。

x:=Cos[θ]
y:=Sin[θ]

とおいて、置換で書いた見たんですが直線になってしまいます。
02752682006/06/29(木) 03:02:00
>>272
<<Graphics`Graphics`
PolarPlot[t,{t,0,2 Pi}];
0276132人目の素数さん2006/06/29(木) 19:31:46
>>275
書き方がわかりました。
ありがとうございました。
0277132人目の素数さん2006/06/30(金) 12:55:28
等高線に対しての流線の描き方がわかりません。
よろしくお願いします。
0278132人目の素数さん2006/07/04(火) 12:10:11
正の整数 N に対して,√N を根号の中ができるだけ小さな整数に
なるように変形して a√b の形に直すとき,この整数 a と b の
求め方を教えて下さい.よろしくお願いします.
0279132人目の素数さん2006/07/04(火) 23:50:17
>>277,278
かんがれ!
0280132人目の素数さん2006/07/06(木) 14:23:48
>>278
root[27]で実行すれば、3root[3]とかなるが、こういうことじゃないのか?
0281132人目の素数さん2006/07/12(水) 22:08:19
あるセルを評価したときに結果がとても膨大になってしまうことってありますよね。
そういうときヘルプの中のセルのように、見出しだけ表示しておいて、
横のところをダブルクリックすると、中身が出てくるようにしたいのですが、
どなたかやり方をご存知の方いらっしゃいますでしょうか。
どうかよろしくお願いいたしますm(_ _)m
0282132人目の素数さん2006/07/13(木) 23:07:49
>>280
その 3 と 3 とを取り出すにはどうしたらいいですか?
■ このスレッドは過去ログ倉庫に格納されています