〓 Mathematica 四 〓

　　　　　　　　_....._{{ 〃
　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-<
　 　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）)
　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ.
　　　',　',　|Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　 ﾞ',
　　　 ヽ､', l:!Kノ}. 　 　 f:_.)iﾞi: ﾘ !　l　ル　
　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉ ﾞ ,ﾘ l ﾘ'´ 　 　　
.　　　　 ',|!!､ 　 　r‐┐ 　 ｀ ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣ 　
　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　 ＜　 Mathematica に関する話題はここに書いてね！
　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　 　 　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､ 　
　　/＼＼　　　　',　　 }　　 //｀ヽ

過去スレ
>>2

〓Mathematica〓
http://science2.2ch.net/test/read.cgi/math/985023298/
〓 Mathematica 2 〓
http://science3.2ch.net/test/read.cgi/math/1078534285/
〓 Mathematica 3 〓
http://science6.2ch.net/test/read.cgi/math/1121413040/

3et

　　 　 　 　 　 ／　　　　 ／￣￣｀ヽ}/´￣￣＼　　　　 ＼
　　　　　　　,ｲ　　　　 ／_　　 -‐――‐‐- 　_　 ヽ 　 　 　 ヽ
.　 　 　 　 / 　　　　,ｲ／ 　 　 　 　 /　　　 　 ｀　 ﾍ　　　 　 ﾍ
　　 　 　 / ./　 　 /　 　 　 /　 /　/　　 　 ∧　 　 ∧　　　　 ﾊ
.　 　 　 / ./　 　 /　 　 　 /　 /　/　　　　　 ﾊ　　　 ﾊ　　　　　l
　　　　 i ./　 　 /　/　　 /　/,'　 {　　 　 　 　 l　∧　 ﾊ　　　　｜
　　　　 |/ / 　 ,'.　 !　　/　/. ! 　ﾊ.　 lﾄ､ ,ﾊ　｜　ﾊ　　ｌ　　 　 ｜
　　　　 |　i　　 i　　|　 /| .ﾊ　|　 | {　 |　i　|l　 ﾄ､ 　}l　｜｜　　 |
　　　　 |　|　|　| 　 l　{Hl {ﾅT爪 |　V |　| 爿　辷} 八　ﾘ　l　　　l
　　　　 |　|　| ,ﾊ.　 ﾄ　|ｨfﾃ示kヽ|　 ヽl　}/ 」/rzｯ/_｜/　,'　 　 l
　 　 　 ヽ ﾄ､ﾄ､ ∧ l! /|l{､__,ﾊ|　　 　 　 　 ´/{ 　 }心/}. /　 　 ,'
.　　　　　 　 〉､　∧l　 l刈刋ﾘ　　　　　　　 {州爪ﾘ }/,'/　　　/
　　　　　　　l　 ＼ ヽ.　ヽ--''　　　　　　　　ヽrzン　 /　　　 /
　　　　　　　|｜ 八 {　 Xxx　　　　　!　　 　 　 xxX / 　 　 /
　　　　 　　 ヽﾄ､ 　 ＼　　　 　 　 ｰ 　‐　 　 　 　 /./　　 /
　　 　 　 　 　 ヽ＼　 ヽ、 　 　 　 　 　 　 　 　 ／,.'　　 ｲ|
　 　 　 　 　 　 | 　 ヽﾊ、＞ ､_ 　　　　　　　,. ィ／　 ／/|i!
　 　 　 　 　 　 |　　　　 ｀¨ヽk}＞ ､ __,. r＜{r-く- ''"　　 ﾘ 　>>1　スレ立て乙です
　 　 　 　 　 　 l　　 　 　 ／￣＼　　 　 ／￣ヽ　　　　 ,'
　　　　　　　 　 ｌ　　　 ／　　　 　 `ｰ‐''´ 　 　 ∧　 　 /
　　　　　　　　　ヽ _,r〈 　 　 　 　 /⌒V　　　　　ﾊ-､_/
　　　 　 　 　 ,.r＜＼　|　　　 　 / l￣l V　 　 　 / ／＼
　　　　 ,. -＜.　　　 　 | 　 　 ／￣.　.￣ﾄ､　　 /　 　 　 ｀丶、
　　　／　　　　　　　　 ｌト-＜ └┐ .┌┘ ヽンl　　　 　 　 　 丶、
　　/{　　　　　　　　 　 l　 　 |l　 └‐┘ 　||　　l　　　　 　 　 　 　 ヽ
.　/ .|　　　　　　　　　　l　 　 ||　　　 　 　 || 　 l　　　　　 　 　 　 　 ﾄ､

欲しいけど高すぎて買えない…

簡単な質問かもしれないけど、教えてください。
y''=3とかの微分方程式って
DSolve[y''[x]==3,y[x],x]
て、入力してもとけないんですけど、どうやったらとけますか？
すいませんが、教えてください。

>>6
解けるけど、、、。
{{y[x] -> (3 x^2)/2 + C[1] + x C[2]}}

まんまコピペして解けるな

質問なんだけど

「
:関数取得( := )(Input入力)
:変数取得( = )(Input入力)
:関数定義( := )(Module定義)
:関数定義( := )(Module定義)

:定義関数呼び出してのメイン演算、条件文込み
」

てのが一行ずつのインタプリンタでは上手く実行されてるとして、これを
通常のプログラムみたいに何かしらの一個の読み込みコマンドを入れたら
全部包括的にやってくれるようにまとめ直したいんだけど、
どうすればいいんですか？
あと、MessageDialogue、Input、ChoidceDialogの実行に明確な順列をつけるには
どうすればいいんですか？

抽象的な質問で非常に申し訳ないんだけど、参照するコマンドとかだけでもわかればありがたい。

str = {12371, 12435, 12394, 26178, 38291, 12395, 12289, 12371, 12435,
12394, 12367, 12384, 12425, 12394, 12356, 12418, 12398, 23455,
34892, 12375, 12390, 12288, 12393, 12358, 12379, 20170, 24180,
12418, 19968, 20154, 12398, 12463, 12522, 12473, 12510, 12473,
20057, 65367, 65367, 65367};
FromCharacterCode[str]

こんな時間に、こん
なくだらないもの実
行して　どうせ今年
も一人のクリスマス
乙ｗｗｗ

str4 = {104, 116, 116, 112, 58, 47, 47, 119, 119, 119, 46, 102, 97,
109, 105, 116, 115, 117, 46, 99, 111, 109, 47, 103, 97, 109, 101,
47, 99, 111, 109, 105, 110, 103, 47, 50, 48, 48, 55, 47, 48, 55,
47, 48, 51, 47, 104, 45, 49, 48, 52, 95, 55, 52, 51, 56, 52, 95,
114, 52, 53, 46, 106, 112, 103, 46, 106, 112, 103};
Import[FromCharacterCode[str4]]

FindRootで
{x -> u ,y -> v}
と求めた後、u,vをそれぞれx,yの値とするにはどうしたらいい?
いまはいちいち手入力で
x=u;
y=v;
とか打ってるんだけど、面倒で。
あほーな質問でｽﾏｿ。

いつも
{x0, y0} = {x, y} /. {x -> u, y -> v}
という感じで解をx0,y0に代入してる

x,yに解を代入しないのは、もう一度x,yで方程式記述するとき
すでに値が入っていると思うように動かないので、それを避けるため

１番目の解なら
{x, y} = {x, y} /. %[[1]]

さらに>>14方式なら
{x0, y0} = {x, y} /. %[[1]]

でおｋ

>>15
わずらわしさが半減しました。どうもです。

Mathematica Playerがダウンロードできないんですがどんな理由が考えられるでしょうか
ttp://www.wolfram.com/products/player/download.cgi
ttp://www.wolfram.co.jp/products/player/download.cgi
を開いても真っ白のままで何も起こらないです
Windows XP+IE6とVista+IE7で両方同じ現象。

>>17
XP+IE7で確認しました。文字コードの影響みたい

右クリック -> エンコード -> その他 -> Unicode(UTF-8)

おおー登録画面出ました
せっかくなので登録画面のご意見欄にこの問題のこと書いておきます

社会人なのですが
業務の関係でmathematicaを購入することになりました
どこでどのようにして購入すればよいのでしょうか？
やはりウルフラムから購入するのでしょうか・・・

test

スレタイ合成のためにage

>>20
はい。

大きなソフト売り場を持つ店舗なら、店頭でも買えますよ。

基本的にMatlabとS-Plusがメインだけど、Mathematicaも便利な道具として手放せない。
積分や面倒な計算も結構すぐやってくれる。記号計算が出来る使いやすいソフトは
やはりMathematicaだと思う。昔はMapleの方が性能が良いと言われていたけど今は
どうなのかな？Mapleは使い慣れないせいか、ほとんど使ってない。

>>1
新スレおめでとうございます

　 　, 　 ＿ ﾉ）
　γ∞γ~　 ＼
　 |　 /　从从) ）
　ヽ　| |　l　 l |〃/ＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＮ＼
　 ｀从ﾊ~ ワノ∩　ハニャウェ━━━━━━━━━━━━ブ　>εε=ヽ( `Д´)ﾉ ｳﾜｧｧｧｧﾝ
　　　　　　 ⊃　.ＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＶＮ／

　 　, 　 ＿ ﾉ）
　γ∞γ~　 ＼
　 |　 /　从从) ）　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　人人人人人
　ヽ　| |　l　 l |〃　ﾆｺﾆｺ 　　　　　　　　　　　　　　　　　　　　　　　　　≪ヽ( `Д´)ﾉ　≫　ﾋﾞﾘﾋﾞﾘ
　 ｀从ﾊ~ ワノ）　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　^Ｙ^Ｙ^Ｙ^Ｙ^Ｙ

　 　, 　 ＿ ﾉ）
　γ∞γ~　 ＼
　 |　 /　从从) ）
　ヽ　| |　l　 l |〃　ﾊﾆｬｰﾝ！　　　　　　　　　　　　　　　　　　　　　　　 　ヽ(・∀・ )ﾉ ﾊﾆｬｰﾝ！
　 ｀从ﾊ~ ワノ∩

>>23
わざわざウルフラムから買わなくても
日本に代理店とかあったんじゃなかったっけ？

ダウンロード版の方が安い。

OSをVistaからXPに戻そうと思っているのですが、
Mathematicaを再インストールするときは、パスワードを再発行してもらう必要はありますか？
MathIDって変わってしまうんでしょうか？

>>28
MathIDってことはスチューデント版だと思うんですが、
再インストールすると確かに変わりますので再発行が
必要です。

>>29
ありがとうございます。
再発行してもらいます。

リモートカーネルがうまくできない。。。
リモートデスクトップからわざわざ立ち上げるという手間がわずらわしい。。。
誰か教えてください。
mathmatica5でコンピュータはXPです。

　

　

>>31
カーネルは超能力者です

ちょっとお尋ねしたいのですが……。

Mathematica6.0Student@Vistaで使ってます。
Export関数でFLV形式がサポートされた、ということで出力しようと思ったのですが、出力されるFLVファイルが再生できません。
だいたいファイルサイズからして7.4KBとビデオの大きさではありません。

入力した式は
Export["MovingPower.flv", Manipulate[Plot[x^a, {x, 0, 5}, PlotRange -> {{-1, 5}, {-1, 5}}], {a, -2, 2}]]
です。

再生環境はffdshow+MediaPlayerClassicです。
別のFLVは再生できるのでコーデックの問題ではないと思います。
「式がちげぇよﾄﾝﾁｷ」「私もできない」などありましたら教えてください。よろしくお願いします。


あと何でさくらﾀﾝなんですか。CCS板住人としてはここで彼女に会えるとは思っていなかったﾊｧﾊｧ。

なんかタグ書いてHTMLから動かさないとだめっぽいね

>>36
できればｋｗｓｋお願いします。

ﾎｼｭ

Import[FromCharacterCode[{104, 116, 116, 112, 58, 47, 47, 119, 119,
119, 46, 108, 117, 99, 107, 121, 45, 99, 104, 46, 99, 111, 109, 47,
105, 110, 102, 111, 47, 105, 109, 103, 47, 105, 110, 102, 111, 95,
101, 118, 48, 55, 49, 50, 48, 50, 95, 114, 101, 112, 111, 114,
116, 95, 116, 105, 116, 46, 106, 112, 103}]]

ごちそうさまでした

もうすぐ卒業だけど買うか悩んでる

>>42
きみは外資金融にでも内定が決まっているのかい？

>>42
英語版ならバイナリと鍵genが至る所に出回っているよ。

まぁ、しかし実際便利な道具だと思う。面倒な計算は大抵やってくれるし、いちいち本をひっくり返す必要がない。

今後一部でもフリーになる予定は？
関数型プログラミングできるし簡易スクリプト言語として便利そう

フリー化は難しそうな雰囲気だけどな……。

10マンなら買う

>>47
つStudent

今の商売が傾き始めたら考えるかもな

学生版って何時までも使えるの？

金がないといい何もしない奴は金があっても何もしない奴に違いない。
ところで、私に通貨をください。

matlabに対してscilabがあるように
ﾏﾁﾏﾁｶにもクローンの類はあるの？

sumathもsutudent版だった

学校で Mathematica を使えるなら、使おう。
しかし、 アメリカ合衆国の公用語に通じていないと厳しいかもしれぬ。

>>54
マニュアル？
豊富な例を自分の手でいじってると
大抵のことは分かる気もする

kingってほら、アレだから。

Mathematica 式の表示につきまして

こんにちは、

下記の式、ｙ１とｙ２があります。
ｙ１を評価させて、ｙ２の通り表示させるには、どうしたら良いでしょうか？



y1 = Expand[(a1*p1 + a2*p2 + a3*p3 + b*mc)^2]


y2 = p1^2*a1^2 + p2^2*a2^2 + p3^2*a3^2 + mc^2*b^2 + (a1*a2 + a2*a1)*p1*p2 + (a2*a3 + a3*a2)*p2*p3 + (a3*a1 + a1*a3)*p3*p1 + (a1*b + b*a1)*p1*mc + (a2*b + b*a2)*p2*mc + (a3*b + b*a3)*p3*mc



ちなみに、ｙ１をそのまま評価しますと、下記になり、ｙ２の通り表示しません。

p1^2*a1^2 + 2*p1*p2*a1*a2 + p2^2*a2^2 + 2*p1*p3*a1*a3 + 2*p2*p3*a2*a3 + p3^2*a3^2 + 2*mc*p1*a1*b + 2*mc*p2*a2*b + 2*mc*p3*a3*b + mc^2*b^2

>>57
これは、はまるべきでないパターン。
単に論文や、教科書が正しいかどうか見たいだけなら

In[18]:= y1=Expand[(a1*p1+a2*p2+a3*p3+b*mc)^2]
Out[18]= b^2 mc^2+2 a1 b mc p1+a1^2 p1^2+2 a2 b mc p2+2 a1 a2 p1 p2+a2^2 p2^2+2 a3 b mc p3+2 a1 a3 p1 p3+2 a2 a3 p2 p3+a3^2 p3^2
In[17]:= y2=p1^2*a1^2+p2^2*a2^2+p3^2*a3^2+mc^2*b^2+(a1*a2+a2*a1)*p1*p2+(a2*a3+a3*a2)*p2*p3+(a3*a1+a1*a3)*p3*p1+(a1*b+b*a1)*p1*mc+(a2*b+b*a2)*p2*mc+(a3*b+b*a3)*p3*mc


Out[17]= b^2 mc^2+2 a1 b mc p1+a1^2 p1^2+2 a2 b mc p2+2 a1 a2 p1 p2+a2^2 p2^2+2 a3 b mc p3+2 a1 a3 p1 p3+2 a2 a3 p2 p3+a3^2 p3^2
In[19]:= y1==y2
Out[19]= True

以上でy1=y2より十分。　

どうしてもこの形に意味があるならだけど、このy2の形はあくまで
「筆者の人間側の都合でこうなった」というのを忘れてはいかん。

Collectとかで係数をまとめればそれらしき形にはなってくるが　規則がわからんので
完全な一致をMathematicaの関数のみで達成するのは難しいと思うぞ。

バージョン6.02にUPした方はいますか？　
またアップデート自体はどのように行いましたか？

学生の頃マテマテカのplot3Dはすごいって思った。
こんな形してるんだって思った。
でもプレステのガンダムseed destinyのほうがすごくない？
動きは早いしまるで3次元空間にいるかのようだ。
plot3Dのf(x,y)曲面よりガンダムのほうが複雑な形してるじゃん。
ガンダムとマテマテカはどっちのほうがコードの量が多いんでしょうか？

比較するものが間違ってますよ。

>>マテマテカはどっちのほうがコードの量が多いんでしょうか？
お前がそういうコード書けばいいんじゃ

ﾏﾃﾏﾃｶは関数型プログラミングで相当短くなる。
ニューラルネットワークの単純なBPなら10行以内でコーディングできる。
対してガンダムは10進BASICだからかなり長くなる。

すいませんPlotで第2軸を利用するにはどうすればよいのでしょうか
一つのグラフのｙデータは100〜200程度で、もう一つは0.1〜1.5のオーダーですが、
見やすく一枚で表示したいです

フロントエンドとヘルプが良いねﾏﾃﾏﾃｶ。
フリーで良いのが出たら乗り替えたいけど、ここまで整備するのは大変そう。

Mathematica book　6版　出るような話はないですかね？

音をExportしてImportすると音が歪みませんか？

例えば、
s = Play[Sin[2 \[Pi] 1000 t] Exp[-9 t], {t, 0, 0.5}]
これで音を再生しても歪みは特にないのですが、
これを
Export["aaa.wav", s]
とし、
Import["aaa.wav"]
すると、えらく歪んでいる。どうなってるのでしょう？私だけ？

確かに歪むね、Exportされたファイルの時点で
もうExp^（-）じゃない変なピークになってる　

ここらで息抜きに一曲…
　"Think of Laura" by Christopher Cross

http://www.youtube.com/watch?v=NKuzZ7wIShw 02:42
http://www.youtube.com/watch?v=nYEWeBarLww 03:02 Live 1998
http://www.youtube.com/watch?v=0OTFLgW_9Qg 04:16　Live @Vancouver 2007/08/29
(1983 hit dedicated to the late teenager Laura)

>>53
もsutudent版だった

　正真正銘な　がくセイ! sutudent の　雛　だから　じゃ　姉ぇーの

（暫く　ロハに近いお値段で　　酷使　しやがル）

アメリカ合衆国の公用語に通じていないと厳しいかもしれぬ
との　こと。

ったくアメリカ合衆国の公用語に通じていないし

　　　 且つ　日本語も覚束ない

の条件達を　完璧に満たす　にも　拘らず

Mathematica を　学セイの分際で　酷使する

ヤツ_______,________,________ を　俎上に

最近Fortranしか使ってなかったけど
mathematicaっていつから暗号作成機能ついたの？

Reply:>>72 Mathematica なら暗号をすぐに作れる。

とりあえず二人とも暗号をうpしてみれば

了解デース
たうえのたきぶんしきょうがもきぬじばたけぬきしたきてまぬすきたよ

(*ワンタイムパッドバーナム暗号*)

In[55]:= key = ToCharacterCode["ルイズかわいいよルイズ"]

Out[55]= {12523, 12452, 12474, 12363, 12431, 12356, 12356, 12424, \
12523, 12452, 12474}

(*一応鍵の長さを確認*)

In[66]:= Length[key]

Out[66]= 11

(*よって2バイト文字11文字までに対して情報量的安全である、続いて平文*)
plain = ToCharacterCode["ルイズの姉も可愛いよぉ"]


Out[58]= {12523, 12452, 12474, 12398, 22985, 12418, 21487, 24859, \
12356, 12424, 12361}

(*一応平文の長さを確認*)

In[64]:= Length[plain]

Out[64]= 11

(*おｋ*)

(*おｋ*)

In[60]:= cipher = BitXor[plain, key ]

Out[60]= {0, 0, 0, 37, 26950, 198, 25515, 20883, 175, 44, 243}

In[61]:= FromCharacterCode[cipher]

Out[61]= "\.00\.00\.00%楆Æ掫冓¯,ó"

(*一応わけのわからない文字列になった*)

(*復号化*)

In[67]:= plain2 = BitXor[cipher, key ]

Out[67]= {12523, 12452, 12474, 12398, 22985, 12418, 21487, 24859, \
12356, 12424, 12361}

In[68]:= FromCharacterCode[plain2]

Out[68]= "ルイズの姉も可愛いよぉ"

(*一応復活しているのがわかる*)

実装から見ればかなりムチャクチャだけど　あくまで一例ということで

Mathematica v6は日本では４０万円以上で販売されてるようですね。
でもwolframのオンラインショップでは2500ドルくらいで買える見たいじゃないですか。
なんだか日本人は日本の代理店経由でないと買えないみたいですが、独占禁止法に違反してませんか？

｢サポートが不要だとおっしゃるのであれば、お客様のご自由にしてくださってかまいません。
なお、認証その他に不都合が発生したとしても、当社としては責任を負いかねます｣

当社ってどちら様ですか？

英語版ならbtに転がっているよｗ
鍵genはver4のがそのまま使えるｗ

今の為替レートだと日本での販売価格は４０００ドル以上。
４０００−２５００＝１５００ドルもサポート料と称してとるわけ？
ボッタクリバーみたいだね。

>>英語版ならbtに転がっているよｗ

btってブータンのことですか？

こんなこともできるんですね。

Resolve[ForAll[p, p > 0,
Exists[q, q > 0,
ForAll[x, x >= 0 && 0 < Abs[x - a] < q,
Abs[Sqrt[x] - A] < p]]], {a, A}, Reals]

Resolve[ForAll[p, p > 0,
Exists[q, q > 0,
ForAll[a, a > 0,
ForAll[x, x > 0 && Abs[x - a] < q, Abs[1/x - 1/a] < p]]]]]

統計処理能力が格段に低くて困る。
データをベクトルとして認識させるのも厄介だし。

Newton-Raphson法も簡単にプログラムが書けない。
EM algorithm にいたっては、てんで弱い。

プログラム形式でalgorithmを書くなら、Matlabの方が向いているのか？

統計処理ならまだイゴールとかのほうがいいかもね。
　それはそうと
データーをベクトルとして認識させるのが面倒って？
それにベクトルの記法にしたって｛｛｝,・・・・・｝たかがこれだけじゃない。　Importは知ってる？
Newton-Raphson法は書かかなくくてもFindRooｔで実装されている。
単にアルゴリズム書きたいならｃかFortranでいいんじゃないすか？　
前提がフリー環境だから文献も多そうだし

データから分布を推定し最尤推定量がダイガンマ連立方程式の解になったりすると、
厄介きわまりないってところかな。
テイラー展開と微分とが交互になると、FindRootでは解に行かれないんだわ。

そういうのはFindRoot使うよりも、NMinimizeで二乗値の最小化とかをした方が良いかもね。

まてぃまてぃか
なんで普及しないんや?

値段とアクチが原因か

普及してるのはなに？

電卓

プログラミング言語で組んでしまった方がはやい。

すみません、改行しようとすると IndentingNewLine と表示されてしまうのですが
これを表示させないためにはどうすればいいでしょうか？

すみません、evaluateって何するコマンドなんですか
調べても、式を評価しますと書いてあるだけで評価するとは何なのかまったくわかりません
具体的に教えていただければ幸いです

>>94
Holdを外す

初心者です。
i*jの行列リストの各行の任意の成分｛i1,i2｝をつかって
2次元グラフをプロットをしたいのですが、方法がわかりません。

ListPlotで作図したいですが教えてもらえないでしょうか。

すみません、どなたか>>93をお願いします。

>>97
使用環境くらいちゃんと書け

>96

行列の任意の要素を抽出
http://reference.wolfram.com/mathematica/ref/Part.html

平面のプロット(ListPlot3D)
http://reference.wolfram.com/mathematica/ref/ListPlot3D.html?q=ListPlot3D

せっかくなのでエロい図形でも書いてみないか?

まず言い出しっぺが女性器くらいは書いてくれ。

アレの数学的構造ってまだ明らかにされてないよね
写真とるか絵師が天才的に描くかしかない以上
数式処理ソフトでアレを再現したらノーベル生理学賞は間違いないと思う

96ですが、n×3のリストデータを使ってListCoutourPlotをしたいのですが，
;
listcoutourpot:ntri ; データが矛盾する三角形分割を生成しています．データを摂動して有効にすることができます．
;

というエラーが出てプロットで来ません．
plotrange->{full,full,automatic},
datarange->automatic
でプロットしています．

どうしたらいいでしょうか．

>103
n=1とか。

>102
http://amatias.cocolog-nifty.com/diary/2003/12/post_14.html

ｎはきちんと入力されていると思います。
１つの列の要素の末尾すべてにバッククォーテーションがついています。
これが問題なのでしょうか。
バージョンは6.0です


5.0で実行すると[2*2以上の長方形行列ではありません」とでます

DataRange->Automaticだと、三角形ができるかどうかはデータに依存する。

大学辞めてMathematicaが使えなくなったのでフリーの数式処理ソフトを使いたいのですが
Maxima以外に有名なソフトはありますか？

最強！！数学ソフトウェア(maple, mathcad etc)
http://science6.2ch.net/test/read.cgi/math/1051275230/

無料ならMapleかなー。

>>109
MAPLEは30万以上するだろ

あれ？じゃあ他の無料の数学ソフトと間違えてるわ。

違法なの覚悟で言えば、
Mathematicaのクラックバージョンが最強だと思う。

おれは Mathematica よりも maple が好きだぁ、
なんてこのスレで言ってみる。

学生版のmathematicaってパソコン買い替えるたびに買い直さなきゃいけないんだっけ？

Mapleのstackmatrixとかconcatに相当するコマンドってなんですか？
因みに、
stackmatrix は行列 A、Bを縦に積んで新たな行列
　　　A
　　　Ｂ
を、
concatは行列A,Bを横に積んで新たな行列
　　　[ＡＢ]

を作ります

>>114
<<LinearAlgebra`MatrixManipulation`
AppendColumns[A,B]
AppendRows[A,B]

>114
あるいは、下でも可。
（http://reference.wolfram.com/mathematica/ref/Join.html）

行追加：Join[A,B]
とか、

列追加：Transpose[Join[Transpose[A],Transpose[B]]]
深さレベル指定すれば、Joinだけでもできる？試してない。

>>113
マジか…
買おうと思ってたんだが止めようかな

>>108
2000年頃、MuPADが学生や研究者には無料で使えて、
簡易版Mathematicaという感じで使ってた。
が、今は有料化へ。30日間の無料試用版があるとはいえ、何だか悲しい。

唯一の希望はMaxima、ライセンス上、常に無料であることと、
LISP（Common Lisp）・Emacs、そしてTEXによる美しい表示が可能なので、
使いこなせれば楽しいのではないでしょうか。

Maximaって会社で業務用に使ってもいいんだっけ？
GPLなら出力物は著作権表示とかなしで自由に使えるんだよね？

bisonとかいう例外は除けば

age

ゆとりでスマないんですけどが、
学校の課題で以下の問題をMathmaticaでやらなければいけないんですが、
やりかったがわからないので教えて下さい

2次方程式ax^2+bx+cの近似解を、解の公式を用いて求めよ
ただし、ひとつの入力セルのみ用いるので入力セルの1〜3行目で各係数a､b､cへの代入を行い、4行目以降から解の公式を用いて結果を出力するように作成する
各係数は自由に設定してよい。NSolveやFactorを用いてはいけない。

ゆとりにもほどがあるだろ…

HELPファイルくらい読めよ…

これ宿題か？
これ30秒かかるか？・・・・瞬殺だろう

>121
MathematicaのHelpは使いやすいので、使い方を覚えるといい。
まず、何のコマンドを使うかを知るのが、一番難しい。

ver5.2？ ver6.0？ 学校なのでまだ5.2かな、5.2ならHelpを立ち上げて、例えば「方程式」とか思いつく単語を入れてみる。
すると、SolveとかNSolveとか、それっぽいものが出てくる。
例えばSolveを見ると、使い方を言葉でごちゃごちゃ説明しているが、一番わかりやすいのは一番下にある「他の例」。
「他の例」を開くと、そのいろいろな使い方の例が挙げられている。
いろいろな例を見ていくと、値の代入の仕方とか、関数の使い方もわかってくるはず。

「NSolve, Factorを用いてはいけない。」が、これもヒント。
NSolve, FactorのHelpの中には、関連する命令についても記載されているので、それらをクリックしてみればいい。

自由なディラック粒子のラグランジアンに、局所ゲージ不変性を課す事によって、電磁相互作用を導き出すことができます。これを、mathematicaを使って計算したいのですが、プログラムはないでしょうか？

物理の解説書だったら数式くらい載ってるでしょ。
それをそのまま書けばいいじゃん。

>>125
�dクス
ちなみに、ver6.0

ところで、馬鹿なオレが確認したいんだが　「解の公式で出てくる」のは近似解か？

厳密解じゃないのか・・・・・・？

式によるだろ。

自由なディラック粒子のラグランジアンに、局所ゲージ不変性を課す事によって、電磁相互作用を導き出すことができます。これを、mathematicaを使って計算したいのですが、プログラムはないでしょうか？

2006年のMathematica Conferenceに行った人いる？

6.0.3
updater 上げ

ライセンス不要のMathematica 4.0.2のユーザーだがvistaで動いてくれますか？

>>121
Mathmaticaなら板違いです。他逝ってください。

こんにちは、
風船から、亜鈴の形になり、また亜鈴から風船になる。
これを繰り返すようなアニメーションは、どうしたらできるでしょうか？

<<Graphics`ParametricPlot3D`
For[a1=0,a1£2,a1++;
r0=1;

ParametricPlot[{{r1=r0+a1*LegendreP[2,Cos[q]];r1 Cos[q],r1 Sin[q]}},{q,0,2Pi},AspectRatio®Automatic];

SphericalPlot3D[r0+Re[a1]*LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint®{-3.382,0.096,0.076},ViewVertical®{0,1,0},ImageSize®300]
];

>>136
すいません。文字化けしました。下記の通り修正します。
こんにちは、
風船から、亜鈴の形になり、また亜鈴から風船になる。
これを繰り返すようなアニメーションは、どうしたらできるでしょうか？

<<Graphics`ParametricPlot3D`
For[a1=0,a1<2,a1++;
r0=1;

ParametricPlot[{{r1=r0+a1*LegendreP[2,Cos[q]];r1 Cos[q],r1 Sin[q]}},{q,0,2Pi},AspectRatio->Automatic];

SphericalPlot3D[r0+Re[a1]*LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint->{-3.382,0.096,0.076},ViewVertical->{0,1,0},ImageSize®300]
];

>>137
すいません。やっぱりまだ文字化けしてました。下記の通り修正します。
こんにちは、
風船から、亜鈴の形になり、また亜鈴から風船になる。
これを繰り返すようなアニメーションは、どうしたらできるでしょうか？

<<Graphics`ParametricPlot3D`
For[a1=0,a1<2,a1++;
r0=1;

ParametricPlot[{{r1=r0+a1*LegendreP[2,Cos[q]];r1 Cos[q],r1 Sin[q]}},{q,0,2Pi},AspectRatio->Automatic];

SphericalPlot3D[r0+Re[a1]*LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint->{-3.382,0.096,0.076},ViewVertical->{0,1,0},ImageSize->300]
];

何度もやればいいじゃん

>138

a1++の部分を、a1+=0.1に。
ParametricPlotか、SphericalPlot3Dかのどちらかにする。

実行すると２０枚のグラフィックスができるので、
右端のひとつ内側の”]”部分をダブルクリックすると、ひとつに統合される。
統合された画像をダブルクリックすると、アニメーションが動き出す。

Mathematicaの素因数分解の関数FactorIntegerは、
100%正しく素因数分解出来ていると考えて良いので
しょうか？約20桁の整数の素因数分解を行っています。

ttp://reference.wolfram.com/mathematica/ref/FactorInteger.html

PrimeQの説明の所に、10^16以上の場合は誤動作の可能性が指摘
されているので、少々気になっています。

ttp://reference.wolfram.com/mathematica/note/SomeNotesOnInternalImplementation.html#6849

>>141
Maple使え

メイプルタウンへおいでよ

教えて君ですみませんが、
Mathematica 4.2 はVistaでも動きますか？

ははは、
動くと思うかね

やってみたまえ

これって認証あるみたいだけどPC買い替えても使える？

ライセンス上では使える
しかし、
まんどくせ

正規ユーザなのに割れ

>>146
使える。

>>141
積を作って比較するとよい。

>>146
Mapleとくらべると使えないソフトだ。

な、ことはない！

こういうソフトって一回出来たらもうそんなにする事がないのでは？
というかソフトって一回出来たらもうすることあんまりないんじゃないの？
マスマティカプログラマってバージョン１の頃より減ってるんでしょうか？
エクセルだって見た目は変わったけど実質はほとんど変わってないじゃん。
セルに２変数関数の数字入れてオートフィルして３Ｄグラフ表示させると
エクセル９７のころよりは結構改善されたとは思うけど、
この程度のバージョンアップが１０年かけてやる仕事なのかなあと思うし。
エクセルの本質は全く同じと言ってもいい。
３Ｄ曲面の表示だけはまあまあよくなったと言える程度。
どうなんでしょうか。
でもエクセルの３Ｄ曲面ってマウスでグリグリ回転できないよね。
マスマティカプログラマがエクセルプログラマに転向したら
すぐに実現するのにと思うのですが。
というか、エクセル作れるぐらいのスキルならマスマティカレベルの曲面の
プログラムぐらい作れるのではないんでしょうか？
どうなのでしょうか。

こういうソフトって一回出来たらもうそんなにする事がないのでは？
というかソフトって一回出来たらもうすることあんまりないんじゃないの？
マスマティカプログラマってバージョン１の頃より減ってるんでしょうか？
エクセルだって見た目は変わったけど実質はほとんど変わってないじゃん。
セルに２変数関数の数字入れてオートフィルして３Ｄグラフ表示させると
エクセル９７のころよりは結構改善されたとは思うけど、
この程度のバージョンアップが１０年かけてやる仕事なのかなあと思うし。
エクセルの本質は全く同じと言ってもいい。
３Ｄ曲面の表示だけはまあまあよくなったと言える程度。
どうなんでしょうか。
でもエクセルの３Ｄ曲面ってマウスでグリグリ回転できないよね。
マスマティカプログラマがエクセルプログラマに転向したら
すぐに実現するのにと思うのですが。
というか、エクセル作れるぐらいのスキルならマスマティカレベルの曲面の
プログラムぐらい作れるのではないんでしょうか？
どうなのでしょうか。

>>152-153
ヒント: π, e

Mathematica v6.0.3 ｷﾀｰ

>>85
一行でかけるけどな・・・
NewtonRaphson[f_, x0_] := FixedPoint[# - f[#]/f'[#] &, x0]

収束判定をコントロールしたいならこうとか：
NewtonRaphson3[f_, x0_] := NestWhile[# - f[#]/f'[#] &,
x0, Abs[#1 - #2] > 0.00001 &, 2]

Mathematicaで計算するくらいなら、Javaでプログラム書いちゃうようになっちゃった。

俺逆だな．昔はC++でプログラム書いてたけど，最近はほとんどMathematicaだ．

オレはFORTRANの方がいい

おれは機械語

Javaで書くような計算なら数値計算だろが。こんなもんはMatlab使え

While を使って1 から整数i までの間の素数の数を数えるプログラムを作
れ。整数i は最初の行にi=4; のように与えること。

これ教えてくださいお願いします

やなこった

Mathematicaだってプログラムは組めるんだし
数学系の関数ライブラリが有り得ないほど充実してるプログラム言語だと思えばJavaなんぞに負ける道理はないはず

勝ち負けって何ぞ

素数の判定はエラトステネスのふるいでロジックを作って、
整数iを1つ1つ判定していけばOK

MathematicaやMapleは基本的に厳密解を求めたり数式処理をするためのソフト
高速な数値解とかシミュレーションするならMatlabでいい。
C/C++を使った数値演算シミュレーションは実装(たとえばハードウェアに載せる)まで見通した場合とか
極限まで高速性を求めたいとかでないと意味がない
Javaによる数値計算なんて何の意味もない

行列ノルムの計算について質問させて下さい。
現在Mathematica5.2を使っています。
A= {{1/2, 1}, {0, 1/2}}という正方行列がある時
k = 1, 2, …　に対して
A^kの1、2、無限大ノルムを求めたいです。

1、無限大のノルムについては、Norm関数を用いて求めることができました。
Norm[MatrixPower[A, k], n] (n = 1, ∞)
(もしくはTotal[MatrixPower[A, k], 2])

しかし、l2ノルムに関しては
"Norm::symb: Normは記号的要素を持つ行列を受け取りました"
と表示され、求めることが出来ません。
また、このエラーに関してのノートは記述されていませんでした。

もしご存知の方がいらしたら、具体的な求め方を教えて頂けないでしょうか。

また別の質問になりますが、以前にMathematicaを用いた際に
たとえば変数xを用いる時に、"xは正の整数である"といった制約条件を
xに対して付与出来たような記憶があるのですが
本当にできたのか、どんな方法だったか、などを失念してしまいました。
併せて教えて頂ければ幸いです。

長文な上に未熟者の質問で申し訳ありませんが
どうかよろしくお願いいたします。

>>168
A' を Aの共役複素数求めて転置したもの として
(A'.A)の固有値の中から最大の引っ張ってくればいいんだから
↑の手順そのまま式にして使えばいいんじゃね？

線形もMathematicaもそこまで詳しくないから参考程度に。

>>168
後半は Assuming などでググれ

ttp://lottery2ch.hp.infoseek.co.jp/cgi-bin/imgf/0020-PICT0027.jpg

>>169
>>170
参考になりました。ありがとうございます。
このAssumingという関数で定義したからといっても
必ずしも解を確実に判定できるという訳ではないようですね。
やはり最後は人間の判断なのでしょうか･･･。

L2ノルムに関してですが
あれから色々調べてみたところ、MathWorldにて
Matrix Norm に関しては Frobenius Norm と Euclidean norm は
同値であるとの記述をみつけました。

Mathematica では Frobenius Norm に関してはサポートされているようですので
とりあえずはこちらを用いてみたいと思います。
･･･とはいえ、kに実数を代入して Forbenius Norm と Euclidean Norm を比較しても
何故か同値にならないので、一体どこまで信用できるのかは分かりませんが･･･
(私の用い方に問題がある可能性が高いですが)

>>172
あなたの Frobenius Norm と Euclidean Norm の定義を書いてみてくれるかな？

mathematicaって，演算子（簡単なものだとx^3+d^2/dx^2)を単独で扱えるようなことってできますか？演算子の二乗とか三乗とかを関数とかに演算したいんですけど．

ぐはははは、
できるんじゃね

>174
Map使えばできそう。
http://reference.wolfram.com/mathematica/ref/Map.html
導関数はD。
http://reference.wolfram.com/mathematica/ref/D.html

D[ ] のスコープにそれらしいのがある。

CやC++使ったほうが速いと言われても
Mathematicaが組み込み関数をどう計算してるのかわからないので
今日も全てMathematicaさん任せなのです

ウェブのドキュメントは日本語で検索できないのですか?

Mathematicaによる結果をコピーしてMapleへペーストして動かすには
コピペをどのような形式にしたらよいですか？
簡単な例も教えていただけると助かります。

そりは無理では

バージョンによると思うけど、maple側からmathematicaノートブックを
変換することは出来た筈。ただし、うまく変換できるかどうかは知らない。

なんか共通形式はないんだっけ？
math XMLとか

カルキングユーザーです。
gooで質問したのですがレスが全くなく、困ってます。
http://oshiete1.goo.ne.jp/qa4245530.html
よろしくお願いいたします。

>>183
z = 4 とか z = -4 とかってなってる時の x と y の値は何になってんの？

>>184
しらね。
そんな事を調べる機能は無い。
でも10年前のマスマティカでもこんな間違いはなかったはずだが。。。
プログラムのことよく知らないんだけど、
３Ｄ表示機能なんてクラス化されて再利用できるんじゃないの？
考えられるバグの理由としてはどう思う？
プロットした点同士を結んでいけば１〜ー１からはみ出るはずがないんだが。。。
どんなアルゴリズムで書いたらこんなグラフになるんだよ。
というかブログに画像アップしといたんだけど、ブログへのリンクが消されてた（ーー；

リンク消されてたから画像はこれね。
http://blog.goo.ne.jp/frozenmelon
ところで、３Ｄのアルゴリズムってそんなに難しいの？

そんな、どこのものとも知れない怪しいソフトを使う方が悪い

なんでMathematicaスレでMathematicaとは全然関係ないソフトの質問するの？

>>185
Mathematicaでやりゃいいじゃん。
ついでに虚数について勉強してみれば？はみ出ることもあると思うぜ。
今回のがそのせいかは知らんけど。

形式的にはカルキングの質問をしているが
実質的にはマスマティカに関係のある問題提起をしているのだよ。
それはplot３Ｄ関数だ。
この関数を作るのは難しいのかと言う問題提起がそれである。

>>190
ゴクリ･･･。とっても･･･奥が･･･深いです･･･。

Mathematica 6 for Linux 日本語版をopen suse10.3 にインストールしたら
まともに動かずpremier supportに質問したらsuse 10.3はサポートしてないと
言われた。どのdistributionをサポートしてるのかと聞いたら
http://www.wolfram.com/products/mathematica/platforms/
を見ろといわれたが、「kernel 2.4 以上の主要なもの」とあるだけ。
（実はこれはすでに見ていて、だったらsuse 10.3はＯＫだろうと思ってた）
誰か「このdistributionでちゃんと動いてる」ての知ってたらおせーて。

ウチは　大学のFedoraで結構ちゃんと動いてますよ

Fedoraも公式サポートはされてない。

相談頼みます！今大学一年数学科で学校でマスマティカの授業あるんですが夏休みの宿題でマスマティカで興味深い絵（線だけででも可）を書けってのが出てるんですがまだなんもやってません(>_<)なんかかんたんな例みたいののってる本ないですか？誰か教えてください。

数IIICの教科書に載ってる極座標で書いてある図形一般

カージオイドやアステロイド

>>195
mathematicaのヘルプみると一杯霊が出てるけど、コピペじゃだめか？
１９６の言う霊もヘルプに出てそうだけど。

>>197レスども。プログラム書いてあるおもしろい図がある本ないですかね？たぶん丸写しでもばれないと思います。誰かと完全にかぶったら終わりだけど・・・

マンデルブロとか、フラクタルでいいんじゃない。
処理に時間がかかる係数をつっこめば、遠目には絵にしか見えない。

よく考えたら、大学1年で「興味深い絵を書け」ってことは、
「数学の情報を調べろ」ということが課題じゃないのかい。

人に聞いて解決するより、自分で調べた事や感じた事をレポートした方が主旨に沿う気がする。

>>198
http://demonstrations.wolfram.com/
みるとさらに山ほど例がでてるよ。ソースもみれそうだし。

みなさんレス色々ありがとなんとか頑張って書いてみます。フラクタルにしよかな・・・

ばかたり!

ここはドラえもんか、アンパンマンの顔を書いて
高評価をねらえ

それは難しすぎて自分にはムリっすよ。書いてみたいけど

「コックさん」のような簡単な図形なら数式で何とかなりそう。
直線と円が描ければいいんだし。

助けて〜(>_<)
ちゃんと自力で教科書読みながらやってるんだがなんかうまくできない。もぅなんかを丸写ししたい

>>206
痔力が理想だが、少しくらい霊におすがりしてもバチは当たらないよ。

>>201みたんですが例ってプログラムものってます？見たんですが英語ばっかでどこ見ればいいのかわからない・・・
誰かたすけて

>>208
気になる図形をクリックして、文中の「SHOW SOURCE CODE」や「DOWNLOAD SOURCE CODE NOTEBOOK」

英語が読めないことはいいわけにならない。

>>209さんありがとうです。ただあんまりレベル高いの提出したら怪しまれるかな？

あんまりレベル高いの提出したら怪しまれるかな？あ
んまりレベル高いの提出したら怪しまれるかな？あん
まりレベル高いの提出したら怪しまれるかな？あんま
りレベル高いの提出したら怪しまれるかな？あんまり
レベル高いの提出したら怪しまれるかな？あんまりレ
ベル高いの提出したら怪しまれるかな？あんまりレベ
ル高いの提出したら怪しまれるかな？あんまりレベル
高いの提出したら怪しまれるかな？あんまりレベル高
いの提出したら怪しまれるかな？あんまりレベル高い
の提出したら怪しまれるかな？あんまりレベル高いの
提出したら怪しまれるかな？あんまりレベル高いの提
出したら怪しまれるかな？あんまりレベル高いの提出
したら怪しまれるかな？あんまりレベル高いの提出し
たら怪しまれるかな？あんまりレベル高いの提出した
ら怪しまれるかな？あんまりレベル高いの提出したら
怪しまれるかな？あんまりレベル高いの提出したら怪
しまれるかな？あんまりレベル高いの提出したら怪し
まれるかな？あんまりレベル高いの提出したら怪しま
れるかな？あんまりレベル高いの提出したら怪しまれ
るかな？あんまりレベル高いの提出したら怪しまれる
かな？あんまりレベル高いの提出したら怪しまれるか
な？あんまりレベル高いの提出したら怪しまれるかな
？あんまりレベル高いの提出したら怪しまれるかな？


どこのヘボ大だよ

お願いがあります。どなたか暇でマスマティカ使いこなせる方ドラエモモン書いてそのプログラム書いて
ここにはっていただけませんか？かなりがんばってるんですができなくて・・・
頼みますｍ（＿＿）ｍ

ませマティか、クレヨンしんちゃん
で検索しろ

213色々検索したんですがでてこなくて(>_<。)

すみません。質問なのですが、
Importでファイルを読み込む際に、01.dat,02.datというような
連番のファイルを順々に読み込ませたいのですが
いい方法はないでしょうか。

ファイル名のパスを文字列として収納すればいいと思うのですが、
その方法がよくわかりません。

>>212
data = Import["http://d.hatena.ne.jp/images/diary/s/survival/2004-11-21.jpg"];
Show[data];

>>212
FromCharacterCode[{12489, 12521, 12360, 12418, 12435}]

>>215
>ファイル名のパスを文字列として収納すればいいと思うのですが、
FileNames[ DigitCharacter .. ~~ ".dat"]
はどう？

ﾏﾃﾏﾃｶ

ちょっと買うのを待てまてぃか

>215
StringFormをToStringするんだ！

質問です。

.csv から、複数列の表データを読み込んだとして、その中のいくつかの
列を入れ替えたりしたいのですが、どうするのが一番いいんでしょうか。

awkで前処理すると吉

>222

Sort
http://reference.wolfram.com/mathematica/ref/Sort.html?q=Sort&lang=en
行列操作
http://reference.wolfram.com/mathematica/guide/PartsOfMatrices.html
http://reference.wolfram.com/mathematica/tutorial/AddingRemovingAndModifyingListElements.html

>215
ToStringのアプリケーション２
http://reference.wolfram.com/mathematica/ref/ToString.html

最近５から６に乗り換えたんだけど、
６のドキュメントセンターって記号の検索できなくね！？
例えば「/.」とか「:>」とかってみんな検索できる？
５のヘルプブラウザはできるんだが

6ってNMinimizeとかの最適化機能が5より落ちてる気がする。

さくらﾀﾝ！！ﾒｯ！！

マティマティか
はVer 4で十分だ

安く売らんかい!

>>229
俺は Ver3.0 だ。文句あるか。

Ver3は
御気の毒としか申しようがありません
ほんとうに御愁傷様です

>>231 同情くださり有難うございます。

で、君は Mathematica をどんな風に使ってるのかね。

ver3.0 では「御愁傷様」だという例を一つで良いから見せてみろ。

バグが多くて
よく落ちた

上司から
承認しましゅって言われた

>>234
おめ。
何系の会社？
どこの販売代理店から買うの？
MapleやMatlabじゃなくてMathematicaの理由は？

マテマテカ4のスツーデントなんだけど、もう学生じゃないんだけど、
２週間たってからＸＰの復元機能つかってインストール前に戻せる？

ご自由に

ライセンス違反だが、目をつむろう。

mapleはともかくmatlabって競合するもの？

mathematicaもmatlabも両方使ってる身としては
これらを比較する感覚が分からない

数式処理と数値計算だからなぁ…

まあやろうと思えばMathematicaで数値計算ができるし、
Matlabで数式処理もできる。結局どっちのウェイトを重視
するかだよ。

Matlab, Maple, Mathematica, Maxima, Mathcad
数学５強ソフト：MA5

質問があります。
ノートブックでの編集時に表示形で表示させる方法を教えて頂けませんでしょう
か。ギリシャ文字をパレットより入力してしばらくは表示名で表示されていたのですが、なにかの
拍子に完全名での表示となってしまいました。形式変換や表示形式をいじくってみたのですがうま
くいきません。

>>243
Ctrl+Shift+tのことか？

Mathematica6に関して、

ParametricPlot3Dなどで描いた曲面をイラストレーターに貼り付けたいのですが、
ただコピペをしただけでは画像として貼り付けられてしまい、
一度epsファイルで保存した後、イラストレーターで読み込まないと
線の集まり（言い方が変ですみません）として読み込めません。
Mathematica5の時は出来ました。

形式を選択してコピーの項目をいろいろ試したのですがうまくいきません。
epsファイルに保存しないでコピペをする方法を
ご存じの方がいらっしゃいましたら教えてください。

コピーの形式を選ぶんじゃないの?
試さずに書いてるが

ああ、やってみたけどダメなのね

正規に効いてみたら
メーリングリストとか

　　　　　　　　_....._{{ 〃
　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-<
　 　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）)
　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ.
　　　',　',　|Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　 ﾞ',
　　　 ヽ､', l:!Kノ}. 　 　 f:_.)iﾞi: ﾘ !　l　ル
　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉ ﾞ ,ﾘ l ﾘ'´
.　　　　 ',|!!､ 　 　r‐┐ 　 ｀ ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣
　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　 ＜　 Mathematica 7リリースおめ！
　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　 　 　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､
　　/＼＼　　　　',　　 }　　 //｀ヽ

奇数バージョンだから　期待しないな

　　☆ ﾁﾝ
　　　 　 　 ☆　ﾁﾝ　　〃　 ∧＿∧ 　 ／￣￣￣￣￣￣￣
　 　 　 　 　　ヽ　＿＿_＼（＼・∀・）＜　6.1まだー？
　　 　 　 　 　 　 ＼＿／⊂　⊂＿)＿ ＼＿＿＿＿＿＿＿
　　 　 　 　 　 ／￣￣￣￣￣￣￣／|
　　　　　　　　| ￣ ￣￣￣￣￣￣:|　:|
　　　　　　　　|　　　　　　　　　　　.|／

Mathematica7・・・アイコンがきしょくわるくなったな（笑）
俺的には4の虹色のアイコンが傑作だったと思う
今回のやつは完全にウイルスだ

ところで7ってwhat's newはあるけど(http://www.wolfram.com/products/mathematica/newin7/)
まだ製品はリリースされてなくね！？

ver3なんですが、PlotVectorFoieldでベクトルの軸を
太くしたいんですが、どうしたらいいか知りませんか？

Field
だろ

ヘルプみりゃ
太くできるんじゃね

見ても線の長さを変えるパラメータはありますが、
太くするのは見当たらないんですが。

typoすまそ

追加。

描いた図をEPSでLaTeXに張り込むと
線が細すぎてきれいに表示されないので困ってます。

うるさい。

まったり汁！

A={
{1,1,-1,-1},
{I a ,-I a, -b, b},
{e^(I ka), e^(-I ka),-e^(I kd)e^(-kb),-e^(I kd)e^(kb)},
{I a e^(I ka), -I a e^(-I ka),- b e^(I kd)e^(-kb),b e^(I kd)e^(kb)}

の行列式を

A[[1,1]] Det[ A[[{2,3,4},{2,3,4}]] ]
+A[[1,2]] Det[ A[[{2,3,4},{3,4,1}]] ]
+A[[1,3]] Det[ A[[{2,3,4},{4,1,2}]] ]
+A[[1,4]] Det[ A[[{2,3,4},{1,2,3}]] ]

で計算しました。答えは出たのですが、最終的には

Det[A] = 8I a b e^{I kd} {Cos[kd] - Cos[ka] Cosh[kb] - ((b^2-a^2)/2a b) Sin[ka] Sinh[kb] }

と整理できるはずなのですが、うまくいきません。
どうすればいいですか？

>>258
ExpToTrigとFullSimpifyでできない？

eは\[ExponentialE]を使ってる？

>>259
ありがとうございます。できました。

＞eは\[ExponentialE]を使ってる？

・・・・使ってませんでした。これが原因か。
４次のままで行列式が計算できなかったのも orz orz orz

初めまして、最近mathematicaを使い始めた初心者です
学校の課題でPrimal-Dual法を用いて最小費用流を求めるプログラムを書け
というのが出たのですが、現在"Residual Flow"の生成が上手くいかない状況です
どなたか詳しい方いましたら教えていただけないでしょうか？
なおVerは6.0です

>>261
うまく行かないといわれても、どういう実装をしようとしていて、
どの部分が分からないのかを言ってくれないと、さっぱりだよ。

>>262
グラフを使って最小費用流問題をPrimal-Dual法で解きたいのですが
残余ネットワークの生成の仕方が分からないということです

>>263
それは >>262 の質問に対する答えに全くなってないよ。

「どういう実装」と聞いているのは、例えばどのようにデータを管理しているか、
とか、そういうことを聞いている（グラフをどう持っているかとか）。

「どの部分」と聞いているのは、残余ネットワークの定義が分からないのか、
残余ネットワークを数式でどう表したらいいかわからないのか、
定義も式も分かっているけど Mathematica での書き方が分からないのか、
とかそういうのを聞いている。

そういう意味だったんですね・・申し訳ありません
一応今の所は下のような感じになっています
g0 = {{1, 2}, {2, 3}, {1, 3}, {1, 4}, {4, 3}, {4, 5}, {3, 5}};
g = FromOrderedPairs[g0];
h0 = Table[Random[Integer, {1, 100}], {7}]
h1 = SetEdgeWeights[g, h0];
h2 = SetEdgeLabels[g, GetEdgeWeights[h1]];
capacity = Table[Random[Integer, {2, 15}], {7}]
h3 = SetEdgeWeights[g, capacity];
costlist = {h0, ToOrderedPairs[g]};
TableForm[costlist];
capacitylist = {capacity, ToOrderedPairs[g]};
TableForm[capacitylist];
flowlist = {Table[flow, {Length[g0]}], ToOrderedPairs[g]};
TableForm[flowlist];

連投失礼します
残余ネットワークのどの部分が分からないかということに関しては
「残余ネットワークを数式でどう表したらいいかが分からない」です。

あともう一つ分からない所があるのですが、mathematicaのグラフの辺に
２つ以上の重みをつけることは可能なのでしょうか？
質問ばかりで申し訳ないですがよろしくお願いします

>>266
数式で書けない限りMathematicaでも書けない。
その部分は教科書なり何なりで勉強してらっしゃい。

二種類以上の重みは、「重みのペア」を一つの重みだと考えればいいよ。

>>267 お返事ありがとうございます。
もう少し自分で勉強してがんばってみます。

%ってひとつ前の出力だけど、
ひとつ前でなくてひとつ上の出力を参照する記号はありますか？
ctrl+shift+Lではなくて%みたいに短いと見やすいんですが、
%で繋いでドンドン積んでいって途中で間違えるとやり直すのが面倒だよぬ

%1など、n番前の出力結果というのがあるよ。

%1, %2などは、n番前ということではなく、n番目の出力結果なので、それでいいはず。
http://reference.wolfram.com/mathematica/ref/Out.html

でも、何度も使うようだったら、変数に記憶しておくほうがいいよ。

みなさんレスありがとう。使ってみます

7でアンドゥは立派になった？

6でフロントエンドのパーサがおかしくなった時は
新しいセルに手で打ち直すしかないのかな？

ドキュメントセンターのもっさり感は７で改善された？

みんな7に興味津々w

だって無料セミナー受けちゃったんだもん。

cpuの片方占領して重くなるのは何とかしてほしい

新しいバージョンだもん
もちろんバグがてんこ盛りだよん!

n番目の完全数を返す関数だと思うのですが4以上を入力すると上手くいきません
関数自体を自分で書き換えるしかないでしょうか？
f[n_] :=
 Nest[If[
        If[# == 0, 0,
            If[
             Mod[#2, #] ==
              0, # + #0[# - 1, #2], #0[# - 1, #2]]] &[#, #] ==
         2*#, #, #0[# + 1]] &[# + 1] &, 1, #] &@n

日本語がおかしい
まずどっからもってきた関数なのか説明しろ

>>280
動かしてないが見た感じ定義をそのまんま書いたっぽい
再帰が深すぎるんじゃない？
再帰の限界を環境変数でいぢればいい

平均m, 項数Nの二項分布をN→∞でポアソン分布になるという計算を
させようとしたのですが、

Limit[Binomial[N, r](m/N)^r (1 - m/N)^(N - r), N -> Infinity]

とすると

Gamma[1/N+1+O[N]^3]で真性特異点が見つかりましたと言われて計算できません
Binomial関数を階乗になおしても同じです。

MathematicaではGamma関数を含むと極限計算ができないんですか？
それともなにか工夫をすれば上記計算ができるようになりますか？

KETpic という、Mathematica の図を picture 環境にできる
ソフトがあるのはこの掲示板でも何回か見たので使ってみた
のですが、横軸の 1 と縦軸の 1 が同じ長さになってしまいます。
Mathematica に入力したのは以下のようなものです。

Get[ "・・・/ketpic6.m"]
Get[ "・・・/ketpicapp6.m"]
XMIN = 0;
XMAX = 5;
g1 = Plot[x^2, {x, XMIN, XMAX}];
setwindow[{XMIN, XMAX}, {0, 25}];
openfile["foo.tex"]
openpicture["0.5cm"];
drwline[g1];
closepicture[1];
closefile[]

・・・のところは本当は ketpic6.m, ketpicapp6.m
のディレクトリの絶対パスです。
TeX では input で読み込んでいるだけなのですが
図は縦横比が 25 : 5 で縦長になっています。
Mathematica では AspectRatio というコマンドで変えるので
ketpic6.m, ketpicapp6.m の AspectRatio->Automatic
という部分を消してみたのですが変わりませんでした。
使っている人がいたら教えてください。

>>283
ver6だと求まったよ

7が出るってメールが来たから何かと思ったら、ベクトルの可視化ってあの程度の事か。
わざわざやる需要あんの？

ver6がいかにも未完成だったから完成度を高めてくれるだけも意味がある。

×　だけも
○　だけでも

いや、
さらにバグ増加です

結局、7はどうなの？ｲﾝｽﾄする価値あり？

まずは買ってから悩めよw

プレミア　だから
無料でくれるらしい

>>292
俺もプレミア入ってるがそんな案内まだきてないぞ？
どこの代理店？日本電子計算？HULINKS?

おれんとこにはwolfram researchから直接日本語のメールが来た。

「ぐだぐだ言わずに買わんかい！」って内容か

いや
「バージョンアップできるのは二つ前までで終わりだ
　二つ前はマイナーを入れて数えるんだぞ」
って奴

7の日本語版でたな。
英語版から１ヶ月経たずに出すとは翻訳チームもなかなかのもの。
ただし、誤字脱字やふさわしくない翻訳も散見されるが。

んで、肝心の内容はというと、俺に言わせれば素晴らしいの一言。
きれいなグラフがサクッと作れるし、並列計算機能が標準搭載されるようになったし。
個人的にびっくりした機能が数列を何項か与えると一般項を返す関数。
数列データベースを参照してるとのことで、例えばｎ次元球の体積だと
５項ではだめだけど６項与えたら一般項をちゃんと求めてきた。すごい。

それはすごい。クイズ解くのも楽になるなｗ

テーラー展開の一般項を求めるのにも使えそうだな。

Reduce関数も強化されてる？

誰か、mathematicaの起動時のメニューバーの長さを変更する方法を教えてくれ。

しごく

>>299
それ専用の関数がある
http://reference.wolfram.com/mathematica/ref/SeriesCoefficient.html

新機能、楽しそうなものが豊富じゃないか。
年度末調整で買うか。

torrentであるだろ。日本語版はどうかしらんが。

少しでも開発者にフィードバックするためにお金を払って買いたいよね。

それは賛成する
しかしアクティベーションがうざい

再インストールしたら使えなくなった。

数十万のソフトだからなあ。torrentでおとして万一ばれたら怖い。

そういえばv7のインストール時に
シングルローカルライセンスかネットワークライセンスを
聞いてこないや。
これが原因か？

これって？

むしろ数十万だからじゃないか？

Mathematica Lightバージョンとか、計算機能を一部制限して、安いのが出ればいいね。
20000円とか

学生なら生協で買えばそれくらいの値段だろ？
社会人？

学生で買う奴いるんだろうか。大抵、学校のパソコンに入ってるだろ。

だったら母校の生協いって買えばいんじゃね！？

パッケージで売ってるんだっけ？

売ってる
要組合員証

v7はかなり良い！

v4-もっさり
v5-やや改善
v5.2-サクサク！
v6-もっさり＆メモリ消費
v7-高機能でサクサク！

これで、5.2&6.0両方使いだったが、v7だけに移行できる。

レポありがとう。やっとbtで落としたﾏﾃﾏﾃｶをインストールする決心がついたよ。
週末に一念発起してｲﾝｽﾄする。

btって日本語あるか？

btのってstudent.verじゃないのであるっけ？

こんどのは日本語版を買う価値がある
ヘルプの日本語はよくできてる

やっぱ落とすやつのが多いのかｗ
そりゃそうだよなｗ

大多数の人（大学・研究所含）は買うでしょ．
落とす奴のほうが多くなったら開発なんて続かない．

落としたやつも開発者への感謝の気持ちは忘れんなよ．

学生と教員は安く買えるからな。
俺もVer3持ってるけどアカデミック版だからこそ買えた。
97年だったけど初めてインターネットのありがたさを実感したわ。
それまでどこで売っているかさえわからなかったから。

というよりもっと簡単な方法無いの？
ドングリだっけ？　ＬＷ３Ｄのやつとか

ドンブリだろ
あんなものはすぐクラックされる

悪血が最狂
アドビをみろ

古いバージョンいれてると
v7入らないっぽい。
OS再インストールしないとだめかも。

Windows7β1にVer.7入れたらきちんと動いた
autorunのファイルでトロイを感知した
拾ってきたやつだからかなｗ

v7の学生価格っていくらでしょうか？
どなたかご存じですか？

生協で聞けば教えて
もらえるんでしょうか？

アクティベーションはクラックされないの？

不正使用に関する話は別のところでやってくれ

＞＞３３３
おれは学生だっつーの！

やはりMathematica使うのやめるわ。
Mathematica使ってるのって
思い出してみると、>>333のような
馬鹿ばっかりだもん。

>329
ver6と一緒に入れているけど、うちの環境では大丈夫。

ver6+画像処理アプリを使っており、ver7で統合されたのは良いのだが、コマンド名が結構変わっている・・・
同じコマンド名のものは２重定義されているとエラー出るし、アプリ部分は結構書き換えが必要。

それでもver6+アプリよりは快適なので、この際ver7に乗り換える！

Maxima使えばいいのに

Mathematicaの6,7 によって
わざわざMaximaなんてものを使う意味は無くなった

金出す価値がある

論文中の式の検証につかったりするので
商用のほうが若干安心できるんだよね．気持ちの問題が大きいけど．

企業がOracle使いたがるのと同じかｗ

>>336
画像処理アプリって何？ウルフラムが画像処理専用のアプリを出してるのか？
リンクはってくれ

>>339
研究者だったら所属する大学か組織が買ってくれないのか？

>>341
もちろん研究室でも買ってるよ。

　　　　　　 l　　l　　　 ｀′　　　　　　　,、　　 `ｰ'´　/　ヽ /`ヽ
　 　 　 　 /　 /　　　　　　　　　　/ﾌ　| |　　r-、　　 `ー'　`ｰ′
　　　　　 L__/　　　　　　　 r"> //　　| |　　 ヽヽ
　　　　　　　　　　　 ,.._ 　 / /　{ {　　 |｜　　　〉〉　,-、
　　　　　　　　　　　//　　{ ｛　　い 　 | |　　／ノ　/ /　 ,、
　　　　　　 　 　 　 { {　　 ヽヽ、 ヽ`'二_二' ／　 / /　 / ﾉ
　　　　　　　　 /7　ヽヽ､ 　 ｀'-ﾆ,/´_=;;;;=_ ヽ｀ニノ　 / /
　　 　 　 　 　 { {　　 ｀ｰ-｀ニﾆ /ｒ'ﾞ_ ヽ /　_`',｀ｰ-‐'´ノ 　 　 r'７
　　 　 　 　 　 ヽヽ、　　　　_,.ノl=・=-.∴.-=・=!'-ﾆ二＿＿／／
　　　　　　｢｀!　 ｀'-｀ニﾆ二-ｧ |　`ｰ'_ !｣ _ー' |ヾー―---‐'´　 ,-､
　　　　　　l　!　　　　　 ,. -'´ri |ｌ　 ‘~`ｰ'~′ ｌ| r‐｀ﾆﾆ'ｰ- ､_／／
　　　　　　ヽ､ヽ--―'ﾆ ‐'´　{ヾ ! ,r;'ﾆｰﾆヽ、!l/}ヽ、　 ｀`ー‐'´_　　
　　　　　　__　｀`'' ''´　 　 ／,ﾄﾐ! iｒ,!-┴-!､ヾ, !'/＼ヽ､_ 　　／ノ　　
　　　　　　い、　　　　 ／／/ハ||′-―-､｀|||' l ヽ ｀ヽ､二.／　　>>61こらまんじゃ夫馬いっちょ
　　　　　　 ＼｀ｰ―'ﾆ´‐',／ !　||!/　　　　ヽl||　lヽ ＼　　　,-､　　おまんじゃえーじゃよなｗ
　　　　　　　　,､￣　,. -'_／ !　 ||!_＿＿＿__l||　 l　＼｀ｰ‐'´ ノ
　　　 　　── ＝ ニ 　 /=､｡｡｡｡｡｡｡｡｡｡｡｡｡｡｡｡ｒ=､ヽ
　　　 　　　─ ＝ニ三　（◎ ヽ-─────（◎ 　)
　　　 　　　　　　ノ◎､　 |＼　 ＼　　 　　　　/　 ／ |　 /◎､　
　　　 　　　　　(_,ｒへ　｀ｿ　 /> ◎) 　　　　 (◎く| 　レ'　,ｒへ )
　　　　 　─ ＝ ニ　 ＼◎'/　／　　　　　　 ＼ ヽ､◎／
　　　　 　　　　　　　　　 ノ ／　　　　　　　　　　＼ ヽ　　　
　　　　　 　─ ＝ニ三 ( ◎(　　　　　　　　　　　　 ) ◎)
　　　　　 　　　　─ ＝　 ｰ､_ら　　　　　　　　　 ⊂､_,ｒ´

大学生協はv7の代理店じゃなくなるんだね。
どこが代理店になるんだろ？
価格上がっていくのかな？

v7からはダウンロード販売になるんだね！

>>344-345
ソースは何？
ほとんどの大学はMathematica導入してるから
生協が代理店じゃなくなったら、他の代理店はウハウハだな
日本電子計算とヒューリンクスぐらいしかないだろ

>>346
ネットで検索してみてください。
どこも代理店契約は今年の３月で終わります。

ただ、代理店じゃなくなっても販売は
継続すると生協のＨＰには書いてありますけど。

やはり今後はＤＬ販売がメインになるみたいです。
StudentバージョンもDLで購入できるようですね。

>>347
学生かどうかの判断はどこでするんだろう？

学生証うｐ？

なんとなくパッケージであのスパイキーを拝まないと買った気になれないよな
生協の代理店契約が切れるからって、ＤＬ販売がメインになるなんてことはないだろ
そもそもＤＬ販売なんてものに代理店なんて必要ないし
ウルフラムとの直接のやりとりなんだから
ということはやっぱり国内の大手代理店と契約することになるんじゃないか？

Ｍａｔｈｅｍａｔｉｃａを全部日本語で購入できるようになればもっと売り上げが上がるだろうに

購入できないんだっけ？代理店通しても
つうか、Mathematica使える人は、購入するぐらいの英語は使えるだろ。

もしそうならわざわざ日本語のヘルプ作ったりしない

あ、購入するぐらいの、か
それでも敷居が低いに越したことは無いだろう
「どうしても欲しい」人は多少のことは何とかするだろうが
「なんとなく欲しいかな」ぐらいの人にとっては障壁は低ければ低いほどいい

生協が去れば
日本電子計算とヒューリンクスだけか
線香はだめだし

試しに拾ってきたものをインストールしてみたのですが、ライセンス番号・MathID を入力する画面になりません。
パスワードは keygen.exe があったので、これで作ると思うんですけど。
私は WAREZ 的なものはあまり詳しくないので分かりません。
これでは安心して購入できません。
教えてください。

>>352
ヨドバシカメラで売ってるよ。
（うちの近くのヨドバシには置いてあった）

あ、悪い。
うちの近くはビックカメラだった。

>>356

落ちているのはニセ物だよ
今回のは買わないとダメなのさ

一般用を一番安く買おうとしたらこの辺？
http://item.rakuten.co.jp/plusyu/4988727009970/
HULINKSは定価しか乗せてないけど、見積もり取ったら幾らぐらいに
なるんだろ。

３元連立２階非線形偏微分方程式ってマスマティカで解けますか？
解けるソフト教えてください。

そんなもの一般には解析的には解けないでしょ。
Mathematicaはある程度の範囲のものは解いてくれるけど、
それ以上は高望みしすぎだよ。

解くって解析的に？数値解だけでＯｋ？それで全然違うぞ

解析的に解いてこそ数学
やっぱそういうのだとマスマティカ最強？

解析厨がわいてきたのか

>363
数値的にです。

昔からmapleの方が得意じゃないかな

数値解ならMathematicaで十分。
ただし計算速度とか精度にこだわるなら別のものを使うべき。

数値解だけでいいなら、
フリーのScilabやMaximaでも十分だよな。

>>369
はい。ScilabやMaximaでも十分な用途の場合は十分です。

俺みたいな素人にはMathematicaの方が使いやすい

Mapleは今市使いづらい。

MapleはグラフィックスがMathematicaより美しいけど、
ヘルプの使いづらさがね。
それに最近の傾向として、
開発元は純粋な数学関連の開発をおろそかにして、
金を儲けやすい工学への応用に走ってる。
トヨタも関係してるのかもしれないが、
これはもう終わりだね。
Mapleさようなら。

数式処理でマスマチカに勝てるソフトなんて世界のどこにもないだろ
数値計算は他にもっとすごいのいくらでもあるけど

http://www.eb.waseda.ac.jp/murata/junichi.mimura/knowledgh.html
こういうふうに行列の微分したいんだけどできる？

>>374
でもごくまれにMathematicaで解けなかった常微分方程式が
HPの数式処理電卓で解けたりするんだよなｗ

>>345
できる。例としてまず行列を用意する。

mat = Table[\[ExponentialE]^(j k x), {j, 1, 4}, {k, 1, 4}]
% // MatrixForm

次にこの行列を格納したリストmatに微分演算をする
ans = \!\(\*SubscriptBox[\(\[PartialD]\), \(x\)]\ mat\)
ans // MatrixForm

微分演算ができている

mathmaticaで計算してでてきた式をjavaやtexに貼り付けたいのですがなにか方法はありますか？
出力をコピペしてもとても読めません。
あと計算結果から、同じ項をなにか文字で置き換えさせたりとかもできるなら知りたいです。
例
out> {Σz^3+y^2　+ Σx+y^4　- Σx^3*y^1　}/{ Σx-6+y^4　+ Σx+y^4　Σ・・}
みたいなのを
{Σz^3+y^2　+ A　- Σx^3*y^1　}/{ Σx-6+y^4　+ A　Σ・・}
といった感じです。
いずれも厳しそうな感じですが、よろしくおねがいします。

まず、MathMLで出力する。
http://www.mathmlconference.org/2002/presentations/smirnova/
これを使う。

なんていうのでどうでしょう?

>>378
TeXへの貼り付けは TeXForm[...] で吐いて微調整すれば何とかなる。
Javaに貼り付けるというのは意味がよくわからんけど、
式をソースコードの一部に使うなら CForm で吐いて調整、
コメントに貼り付けるなら InputForm かね。まあこのあたりで調べると色々あるよ。

初心者ですみません
複数の値を覚えさせる時、関数を使うのとリストを使うのとどちらがいいのでしょうか？

リストです

ありがとうございます

順番に値を覚えさせて行く時でもa[5]=48、a[6]=115、……、とするよりAppendToを使った方がよいですか？

>>379>>380
�dくす！研究がいっきにすすみそうです。どうしても手計算じゃすまないところを
mathematicaに頼ってみたものの項が多すぎて詰んでいました。

f,h,i,j,k,lはそれぞれx,y,zの関数で非線形であり、x,y,zがそれぞれu,vの関数。
∂f/∂u+∂g/∂v＝０
∂i/∂u+∂j/∂v＝０
∂k/∂u+∂l/∂v＝０
を境界、初期条件を入れてx,y,zについて解きたい（数値解析）したいんですが、
差分してニュートンラプソン方が良いのでしょうか？
マスマティカで可能でしょうか？
フォートランとかがいいんでしょうか？

ふぉーと　らんらんらん
ふぉとらんらんらん
ふぉーと　らんらんら〜ん

リストとして
{a,b,c,d,e,f}
があってそれを
{ {a,b}, {b,c}, {c,d}, {d,e}, {e,f} }
のようにして返す関数を作りたいのですがどうすればいいでしょうか？
組み込み関数にあるかと思ったのですが見つからなくて

Table[{#[[i]],#[[i+1]]},{i,Length@#-1}]&[{a,b,c,d,e,f}]

>>389
ありがとうございます！
できました

質問です。
data=Import[original];
data2=Drop[data,*];
data3=オリジナル関数[data2];
ていうのを手順を何回もくりかえしたいのですが、
この一連の処理をさらに関数化したいのですが、うまくできません（ていうかできない？）
どなたか教えてください。data3をrerturnしたいです。
NewData=新しい関数[orginal]のようにしたいです。

新しい関数[orginal_]:=(
data=Import[original];
data2=Drop[data,*];
data3=オリジナル関数[data2]
);
みたいな感じ？

392>>
391です。有難うございます。
新しい関数[orginal_]:=( data3=オリジナル関数[Drop[Import[original] ,*]; ]);
みたいに入れ子にすればできるみたいです。関数宣言のなかに変数(dataとかdata2とか)使うとむりみたいですね。
仕様？

>>388
Partition[{a,b,c,d,e,f},2,1]
でできるだろ？

>>393
Block,Module,With
あたりを使え

386です。
DSolveでとけますかね？
解けるやつ教えてください。

>>396
なんでマルチしちゃったの？

組合せ Combination nCr: Binomial[n,r]
順列 Permutation nPr: FactorialPower[n,r]

この事実に俺は日本で一番最初に気づいたので、ここにメモしておく。

Mathematica使い始めた頃を思い出すが、
ヘルプでCombinationとして探すも出てこず、
色々調べた結果、２項係数として関数になっていることがわかった。
Permutationは探せど探せど出てこず、
結局自分で定義して使ってた。
後に分かったことは、CとかPとか使ってるのは日本だけということ。

大学行くと習うけどね。2項係数。[2,n]。

nCr、nPrがそれぞれ特殊なケースでしかない。

Solveで、解の範囲を制限して答えを導くことはできますか？
例えば
Solve[2 x^2 +x-2==0, x]
とすると、
{{x -> -1.28078}, {x -> 0.780776}}

とでますが、あらかじめ、x は正の範囲に限った答えを出したい（つまり後者の答えだけ出したい）んです。

マニュアル、ヘルプを見ているのですが、よくわかりません。
初歩的な質問で申し訳ありませんがよろしくおねがいします。

＞CとかPとか使ってるのは日本だけということ

マジで！　どういう経緯でできたんだ

>>401
日本だけ，ということはないよ。
nCkも，一応海外でも通用する記号。

ただし，国際的に主流の記号ではないので，書物で見る頻度は，列ベクトル記法の方が圧倒的に多い。
かの Donald.E.Knuth の Concrete Mathematics（訳書：コンピュータの数学） には，
「熱心な査読者を除いて，呪わしい記法 nCk は，誰も使いたくない代物である」
と書いてある。

nCk という記号が海外で嫌われる理由は，Cというシンボルだけが目立っていて，
それより大切な n と k が”虐待されて"小さく書かれているから。

ちなみに，C_n^k という記法（Cの右下にn，右上にk）もある。
中国では，この記法が主流だそうだ。

nCrが列ベクトル記法で
(n)
(r)
と表されることはわかったんだけど、
じゃあnPrとかnHrとかは海外ではどうやって表現されてるわけ？
nCrは「ｎ個のものからｒ個とった組合せ」を表現しているわけだし、
nPrは「ｎ個のものからｒ個とった順列」を、
nHrは「ｎ個のものからｒ個とった重複組合せ」を表現してる。
もしたとえばnPrが
(n)
(r)r!
とかで表現されてるんだとしたら、これは本来の意図が伝わりにくくて、nPrの方が優れた記法だと思うんだが。
nHrもしかりで、
(n+r-1)
(r)
としか表現できないんだとしたら（単一の表現がなくて）、nHrの存在ってありがたくね！？

ところで
「ｎ個のものからｒ個とった重複順列」
ってシンボルとして存在してたっけ？
なんかあったような気がするんだがど忘れしてしもた

nHrって知らないぞ。

自分で使う気が無いから覚えてないけど、二項係数の (n r) の ( ) のかわりに
< > や [ ] や { } を使うと重複組合せなり何なりを表すルールだったはず

>>405
Homogeneous Product だから nHr という記法が使われる。
高校の教科書には掲載されてないけど，「大学への数学」なんかでよく使われる。

>>403
Knuth の本だと，
(n)_r で nPk
(( n r )) で nHk
[ n r ] で第一種スターリング数（n個のものを r 個の円形テーブルに並べる方法の数）
{ n r } で第二種スターリング数（n個のものを r 個の組に分ける方法の数）
を表している（本当は行でなく列記法で）。

>>404
ただの n^r では？

ひょっとして大学とか未だにC記号使ってるのか？
とっくに絶滅したと思ってた…

質問です。
マスマティカで多重総和を計算させるとき, 外側の総和の添字の記号って内側の総和の添字、
中の関数にも影響しますかね？
うまく文字で説明できてないと思うんで式を書いておくと、例えば
Σ_{i=0}^t Σ_{k=0}^i F(k)
っていう式があったときに、外側の総和のiが2だったとすると、
内側の総和のiも2になって順次計算されていくか？ということで。

C言語でいうfor文みたいになっているのかなーと。

そう

>>407
あれっ、高校の教科書に載ってたと思うけどなぁ、重複組合せ。ちなみに俺２５。

＞[ n r ] で第一種スターリング数（n個のものを r 個の円形テーブルに並べる方法の数）
多項式の係数としての意味は知ってたんだが、こんな解釈ができたのか。知らなかった。

＞ただの n^r では？
いや、そうなんだけど、シンボルがたしかあったなぁと思って。
調べたらΠだった。nΠr

>>411
25歳だと，複素数平面があった旧課程ですね。
当時の学習指導要領
http://www.nicer.go.jp/guideline/old/h01h/chap2-4.htm
によれば，
>(イ)　組合せ
>［用語・記号］nPr，nCr，階乗，n！
なので，少なくとも高校範囲内ではなかったはず。
ただし，出版社によっては教科書の「発展事項」コーナーなどに掲載されていた可能性も否定できない。

>調べたらΠだった。nΠr
そういえば，そんな記号があったね。思い出しました。

>>408
1年生の時に高校との橋渡しとして習うけど、
すぐに2項係数のベクトル表記を教えられ、「なるべくこっちを使ってください」と言われた。

授業では2項係数をほとんど使わなかった。

質問です。

z[x_, y_] := x^(1/2) y^(1/2);

Plot3D[z[x, y], {x, 0, 10}, {y, 0, 10}, PlotPoints -> 50,
ViewPoint -> {0.70, -3.0, 0.75}]

z=3 で xy平面に平行な平面を描きたいのですが
どうやって描けばよいですか？
どなたかお願いします。

{z[x, y], 3}
とするのです。

>>415
完璧です。 綺麗にお絵描きできました。
ありがとう。

もうひとつお願いします。
z=3 で xy平面に平行な平面のみをプロットする方法も
お願いします．

自己解決しました。

{z[0, y], 3} とするのですね．

名前間違えました．
>>416
でした．

Plot3D は 2 変数の関数のプロットです。
z[x, y] = 3 というのは x, y の関数なので
z = 3 は
Plot3D[3, {x, -1, 1}, {y, -1, 1}]
でいいのです。

>>419
完璧です。
ありがとう．

5.2のPlot3Dは関数の複数指定ができなくていちいちグラフィックスオブジェクト
として保持させてからShowする必要があった。
最近のMathematicaはかなり改善されてるが、まだまだ改良の余地はある。
がんばれウルフラム！

質問です。

>>414 or >>416 に関連して
ｙ=5 で xz平面に平行な平面のプロット方法
をどなたかお願いします．

そういうのは ContourPlot3D を使ったほうがいい。
Mathematica 持ってないけど。

>>423
感激です。
素早いレスありがとう。
調べてみたいと思います。

曲面ぐりぐりするためのおまじない忘れた。
なんだっけ？

<<RealTime3D`

>>425
さっさとバージョンアップしろ
７はすごいぞ
3D回転もデフォルトだ

4*2=8

>>426
ありがと。
>>427
すごいって、細部まできれいってこと？

すごいなんてもんじゃない
細部まで綺麗だし

そのうえに
あんなことも
こんなこともできるんだ!

なぜか数学のソフトの説明に見えない不思議ｗ

株価のチェックとかできるし、世界各地の気温や果ては人ゲノムの情報まで
全部わかる。確かに何のソフトかわからんｗ

v7を入れたら、v6はとっとと捨てたらいいが、v5.2は残しておいた方がいいかも。
v7はv5.2よりはメモリを消費するので、v7で計算できないことがv5.2ではできることもある。

グラフィックの質が格段に良くなっているので、計算点数減らしてでもv7を使いたくなるが。

>>433
マジか
下位ヴァージョンのカーネルも選べるようにしてくれたらいいのに

JavaのIDEみたいな機能だな。
確かに欲しいけど。

メモリがネックになって計算できないって
おまえらどんだけ少ないメモリ積んでんだよ

メモリを使うようなコードを書くと計算終わらん

方程式をSolveしながら
Manipulateしても
ちゃんと着いてくる

なかなかやるな　7

最近、4年ぶりにMathematicaを使い始めたが、ver 7は
３Dのグラフィックスが使いやすくなった。
マウスで3Dの視点をお好みの位置に決めた後、
その位置がどういう Viewpoint なのかを教えてくれる
機能があれば良いのだが．．

viewpointって(1,1,1)でも(5,5,5)でも同じじゃない？

ちがうよ。

core i7のPCに
windows 7のOSで
Mathematica7これ最強

ListStreamPlotって version7 の売りだと思うんだが、ものすごくメモリ食うのな。
↓はWeatherDataで取得した風速ベクトルのデータだけど、たかが9点を
元にした流線を描かせただけでメモリを1.8GBも消費して（俺の環境では）
落ちる。なんとかならんもんかね。

ListStreamPlot[{{{139.767, 35.683}, {0.740407, -4.74255}},
{{139.833, 35.633}, {12.5036, -7.3423}},
{{139.78, 35.552}, {8.58425, 7.34851}},
{{139.662, 35.837}, {0.349356, 1.35571}},
{{139.528, 35.672}, {4.63184, -1.2594}},
{{139.401, 35.737}, {0.896234, -1.07553}},
{{139.348, 35.748}, {4.86188, 4.16199}},
{{139.384, 35.428}, {-0.593812, -1.26783}},
{{140.386, 35.765}, {-9.48151, -6.14744}}}]

>>442
>ListStreamPlot[{{{139.767, 35.683}, {0.740407, -4.74255}},
>{{139.833, 35.633}, {12.5036, -7.3423}},
>{{139.78, 35.552}, {8.58425, 7.34851}},
>{{139.662, 35.837}, {0.349356, 1.35571}},
>{{139.528, 35.672}, {4.63184, -1.2594}},
>{{139.401, 35.737}, {0.896234, -1.07553}},
>{{139.348, 35.748}, {4.86188, 4.16199}},
>{{139.384, 35.428}, {-0.593812, -1.26783}},
>{{140.386, 35.765}, {-9.48151, -6.14744}}}]

確かに落ちる

２群に別れてるからかな

たとえばf(f0,t)=sin(2*Pi*f0*t)というような定数f0を含む変数tの関数をグラフに描くときに
f0に具体的な数字を入れずに描く方法はありませんか？
f0はそのままグラフ上に反映されて欲しいです．

なにをいってるのか意味不明

0< f_0 < 1 ∈ R
みたいな定数を表示させたいんじゃない？

Mathematica Home Edition
http://www.wolfram.com/products/mathematicahomeedition/
・価格は295ドル
・機能はプロフェッショナル版と全く同じ
・商用利用・授業での利用等は不可（宿題に使うのは可）
・ダウンロード販売のみ
・現時点では英語版のみ。アメリカとカナダからのみダウンロード可

個人でライセンス持ってる俺の立場は…

>>447
>（宿題に使うのは可）

宿題に使ったら
バカになるじゃないか

>>449
アメでは大学院のクラスでもhomeworkが出るよ。日本で言うレポートのこと。

高価なMatlabでも安価な学生版があるわけだし、Mathematicaにあってもおかしくない。

>>448は、スチューデント版とは別な。スチューデント版はもっと安い。
学生でもなく、企業で使うのでもなく、家で趣味として数学を楽しむ人
がターゲット。ってどれぐらい需要があるんだろ。

大学で使ってて、卒業後に家でも使いたいって層？

正規ライセンスの Professional が職場にあり、自宅でHome
を使用し、その結果を職場のProfessionalで走らせる。
これライセンス違反になるの？

>>448
これって日本語版もあるの？

>>454
まだ今月の5日に発表されたばかりで、日本語版は無い。
そのうち出るかもしれんが。

日本語版は5万円と予想してみる

>>453
ライセンスがProfessionalとHome分あるなら何も問題無いでしょ。

>>456
円高メリットが、、、

逆にサンプルの数を増やしたら上手く行ったわ。やっぱり場所が
偏在してると駄目なんだろうか。

data = Cases[{{WeatherData[#, "Longitude"],
WeatherData[#, "Latitude"]}, WeatherData[#, "WindSpeed"]*
({Cos[#], Sin[#]} &[WeatherData[#, "WindDirection"]])} & /@
WeatherData[{"Tokyo", 50}], {{_Real, _Real}, {_Real, _Real}}];
g1 = ListStreamPlot[data];
g2 = Graphics[{EdgeForm[Black], LightGray,
CountryData["Japan", "FullPolygon"]}];
Show[g2, g1, ImageSize -> 400, AbsoluteOptions[g1, PlotRange],
AspectRatio -> Automatic]

書き方が悪かったみたいですね．
たとえばy=axのグラフを書こうとするとaに何か値を入れなければなりませんよね．
このとき，aが未定定数のまま，プロットができないかという質問です．

mathematicaにa,2a,3aと目盛りを振ってもらいたいわけですね？

>>461
違うと思います。
>>460
a が分からないと右上がりかどうかも分からないので、グラフの形は分かりません。
a に具体的な値をいろいろ当てはめて同じグラフにプロットすることはできます。

>>460
Manipulate[Plot[Sin[a x], {x, 0, 2 \[Pi]}], {a, 1, 10}]

ttp://commons.wikimedia.org/wiki/File:KleinBottle-01.png

このページのコードを実行してみたんだけど、
切れ目が入ってしまって、画像のような出力は得られない。
どうやらu=πのところっぽいんだけど、なんでだろ？

うちの大学がMathematicaからMaximaに乗り換えたんだけど氏ねよ

>>464
UnitStepをIf文に直してuの範囲を0〜4Piにしたらそれらしいものに
なったけど、まだ何か違うな

あ、if文に直すというのは、uの判定を元の式どおりにするってことです。

>>465
定員割れしてフトコロ寂しいんだろ。察してやれ

>>461
そうです．
>>462
言われてみれば，そりゃそうですね．
>>463
こりゃいいですね．ありがとう．

>>465
さようならｗ
生協とか>>448の日本語版を個人で買ってこのスレに戻って来ることを期待しとるぞ

>>464
最初問題なく出力されたけど、
起動しなおしてもう1回やってみたら切れ目が・・・
なんじゃこりゃ

Home Editionいいな
でもさ、通常より格段に安いとは言ってもそこそこの値段がするものなんだからパッケージが欲しいなあ
カスタマーサービスに連絡しろって書いてあるけど実際もらえるのかな

すいません、どうしても判らないんで質問させてください
行列式とか、複数の数式をまとめてカウントする方法って無いでしょうか？
\!\(\*
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{"8", "5", "0", "6", "0"},
{"3", "0", "0", "2", "0"},
{"0", "0", "4", "0", "0"},
{"2", "0", "2", "2", "2"},
{"7", "0", "0", "0", "0"}
}], "\[NoBreak]", ")"}],
Function[ BoxForm`e$,
MatrixForm[ BoxForm`e$]]]\)
こんな感じの行列式の１つの数をカウントする方法をご教授願いたいです
CountだとList要素で、しかも１つのlistしか出来ないみたいなので

意味がよくわからないけど、例えば0より大きい要素の数を数えるとかだったら、


Count[\!\(\*
TagBox[
RowBox[{"(", "", GridBox[{
{"8", "5", "0", "6", "0"},
{"3", "0", "0", "2", "0"},
{"0", "0", "4", "0", "0"},
{"2", "0", "2", "2", "2"},
{"7", "0", "0", "0", "0"}
}], "", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]\), a_ /; a > 0, 2]

すいません、わかりづらかったですよね。。
例えば
{5,0,0,0,0}という数字の羅列があったとして
０の数字の個数が知りたい場合は
Count[{5,0,0,0,0},0]にすればいいですよね？
それでは、{5,0,0,0,0}、｛6,2,0,0,3｝、{5,2,0,0,0}などの乱数リストが
１００個あった場合、０の数の合計を知りたい時どうすればいいのか判らないのです

>>475
リストのリストなら、例示したように、Countの第３引数のレベルを２にすればいいんじゃないの？
Count[リストのリスト,　a_/;a==0,　2]
で０の数がでるでしょ？

>>475
Countの説明をよく見た方がいい。第３引数のlevelspecとか。

>>477
解決しました！ありがとうございました！

みんなどれ落としてるの？いくつかあるけど容量が微妙に違うような？

買ってるよ

うむ、この機能で30万円なら安い

ListPlot で自動で凡例はつけられないのですか？

> うむ、この機能で30万円なら安い
orz

質問ーー。
Mathematica 4 をVistaでも使ってる人いる？

age

一体いつになったら日本でhome editionが買えるんだよ？

英語で問い合わせる作業に戻るんだ

コンセントの形で国を分類するとか、ver7の例題が楽しすぎる。Mathematicaはどこを目指してるんだ。

MS-Officeじゃない？

Mathematicaはスゴイよ。最適化の問題を解こうと思ってMatlabでmpsフォーマットのファイルを
読み込むライブラリか何かネットにないかな？と思って探したけどなかった。しかし、
Mathematicaはデフォルトでサポートｗ

ListPlot で簡単に凡例できませんか？
ver6 です
Mathematica のサイトにある

Needs["PlotLegends`"]
ListPlot[Table[{x, f[x]}, {f, {Sin, Cos, Log}}, {x, 0, 10, 0.5}],
PlotLegend -> {"Sine", "Cosine", "Log"},
LegendPosition -> {1.1, -0.4}, Joined -> {True, True, False},
PlotMarkers -> Automatic]

↑すらできません
PlotLegend は Plot, ListPlot で使えるらしいです

http://reference.wolfram.com/mathematica/PlotLegends/ref/PlotLegend.html

Mathematicaの凡例って無駄にでかい余白と言い無駄にでかい影と言いあんま
使う気しないんだけど、なんであんなデザインなんだろう。

>>491
では１つのグラフに複数のデータをプロットする場合は
どのように区別させているのですか？

>>491
確かに。
俺は Ｔｅｘｔ とか Ｌｉｎｅ とか使って手動で凡例作ってる。

理論物理学者で数学ソフトを扱う米Wolfram ResearchのCEO、スティー
ブン・ウルフラム氏が、新たな検索技術「Wolfram|Alpha」を限定公開した。

Wolfram|Alphaは5月に公開予定で、現在は少数のユーザーにのみ公開している。
このプロジェクトは、アルゴリズムや自然言語解析などを活用して、「円周率の300
けた目は？」といった事実に関する質問に回答する「計算型ナレッジエンジン」を目
指している。ウルフラム氏は具体的な仕組みを明らかにしていない。

Wolfram|Alphaのデモに参加した検索エンジンの専門家ノバ・スピバック氏は、この
検索エンジンは「Googleのように答えを含む文書を返すのでもなく、Wikipediaのよう
な巨大なデータベースでもなく、Powersetのように単に自然言語を解析して文書を
検索するのでもない」としている。

Wolfram|Alphaは、「ティンブクトゥはどこの国にあるか」「水素原子の陽子は幾つか」
「今月のシアトルの平均降水量は？」「ISSはどこか」「GOOG（Google）の株価が300ド
ル以上だったのはいつか」といった事実に関する質問に回答してくれるという。

スピバック氏は、Wolfram|Alphaは幅広い質問への答えを提供し、データベースから
答えを探すのではなく回答を計算するという点でGoogle よりもずっとスマートだとし
ている。「WebにとってGoogleと同じくらい重要かもしれない」「Googleの人がこれを見
たら驚いて、欲しくなるだろう」と賞賛している。

http://headlines.yahoo.co.jp/hl?a=20090309-00000039-zdn_n-sci

やりたいことは分かったが、サービスとして定着するかどうかは微妙だな〜

本当にそんなにすごいものであれば普通にGoogleに買収されそうで

スティーブン・ウルフラムって天才科学者だぞ

どこが天才?
物理に寄与した?

天才がいれば儲かれば誰も苦労しない。

スティーブン・ウルフラムが宗教を創ったとしたら
俺は入教してしまうだろう

数学なんて一種の宗教じゃんw

人類への念の無許可見による介入を阻めば、哲学もよりよくなるだろう。

折れのいたところの数学教員はカルトの教祖みたいなｷﾁｶﾞｲばかりだった。

ウルフラムの業績を列挙しようぜ

Mathematica を開発した

以上

他には

Mathematica で儲けてウハウハ
弟も役員でウハウハ
一家安泰

やせました

http://ja.wikipedia.org/wiki/%E3%82%B9%E3%83%86%E3%82%A3%E3%83%BC%E3%83%96%E3%83%B3%E3%83%BB%E3%82%A6%E3%83%AB%E3%83%95%E3%83%A9%E3%83%A0

ver6からデバッガがついたのはいいんだけど、使い方がいまいちわからない。
どっかにチュートリアルとか使い方の解説はないかな。
ヘルプ見ても個別のコマンドの短い説明が載ってるだけで、ノートブックでの統合的な使い方が
わからない。Weｂで探したけど、不思議なほど見つからないんだよね。

は?

>>510
Mathematica Navigatorという本に入門程度は書いてあった。
日本のMathematica本はみんな古くて、ver6以上扱ってるのないな。

ヘルプが充実したから
もう本は出ないか

ヘルプは一通り知った後辞書代わりにするにはいいんだけど、
一から勉強したいとか、面白い機能がないかパラパラ探すにはまだ
本のほうがいいな。その点」Wolframの本はよかったんだが。

ウルフラムの本も
ソフトに入ってるんじゃなかったっけ

virtual bookというのはついてるけど、ウルフラムの本がそのまま入ったという感じはしない。、
トピックの寄せ集めって感じ。実質的な中身は同じなんだろうけど、本に比べると読みづらいんだな。

ところで、home editionは日本ではでないの？

Mathematicaって使えたら就職には有利？

あんまり。Mathematica使えることを有利と評価する会社なら、
論文一本でも多く出しておくことのほうが有利に評価される。

>>518
全然有利にならない。RDBMSのSQL使える方がよっぽどウリになる。

なるほど…

うちは業務でMathematicaを使うことが相当な頻度で
あるから少しは有利になるかもしれんけど、操作法なんて
すぐ覚えられるからなあ。むしろ使いこなすだけの数学的
知識が必要。

今年度からver7になったんだが、windows版はウチの研究室のどのマシンでもがまったく動作しない
調べて見ても同様のケースについては全然報告がないので当惑している。
今度の学会に、今から他言語で一から組み直すにはもう時間がない・・・・
かといって本学のVer6.0の契約は今月いっぱいまでだろう、・・・・・正直困った。

具体的には

・フロントエンド、カーネルが起動するが1+1すら演算中のまま答えを返さない、ずっとそのまま。
・ドキュメントセンターを開こうとする、あるいはネットワークのテストを行おうとするとうまくいかない。
・関数自体の動作は起動のたびにjavaw.exeというプロセスを止めると回復する、
　これでもドキュメントセンター、ネットワーク系ははおかしいまま。
・Administratorでインストール、起動した場合はすべてうまく動作するようである。
・カーネルだけの動作ならどの権限でもうまく動作する。

何か指針がありましたら、ご教示ください・・・・。よろすく・・・・。

人への念の無許可見による介在を阻止せよ。

>>523
インストールして一回目の起動ではそういう現象も有った。
>1+1でずっと演算中…

pc再起動したらなくなったよ
ちなみにWindowsXP 64bit

確かver6は平行して使えると思ったけど??

ありがとうございます。

>pc再起動したらなくなったよ
インストール以来再起動することは何度もあるんですが
やはりうまくいかないままです。

>>526
向こうも理由が分からないならVer6.0をもう少し使わせてもらえないか掛け合ってみるとか？

ごめんなさい、解決の報告をしてませんでした。
システム課の方と戦って、なんとか動くようになりました
プロキシの設定をうまく参照できていないようですねぇ・・・。

ネットワークライセンスか
起動のたびにいちいち確認するのな
シングルアロンならそんな問題はない

スタンドアロンね。

(* 二次元のOrdering[ ] *)
Ordering2D[mat_] :=
Transpose@{Quotient@##, 1 + Mod@##} &[
Ordering@Flatten@mat + # - 1, #] &[Dimensions[mat][[2]]]

(* テスト *)
MatrixForm /@ (test =
Flatten[Table[RandomInteger[10, {row, col}], {row, 10}, {col, 10}],
1])
Table[
{# == Sort@#} &[test[[n, Sequence @@ #]] & /@ Ordering2D[test[[n]]]],
{n, Length@test}]


変則的なやつ{{1,2},{3}}とか、3次元以上{{1,2},{{3,{4}}},5}
にも使えるOrdering関数を寛大でエロい人下さい!

これだべ
Flatten[Position[list, #] & /@ Sort@Flatten@list, 1]

>>532
すいません返事が遅くなりました。
こんなコンパクトな方法があるとは！ありがとうございました。

7.0.1

自分でパッケージ作ってそれをみんなに使ってもらいたい時は
自分のホームページで宣伝するしかないのかな？

MATHEMATICA 7.0.1
* 中核となる画像処理関数の性能アップ
* 基本数学，数学授業，文章作成のアシスタントパレットの日本語化
* 素早い画像処理を可能にする右クリックメニュー
* 新しいチュートリアル，「How to」ガイド，スクリーンキャスト
* ドキュメントに含まれる数千に及ぶ新規例題
* 向上したドキュメント検索
* Windows 7の手書き数式認識機能との統合
* 新しいgridMathematica Serverとの統合

こんにちは、

下記ＨＰに記載の計算をwindows PCで評価するのですが、
メモリーが足りなくて途中で終了してしまいます。環境は
次の通りですが、最後まで計算させる方法を教えてください。
また計算結果（グラフ）がどのような結果なのか？教えて下さい。
なお、計算時間は５時間以上必要です。
計算プログラム
http://www.geocities.jp/etctransformation/mathematica1-7.html


ＰＣ環境
Dell
Inspiron 530
Windows Vista Home Basic
メモリ　100GB
プロセッサ　Intel Core 2 Duo CPU 2.2GHz

>>537
>>メモリ　100GB
　　 。　。
　 / ／ ﾎﾟｰﾝ！
( Д )

教えてやる

Windows Vista などというクソは
100GBのメモリを乗せても
3GBとしか使えないのだ

出直してこい!

確かに「その通り」で、それは何もビスタに限らず、マイクロソフトのＯＳは
クソばかりだが、リナックスに入れ替えようとするとデバイスドライバが無い
から困りますよね。この「糞みたいな状況」は何とかなりませんかね。

猫の希望はソラリスですが、もしダメならデビアンでもエエのですが。
誰かどうしたらエエか教えて下さい。
（恐らくクノーピクスで解決なのではと思っていますが。）

おまえはマックがお似合い

マックも確か中はＢＳＤですよね。だからＸウィンドウが使えるから好きですねｗ

今、Linuxでデバドラがないって、どんな機器を使ってんだろ・・・

Win以外はエロゲができんからダメだ

いや、ノートとかだったらＮＩＣのチップさえ判らないのが時々ありますよ。
今だったらプローブとかもあるんだろうけど。
とにかくノートって安く小さくする為に変な部品を使ってますよ。

Slackware全盛のころは大変だったが、最近は無線LANでも
あんまり苦労せんがのお・・・

ああ、そうなんですか。ほんなら明日にでもやってみますわ。
そやけど本当はソラリスが一番好きなんですが。

ネットで落とせばいいけど、最近の↓

http://www.amazon.co.jp/dp/486190465X
Linux Perfect Vol.3

http://www.amazon.co.jp/dp/4048679244/
週刊アスキー別冊 さくさくUbuntu! 新バージョン9.04 CD-ROM付き

いや、金が勿体無いから何時もネットから落としていますｗ
そんでソラリスだけは今でもかなり「相性」があるみたいで、
ちょっとな〜　って感じなんですが。

でもドライバが揃い過ぎて複雑怪奇になったヤツよりかは
マシなんですかねｗ　（一番酷いのがマイクロソフトでしょ。）

>>537
このレスを見てから、2日間ぐらいMacProをぶん回して
計算しているんだけどいまだに計算おわんないよ?
まぁ計算は進んでいるみたいだから
2週間ぐらいまでなら最後まで計算させてみるけれど。

ちなみにMacProのスペックは
Core2Duo Xeon 2.66MHz x2
メモリ6GB

>>549
>Core2Duo Xeon 2.66MHz x2

ばーか
そんな石ねーよ

すまん。
Core2アーキテクチャのデュアルコアXeonって意味ね。
具体的にはWoodcrest。

良く判らないけれど（コメント文が少ないので）、散乱振幅とかの
数値計算ですか？

mathematicaでコードを高速にするコツはなんでしょうか？
Cっぽく書くと速くなる位しか分かりません

どなたか学生版のプレミアサービス価格ってご存じだったら教えてください。
HPはサービス内容の宣伝はいいだけど、肝心の価格が見当たりませんでした。

>>554
カビのはえたような時代遅れのＣっぽく書くいたって遅くなるだけだろｗｗｗ
４ＧＬライクに記述できる部分は４ＧＬライクに記述汁。

>>554
コード晒せば

>>550

長時間の計算、有難うございます。申し訳ございません。一部コードが間違っていました。修正しましたので再度計算頂きましたら幸いです。
http://www.geocities.jp/etctransformation/mathematica1-7.html
但し、５日間くらい必要だと思われます。表示した図を教えてください。

３D曲面をマウスでグリグリできるHPを作りたいと思ってます。
どうしたらよいでしょうか？
フラッシュというのでできますでしょうか？
それともjavascriptでできますでしょうか？

できるできる
なんでもできるぞ
かんがれ

ウィンドウズでも「ダイレクト何とか」みたいなヤツがあるでしょう
そんでジャバやったら関数でもライブラリでもいっぱいあるでしょうが

マウスグリグリの動的HPは
サーバーサイドでしょうか？
クライアントサイドでしょうか？
マスマティカをhtmlに埋め込むか、サーバーに置くか
という選択肢もありでしょうか？

トラフィックが重くなると困るので、ローカルのスクリプトで処理
した方がベターじゃないんですか
だから使用するポートは８０番その他だけにして、
後はセッションを張った時だけに
ローカルのＣＰＵパワーを喰う様にした方が
良いのではないでしょうか

すみません質問です。
ある陰関数をContourplotでPlotして得られる曲線があるのですが
曲線上のある点（x,y）を境に2つに分けて、片方を点線、片方を実線で
プロットしてやりたいのですが、どうすればいいでしょう。

Range[10,Infinity,10][[3]]で30を返すmyRange関数はどう書けば良いでしょうか？

NumericQ@Infinity が False だから、そこに持ってこれないんじゃないの

どなたか教えてください。
微分方程式の解き方ですが
yl''[x]==-k x (0<=x<=L/2)
yr''[x]==-k(L-x) (L/2<=x<=L)
B.C. y[0]=0,y[L]=0,yl[L/2]=yr[L/2],yl'[L/2]=yr[L/2]
これをDSolveを使って解く方法を教えて下さい。
あるいは、DSolveを使わない方法でもかまいません。
どうか、お願いします。

ユーザ会は解散と言うことですが
どうなるでつか?

なんでまてまてぃかっていつまでも日本語購入サイトが無いの

v3で作ったプログラムを、v7の体験版（Ｅｎｇ）で実行したらエラーが出まくり。
どっかに日本語ヘルプが読めるv7体験版ありませんか？
ベクトル場グラフィックをがどの程度なのかきになるんです。

web のリファレンスを使えば良いんじゃないですか

>>571
v6と違ってちゃんと日本語だったんた。ありがとう。
使用期間もあとちょっとだし、てきとーにいじって遊んでみる。

いや〜　流石だね！
ウォルフラムのヲッサンはお商売が巧いなぁ

なんかあった？

alpha?

超初歩的な質問で申し訳ないです・・・．
実験データを片対数グラフにプロットし，
ListLogLinearPlot[・・・・・,Joined->True] でデータ同士を線で結んだんだけど，データとデータの間が直線で結ばれてしまう('・ω・`)．
Excelの散布図の様に，なめらかな曲線でデータ同士を結ぶにはどうしたらイイ？

二次元の極座標成分のベクトル｛Aｒ（ｒ、φ）、Aφ（ｒ、φ）｝を
｛ｒ、0、10｝、｛φ、-π、π｝の範囲でグラフィックにしたいんですけど、

VectorPlot[]

だと、デカルト座標（ｘ、ｙ）として読みこまれてしまって、
ほしい図がでてこないんですが、どうしたらいいですかね。

PolarPlot という関数がありますが違うんですかね

普通のベクトルの図示ができんでごわす

>>577

これでどうでしょう。

rotate[x_, y_, t_] := {x Cos[t] - y Sin[t], x Sin[t] + y Cos[t]}

Ar[r_, t_] := 0
At[r_, t_] := 1;

With[{r = Sqrt[x^2 + y^2], t = ArcTan[x, y]},
VectorPlot[rotate[Ar[r, t], At[r, t], t], {x, -10, 10}, {y, -10, 10}]]

>>578
PlorPlotだとベクトルが描けないような気がします。

>>580
できました！ありがとう！
当初、極座標成分ベクトルをデカルト座標成分ベクトルに直してプロットするところまでは
思いついたんですが、逆三角関数の処理がうまくできなくて困ってました。

自分の書いたプログラムの方も
ArcTan[y/x]→ArcTan[x,y]
に直したら正しい図が出てきましたｗ

>>567

もう解決してしまったかもしれないが、一応。

DSolve[{yl''[x] == -k x, yr''[x] == -k (L - x), yl[0] == 0,
yr[L] == 0, yl[L/2] == yr[L/2], yl'[L/2] == yr'[L/2]}, {yl[x],
yr[x]}, x]
Block[{L = 1, k = 1}, Plot[If[x <= L/2, yl[x], yr[x]] /. %, {x, 0, L}]]

>>576

やりたいことはこういうことでしょうか?

data = {{1, 1.0}, {10, 2.2}, {100, 3}, {1000, 4.4}};
Fit[data, {1, Log[x]}, x]
g1 = LogLinearPlot[%, {x, MapThread[Min, data][[1]]/10,
MapThread[Max, data][[1]]*10}];
g2 = ListLogLinearPlot[data];
Show[g1, g2]

米株やってるけど、これほど安いソフトは無いと思う

高度な場所っぽいのにごめんなさい

y^2/b^2+z^2/c^2-x/a=0
この関数をxyz平面にプロットしたいんだけどどうすればいいんでしょう

sage忘れた・・・
ごめんなさい

a,b,cはどうするつもりだ
xyz平面とは何か

>>587
まあまあ、いいじゃないですか。

>>585
バージョン6以降が使えるのだったら、

Manipulate[
ContourPlot3D[x^2/a^2 + y^2/b^2 - z/c == 0, {x, -2, 2}, {y, -2, 2}, {z, -10, 10}]
, {a, 1, 4}, {b, 1, 4}, {c, -1, 1}]

などとすればよいです。

フーリエ変換について


こんにちは、
y=t^4 Exp[-t^2];をフーリエ変換すると、

y2=FourierTransform[t^4 Exp[-t^2],t,w];
Print[y2];
になります。
この計算を、まともに下記で確認したいのですが、答えが一致しません。
下記のどこを修正すれば良いのでしょうか？


y=.;
y=t^4 Exp[-t^2];
y1=Integrate[y*Exp[-I*w*t],{t,-Infinity,Infinity}];
Print[y1];


下記HPのNo2の計算を、MATHEMATICAで確認したいのですが、プログラムを教えてください。

http://www.geocities.jp/dirac_equation/img010.pdf

>>589
うーん、マジレスしないといけないんだろうか?
y1 = FourierTransform[t^4 Exp[-t^2], t, w];
y2 = 1/Sqrt[2 Pi] Integrate[t^4 Exp[-t^2] Exp[-I w t], {t, -Infinity, Infinity}];
y1 - y2

とりあえずここで止めとく。反応があれば考えるかも。

Mathematica7のWindows版で永久実行中の不具合が出た人いる？

>>590
お返事有難う御座います。

y2 = 1/Sqrt[2 Pi] Integrate[t^4 Exp[-t^2] Exp[-I w t], {t, -Infinity, Infinity}];
の積分計算がうまくできません。if==何とかが出ます。設定が悪いのでしょうか？

>>592

if==何とかじゃわからん。エラーメッセージを全部書きなさい。Mathematica のバージョンも書くこと。
もし初心者だったら、掛け算にはスペースを使わず、* を使ったほうがいい。

>>593
下記の通りです。Ver4.2です。

y2=1/Sqrt[2 Pi] Integrate[t^4 Exp[-t^2] Exp[-I w t],{t,-Infinity,Infinity}]

General::stop: 計算中，Series::esssのこれ以上の出力は表示されません．
SeriesData::argbu: SeriesDataは1つの引数で呼ばれました．3個と6個の間の引数が想定されています．

If[Im[w] == 0, ((1/16)*Sqrt[Pi]*(12 - 12*w^2 + w^4))/E^(w^2/4), Integrate[E^(-t^2 - I*t*w)*t^4, {t, -Infinity, Infinity}]]/Sqrt[2*Pi]

>>594
問題なく答えが出てるじゃないか。w の虚部が零か非零かで場合わけされてるんでしょ。

Ver5.0で

Plot[x, {x, -1, 1}]

といった具合にPlotでグラフ描画すると全部シャープ記号(#)で描画されて
めちゃめちゃ荒いアスキーアートになってしまいます。
普通に細い線でグラフを描くにはどうすればいいのでしょうか

>>596
実行環境を詳しく

>>597
kernel上でやってたのが原因でした・・
お騒がせしてすみませんでした；；

実行環境はWindowsXP，Mathematica5.0 kernelでした

市販のソフトを使わずにMathematica5.0で書いたグラフをExcelのセルにコピーする方法を
教えて下さい｡

普通にコピーしてテキストファイルに貼り付けて、いらない部分をデリートした後、
Excelでスペース区切りで開くと確かにデータは読み込めて、グラフにもすることが出来るのですが、
縦軸と横軸に何か縮尺がかかっており、また縦軸にオフセットが乗っています。
どうすれば正常に読み込ませることが出来るのでしょうか？

そんな「ソフトの使い方」なんか覚えても
数学は出来しませんねん！
それよりも自分の頭を使って何か考えたらどうですねん！

>>600
てめー、関係ない所まで首突っ込んで
何も悪くない人間まで罵倒するようになったか。
Mathematicaが仕事で必要な人間もいるんだよ。
自分と他人、他人Ａと他人Ｂが、それぞれ違う立場の
人間であることを理解できないなら人間やめろ。

あ、もうやめてるか。

オマエ何言ってんねん！
意見したらアカンのか！

ワシ、ちょっと考えたんやけどなァ、
自分と他人のそれぞれ違う立場が理解出来ひんんは
どっちやねん！
そういう事を何時もしとんはオマエとちゃうんか！
しかも此処は誰でもカキコ出来る自由掲示板やないけ！
オマエ、口に気を付けへんとエラい事になるゾ
判ったか！

とりあえず、600さんの名前が気に入らない。

>そんな「ソフトの使い方」なんか覚えても数学は出来しませんねん！
それは使う人に依ると思う。あなたが必要ないなら覚えなければ良い話。

>>599

こんなんじゃだめですか?

Table[{x, Sin[x]}, {x, 0, 2 Pi, 0.1}]
Export["graph.xls", %]

Ver 5.0 が xls 形式をサポートしていないようだったら、

Table[{x, Sin[x]}, {x, 0, 2 Pi, 0.1}]
Export["graph.csv", %]

とすればいいのかも。(Ver 5.0 が手元にないのでわからない。)

ｘ二乗検定のプログラムをmathematicaで作成せよ。

mathematicaでのカイ二乗統計量

ｘ二乗検定のプログラムをmathematicaで作成せよ。
(1) データ２項目のカテゴリカルデータ
1: 3 2
2: 2 1
3: 1 1
(2) クロス集計表を作成せよ。
(3) 独立性検定に用いるx二乗統計量を計算せよ。
(4) p値を求めよ。(x二乗分布の分布関数から)

(3)は、d1 = RandomChoice[{0.8, 0.2} -> {1, 2}, 100];(*Ａ大 1:男、2:女*)
d2 = RandomChoice[{0.4, 0.6} -> {1, 2}, 100];(*1:大学院進学 2:その他*)
d0 = {d1, d2};
d = Transpose[d0]と
入れた後に、
MeanTest[d, 1, FullReport -> True]
としました。
この結果、
{FullReport ->
TableForm[{{{1.15, 1.67}, {4.1797832761154154`, 14.177446878757825`}, StudentTDistribution[99]}}, TableHeadings -> {None, {"Mean", "TestStat", "Distribution"}}], OneSidedPValue -> CDF[StudentTDistribution[99], {-4.17978, -14.1774}]}
と出ました。

検定統計量は、mean1.15に対し、teststatが4.17978、mean1.67に対し、teststatが14.1774となりました。
検定統計量は２つでることは間違いではないでしょうか。

早急にお返事ください。

また、(4)はNeeds["HypothesisTesting`"]
ChiSquarePValue[0.3, 1]でいいんでしょうか？

こんにちは、
mathematica4.1,windows2000を使用しております。
BaseForm[200,2]を評価すると、そのままBaseForm[200,2]と表示されて答えが出ません。どうしたら計算するでしょうか？

>>607

1+1 などの他の計算は実行されるでしょうか?

>>608
>1+1 などの他の計算は実行されるでしょうか?

大丈夫です。

>>609
ってことはカーネルは動いているってことなのかな?
情報が少なすぎてよくわからん。スペルミスとかではないですよね?

すいません質問です。
Clear[dt]
Hold[dt[7]] /. dt[hikisuu_] :> hikisuu

がHold[7]になりますが、なぜ

Hold[Dt[7]] /. Dt[hikisuu_] :> hikisuu

でHold[7]とならず、Hold[Dt[7]]となってしまうのでしょうか？

Hold[Dt[7]] /. _?(Dt == # &)[hikisuu_] :> hikisuu
で無事Hold[7]になりましたが、上の原因が知りたいです。
(Windows VistaでMathematica 7.0.0です)

なんでDt[7]を7に置き換えたいんだ？
全微分ってのは、中で使われてる変数の名前に応じてその値が異なるから、
Dt[x]:>dx みたいにして使うことはあるけど、
例：http://reference.wolfram.com/mathematica/tutorial/TotalDerivatives.ja.html の最後
Dt[x_]で定義するのがそもそもおかしいな。

>>612
レスありがとうございます。例が不適切＆説明不足でした。
Dtは定義済みの関数の一つという意味で使いました。
パターンがマッチしてるのになんで置き換えてくれないんだ〜と

>>613
定義済み関数の中でも、Dtは特殊なんだ。変数の名前に依存するから。
Sin とかでやってみれば違いがわかる。

>>614
すいません確認が遅れてしまいました。レスありがとうございます。
Sinで試したらおっしゃる通りうまくいきました。定義されてるかどうかで決まるわけではなかったんですね。
Hold[Dt[7]] /. Hold[Dt[hikisuu_]] :> hikisuu
にするとどの関数でも変換するようなので、とりあえずHold付けておきます。
重ねて、レス下さった方々有り難うございました。

ForAll[x, 1/x == 1/x]はTrueとなりますが
xが0のときはどういう解釈なんでしょう？

>>616
1/0はComplexInfinityと一応定義されてるからじゃ？とおもったけど
単純に1/x == 1/xがTrueだからみたいね
ForAll[x, x == 0, 1/x == 1/x] // Trace

ComplexInfinityってDirectedInfinity[]の引数無しなんだな。
そんで方向が分からない時はイコールは未定義と。
全然知らんかった。

学生版を買ったとして、
卒業後もそれを使えるの？

読み込んだ画像の一部(○×○[pixel]）を切り出したいのですが
どう書けばいいのでしょうか。

imagecropでできるはず

一様な電場中に導体球がおかれたときの電位・電場の様子を可視化せよ
ってわかる人いない？

流体力学の分野かな。

>>622　二次元

f[z_] := 1/Log[c] Log[Abs[(z - c)/(c z - 1)]] /. {c -> (3 - Sqrt[5])/2}

ContourPlot[Abs[f[x + y I]], {x, -1.2, 1.2}, {y, -1.2, 1.2},
PlotRange -> {0, 1}, PlotPoints -> 50, Mesh -> False,
Epilog -> {{PointSize[0.05], Hue[0.5],
Circle[{0, 0}, 1]}, {PointSize[0.05], Hue[0],
Circle[{1/3, 0}, (3 - Sqrt[5])/2]}
}]

数字を文字列に変換するにはどうしたらいいでしょうか？

わかりました。

z=f(x,y)
z=g(x,y)
の２曲面の交わる線がどんな形してるかを調べたいのですが
やり方が思いつきません。


z=f(x,y)
y=2x+3
のように曲面と直線の交わる線ならやり方はわかります。

x=t
y=2t+3
z=f(t,2t+3)
とやってParametricPlot3Dすればいいと思うのですが、
曲面と曲面の交わる線の表示がわかりません。
よろしくお願いします。

>>627
ParametricPlot3D[{x,y,z}/.Solve[{z==f(x,y),z==g(x,y)},{x,y}],{z,-1,1}]
とか

Mathematica で　Ｔｅｓｌａとかのgpgpuは使えるようになった？

ぐぐると去年の１１月ごろの記事で、「Mathematicaは来年初頭に対応する」というのが
ぞろっと出るけど、それ以降のがでてこない。wolframのサイトでも見つからないし。
うまくいかなかったって感じなのかな。

ParametricPlot3D[{x,y,z}/.Solve[{z == Sin(x+y),z==x+y},{x,y}],{z,-1,1}]
とやってみましたができませんでした。。

>>630
？？？

すみません、マスマティカいじったのがひさしぶりなので
お願いします
ver4です。

Pe[z_] := Erfc[-z/Sqrt[2]]/2
LogPlot[Pe[-Sqrt[R]], {R, 5, 15}]

がほとんど直線になります
QPSKという変調方式の誤り率なんですが
教科書では曲がっておりました
なぜでしょう？？？

Arrogance goes before destruction.

>>630
何をしたいのかさっぱりわからん
ソルウ゛にz要らんやん

例えばですね、z=sin(x+y)とz=xyの２曲面の交わる線が
どういう形をしているのかを知りたいというわけです。
よろしくお願いいたします

>>636
ちょっとズルですが、以下のようなやり方はどうでしょう?

f[x_, y_] := Sin[x + y]
g[x_, y_] := x y
h[{x_, y_}] := {x, y, g[x, y]};
points = ContourPlot[f[x, y] == g[x, y], {x, -10, 10}, {y, -10, 10}][[1, 1]];
ListPointPlot3D[Map[h, points]]

わざわざありがとうございます。
しかしできませんでした。。

>>638
ContourPlot[Sin[x + y] == x y, {x, -3, 3}, {y, -3, 3}]
では？

初歩的な質問かもしれませんが、
mathematicaで出力をベクトルで表示することって可能ですか？
a={ax,ay,az}
b={bx,by,bz}
c={bx,by,bz}
d={bx,by,bz}
って入力したとして、
Dot[a,Cross[b,c]]+Dot[a,Cross[c,d]]
とかの出力は
az(-by cx +bx cy)+・・・
って感じでベクトルの成分で出力されてしまうんですが
a・(b×c＋c×d)
みたいにベクトルのままで出力されてほしいんです。
これってできますか？

>640
最後に、//TraditionalForm

よく思うんですが、
ベクトルとマトリックスをリストによって表現するのはwolframのミスではなかったか。
ベクトルの集合を作って、この集合に何か作用させようとするときに、リストを使うと簡単に行えない。

解析解を持たない微分方程式と実験データとを(最小二乗法等で)フィッティングしたいのですが
どのようにしたらフィッティングできるでしょうか？
FindFitとNDSolveを組み合わせるのでしょうか？
ご教授よろしくお願いいたします。
ちなみにver7を使っています。

>>642
明示的に列行列や行行列にするのはだめ？
{{1}, {0}} とか{{1,0}}とか。

>>644

ドキュメントセンターのFindFitの項目に例題がいっぱい載っている。
例えば

data = {{1.0, 12.}, {1.9, 10.}, {2.6, 8.2}, {3.4, 6.9}, {5.0, 5.9}};
model[a_?NumberQ, k_?NumberQ] :=
model[a, k] =
First[x /. NDSolve[{x'[t] + k x[t] == 0, x[0] == a}, x, {t, 0, 6}]]
fit = FindFit[data, model[a, k][t], {a, k}, t]
Show[Plot[model[a, k][t] /. fit, {t, 0, 6}], ListPlot[data]]

のようにすればよい。

>>645

間違えた。643への返信でした。

>>645
あまりの優しさに感動しました。
本当にありがとうございました。
またなにかございましたらよろしくお願いいたします。

あのーマスマティカって２変数関数の微分ってできますか？たとえばxa＾yをxで微分するとか

単なる偏微分だろ

マスマティカは偏微分できますか？xa＾yだったらどう入力すればいいんですか？普通の微分と同じですか？

普通の微分と同じ
D[x a^y,x]

>>651ありがとうござます。あと（y＾2）（a＾x）＋（x）（a＾（y＾2）） （a＞0）の点（0,0）における３次のテイラー多項式を出したいんですがこれってどう入力すればいいんですか？昨日やったんですがエラーばっかでできなかったんです・・・教えていただけませんか？

少しは調べろよw
Series[(y^2)(a^x)+(x)(a^(y^2)),{x,0,3},{y,0,3}]

どこかの質問スレと違って大違いだ

>>653ありがとう。ずっと｛（x,y）,0,3｝でやってました笑

マスマティカって関数の極値をいきなり出すことってできますか？

FindMaximumとかあるみたいだが、、、。
自分でもっと探してみな。

グラフにラジアン表記で目盛りをつけるにはいちいちTricksで指定しないと無理なんでしょうか？

いくつもある数を指定した数にくくりたいです。
たとえば、
(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o)
とならんでいるものを、
(a,b,c)(d,e,f)(g,h,i)(j,k,l)(m,n,o)
みたいに3つずつのくくりにしたいです。
お願いします。

Partition[{a,b,c,,,},3]

Mathematica Home Editionが日本で買えないのは、やっぱり日本の代理店とか商社がじゃましてるから？
せっかくの円高なのに、、、。

MicrosoftMathかなりいいね
カルキングはクソ。買って損した

ホームエディションて普通のと何が違うの？

機能はプロフェッショナルと同じらしい。
しかし値段は無限と有限の違いだよ。
アメリカでは２月から売られてるようだが、日本では気配の気の字も感じられないな。

http://www.wolfram.com/products/mathematicahomeedition/qa.html

Q: How is Mathematica Home Edition different from the professional version of Mathematica?
Mathematica Home Edition is a fully functional version of Mathematica Professional with the same features.

Ver．6です。
図を見ている場所を知る方法ってありますか？
ViewPointで設定できることはわかるのですが、逆に見ている状態がどこにいるのかが知りたいです。
比較検討のために必要なので、わかる人がいましたらお願いします。

>>664
マジで！？
これって知人に海外旅行のついでに買ってきてもらうとかいうわけにはいかないの？

>これって知人に海外旅行のついでに買ってきてもらうとかいうわけにはいかないの？
どうなんだろうね。
しかしダウンロード販売してるのに、わざわざそんなことしなきゃいけないのかね。
日本の代理店とかに聞くと、ウルフラム社が日本とかに売らないようにしてる、とか言うが、
日本の代理店の方が高い方しか売らないでくれとお願いしてるんじゃないかと疑っている、、。
これは公正取引法違反じゃないのか？

Mathematica も実質無料で使える方法が開発されてきたからなあ。
Home Edition すらいらないじゃんってなるかも。

http://www.hirax.net/keywords/log/IronRuby/latest
これと、そこから貼られているリンク先の記事参照。

これは本にでもまとめてほしいところだな。でもwolfram社からすれば想定外の使い方だろうから、
そのうち手が打たれそうな気もする。

Mathematicaがすぐにフリーズするようになったのですがどうしたらいいですか？
何かを計算させている最中ではなく数式を打ち込んでいる最中に止まります

HDDに不良セクタが出来てるんじゃないの？

ディスククリーンアップとデフラグで治りますか？

クリーンアップとデフラグは意味ない。

確実に不良セクタの確認をするには専用のツールを使うしかない。
スタンダードなのはHD Tune。無料版と有料版があるが無料版でOK。
http://www.forest.impress.co.jp/article/2007/02/22/hdtune.html

もしこれで不良セクタが見つかった場合はHDDをフォーマットしてから
(クイックフォーマットはダメ)再インストールするのが確実。

一応OS標準の機能としてHDDのエラーチェックがありこれで修復もできるが
これは不良セクタで読めなかったファイルを勝手にダミーで置き換える上に
どのファイルを置き換えたかレポートしないのでほとんどの場合解決にならない。

microsoftMathって曲面を回転させるときになんで縦と横の比率が変わるんでしょうか？

こんにちは、

下記は強制振動のHPです。In[4]の変数を変えて、振幅をどんどん増加させるには、変数をどのようにすれば良いでしょうか？
http://www.eng.toyo.ac.jp/~nwada/AppliedPhysicsMathVer6/AP_ForcedOscillation_No7/5ForcedOscillationVer6.html

Mathematica Home Edition の紹介ページ
http://www.wolfram.com/products/mathematicahomeedition/
の一番下に，
Mathematica Home Edition is not licensed … …
… in the U.S., Canada, U.K., and Europe.
For other territories, please contact Wolfram Research.

とあるけど，「日本のユーザだけどHome Edition を買いたい！」
と尋ねてみた人いる？

>>676
半年くらい前だけど、日本でいつ売るのか聞いたことがある。
いつからという質問には答えてくれず、...not available in Japan と
そっけなかった、、。

みんなでもっと騒げば売り出すかも、、。

そうなの。 売るつもりないのかな？

俺は、先生に買ってもらったMathematica がノートPCに
入っているから Home Edition は買わないけど、日本から
買った人いるのかどうかには関心がある。
もし、俺が本当に買うつもりなら、ダメもとで他の連中と同じ
値段で、俺にも売るようにプッシュしてはみるけど。
相手も商売だから、こちらからプッシュすれば、売るような
気がするけど。 そのために、please contact Wolfram Research
っていう一文があるんじゃないのかなぁ．

>>675
ただ拡大するだけなら線形だから入力信号のAを大きくするのとちゃうの？

共鳴に着目するなら Ω だけどね

>>670

もしかしたら、変なものをつけない？
マウスとキーボードだけにしても同じなの？
それからマウスもとりかえてみて

複雑な連立方程式をReduceで解いてRoot[…]であらわされている場合、
厳密解ではなく、近似解である可能性はあるのでしょうか？
RootをさらにToRadicalで計算した結果、無理数などで解が表されたら、それは厳密解ですよね？？

>>681　ん〜難しい問題だな。ってか質問の意図がわからん。
Reduceは厳密解を求めるものであって近似解を求めるものじゃないから解が表示されたら、それは厳密な普通の解じゃない？

具体的な問題

複素数の極限なんか素通りしていることがある？

グラフの領域の大きさを指定する方法はありますか？
縦軸の桁が増えたりすると大きさが変わるじゃないですか

Plot[x, {x, 0, 1}]
Plot[100000 x, {x, 0, 1}]

Windows版のver6を使用しています
↓にはプロットから座標を抽出するためにはCtrl+マウス左と書いてあるのですが実行しても何も起こりません
ttp://library.wolfram.com/howtos/coords/
ver5.2ではうまくいくのですが、何か操作法が変わったのでしょうか?
そうであれば正しい操作を教えて頂きたいです
よろしくお願いします

>>685

Plot[x, {x, 0, 1}]
Plot[10 x, {x, 0, 1}, PlotRange -> {{0, 1}, {0, 1}}]

うわつまんね

そういうなら面白い答えを見せてくれ

もうすぐお正月やね。

ある関数 f[x_]:=
をつくってそれを呼んだときにアンダーフローがおきてしまうのですが
なにか解決方法知りませんか？　アンダーフローがおきたら小数点５位まで丸め込むとかできないですかね？

x (1+y^2/x^2)^(1/2) を　(x^2+y^2)^(1/2) に変形するコマンドはありませんか？

Simplify[x (1 + y^2/x^2)^(1/2), Assumptions -> x > 0]

Simplify[x (y^2 + x^2)^(1/2) /x^(2*(1/2))]

Simplify[x (x^2+y^2)^(1/2)/(x^2)^(1/2),x>0]

>>693
ありがとうございます。
古いバージョンのせいか
Simplify[x (1 + y^2/x^2)^(1/2), x > 0] でうまくいきました。

取り込んだイメージを３次元的に回転できる方法ありませんか？
2次元上ではImageRotateでできるらしいんですが...

座標視点をかえる
アニメーション機能をつかえば

Mac版って動作速い？

動作はメモリとCPUできまるから、おなじCPUをつかっているなら差はでないはずだが
表示その他ですこしかわるがあまり差は出ないと信じる

mathematica home editionは６４ビットバージョン？

PlotなどしてX軸に表示される数値(刻み値)を消したりできるコマンドありますか？

>>702
Plot[x^2, {x, 0, 1}, Ticks -> None]

サンキュ

二年三十七日十五時間。

ココでちょっとしたメッセージや
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★
★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★

小沢先生、頑張って下さい。私は最後まで味方になります。

猫

ご教授下さい！
4次元データのcontourを書きたいのです．
ControurPlot3Dを使えば良さそうなのは分かるのですが，使い方がわかりません．
陰関数f(x,y,z)=0の球面上にg(x,y,z)の大きさでcontourを描きたいんです．

http://tech.ckme.co.jp/img/sphere8cdata_matlab.png
こんな感じです

よろしくお願いします！

Mathematicaで

x J0(x) = 2 J1(x) [J1(x), J0(x)はベッセル関数BesselJ[n,x])

を満たすxを求めるにはどうすればいいですか？

NSolveは使えませんでした。

>>708
FindRoot[x*BesselJ[0, x] == 2*BesselJ[1, x], {x, 5}]
FindRoot[x*BesselJ[0, x] == 2*BesselJ[1, x], {x, 8}]
etc.

>>709
出来ました。ありがとうございます。助かりました。

Solveで実数解だけ取り出すにはどうすればよいでしょう？

たとえば Solve[{3*x^2+y==0,x+3*y^2==0},{x,y}] だと虚数解も出ます
実数解だけ表示したいのです

Solve[{3*x^2+y==0,x+3*y^2==0},{x,y}, Assumptions->{Element[x,Reals]}]
とやってみたけど Solve は Assumption なんてオプションは知らない
と言われました
よろしくお願いします

Solve[{3*x^2 + y == 0, x + 3*y^2 == 0}, {x, y}]でだめなの？

Solve[{3*x^2 + y == 0, x + 3*y^2 == 0}
Select[x/.%,Element[#,Reals]&]

>>714 どうもありがとうございました！

Solve[{3*x^2 + y == 0, x + 3*y^2 == 0} ,{x,y}]
Select[{x,y}/.%,Element[#,Reals]&]

AVIとかの映像ファイル全般がAnimationのインポートができなくて

Import["Capture.avi", "Animation"]

Import::fmterr: データをvideo形式でインポートすることができません． >>

と出てくるんだけど原因に心当たりある方いませんか？
使ってるver6のドキュメントでは当然のように読み込めてるみたいだし

もしくは読み込めるvideo形式に何か制限があるのかもしれないけど

>>52
Maxima
歴史的には、こっちの方がふるい。

機能的としてはmathematicaの圧勝

のはず。あんだけ高いんだから。

>>719
>>あんだけ高いんだから。
Home Editionは３００ドル弱。そんなに高くない。

日本では買えないが、、、。

3次元y軸まわりの回転行列を生成する関数はありませんか？
RotationMatrix3Dだとz-x-zのオイラー角らしく、y軸回りにならないのですが。
z軸90度回転行列で前後をはさめばいいことですが、見難くなるので。

>721
RotationMatrix3Dは無い。

RotationMatrix[x,{0,1,0}]
ではどうか？
だめなら1の位置を変えて試すべし。

>>722
Ver.3.０なので、RotationMatrixはないみたいです。
RotationMatrix3Dで何とかするしかないですね。

ヒント：わざとMathematicaの形で書いていない

平面上の回転行列
{{Cos x,-Sin x},{Sin x,Cos x}}
を(u,v,w)空間の(u,ｗ)だけに作用させる。

あとは自分で考えよう！

>>724
それって自分で関数を書けってことですよね？
引数をとって行列を帰す関数なんて、どうやって書くんでしょう？

>725
行列の積だけなんだけどな、
これで分からなければMathematicaやめた方がいい。

>>725
行列たって、ListのListだろ？それを返す関数なんて別に難しくはないだろ？

>>726
>>727
どうも。単純にこれでよかったんですね。

Ry[s_]:={ {Cos[s], 0, -Sin[s]},{0, 1,0}, {Sin[s], 0,Cos[s]}};

Cだと行列は帰せないから配列を変数で渡さないと、
とか考えてしまって。

mathematicaは思い出したようにたまに使うだけで、
で行列を扱うのがほとんど初めてだったので、
だいぶまごつきましたけど、おかげさまで解決しました。

b-aと入力すると自動的に-a+bになってしまうんですが、
デフォルトの順序を変更するなどして、b-aのままにしておく方法があったら教えて下さい。
やりたいのは、Cのプログラムを出力することで、
-a+bと書くよりb-aと書いた方が速いと思うからです。
イメージ的には
Plus[b,-a, Order->f]
みたいにして項の順序を決定する関数を指定できたら、と思うのですが。

今のCコンパイラーならそれくらい気にしなくても大丈夫じゃないか？
試してみれば？

Mathematicaを勉強しようとおもっているんだよね。でもMapleにしようかと迷ってもいます。
その二つでどちらが使いやすいとおもいますか？

日本だとMathematicaしか使ったことない人が多いと思う。
使いやすさとは関係ないけど、日本語の情報量の多さで言えば断然Mathematicaということに
なる。

だね。メイプルは日本のサイバーネットが子会社にしてしまったんだね。
自分がもっているメイプルは入力がマセマティカみたいにwiyswigでないから使いにくい

新しいバージョンはwysiwigみたいだが

なるほど
メイプルは
だからカナダなんだね

MapleとかMathematicaとかを使えるようになったとして、それがどう数学の勉強に
役にたつとおもいますか？

いつも紙に書いていた式変形が超楽になる

紙に書いた方が、自分のあたまに残るし、予想した通りの式がでたときに
うれしいでしょう？

linuxで使うとなると実質選択肢がmaximaしかなくなるのが痛いが…linuxでは数式処理ソフトどうしてますか？

tex

>>738 linuxでもMathematicaを使ってるぞ

linux版を持ってないんで…

数学の勉強には余り役にたたないと思う。
むしろ科学への応用がメインなのでは。

>>742
英語が読めない旧来型日本人の勉強法だとそうでしょうね。

グラフがホイホイ作れるのがどれほどありがたいか
そのほか思い浮かんだアイデアをほとんどその瞬間で検証できるんだから
勉強が進みますよ

単純計算なんて計算機の仕事
何をどう計算すべきかを考えるのが人間の仕事

計算機にどう計算させるかは人間が与えている訳で、
使っている人間が賢ければ計算機も賢くなる。
単純計算しかできないなんてことは無いと思うよ。

俺は、家計簿の集計から地震の研究までMathematica一本でやってるぜ。

固定したnに関して、

　　すべてのn次（対称）行列Aの対して、性質expr(A)が成立つ

ことを、ForAllを使って調べたいのですが、
「全ての行列」ってどう表現すればよいか
どなたか教えてもらえますか。

3Dのグラフを書きたいのですがy軸のみ対数プロットしたいです。そのようなコマンドは存在しますか？
2DではLogLogPlotやLogLinearPlotは存在するようですが、3Dで見つけることができませんでした。
ver7を使っています。
よろしくお願いいたします。

yの値だけ Log を取れば。

>>748
こんな感じ？「実対称行列の固有値は実数」

n = 2;
matrix = Table[a[i, j], {i, 1, n}, {j, 1, n}];
alist = Flatten[matrix];
expr[x_] := ( Element[Eigenvalues[x], Reals]);
ForAll[Evaluate[alist],
Element[alist, Reals] && matrix == Transpose[matrix],
expr[matrix]] // FullSimplify

n=2までしか解けなかったorz

パラメータ表示をした関数の微分はどうしたらいいでしょうか？

具体的には、
x = t - sin(t)
y = 1 - cos(t)
で書かれたサイクロイドの微分をしたいと思っています
よろしくお願い致します。

>>750
ｙの値をLogで取ってしまうと軸に表示される値はLogを取った値が表示されますよね？
純粋なyの値を軸に表示させることはできるのでしょうか？
やはりyの値だけLogをとらないとだめなのでしょうか・・・

参考書買え

>>753
軸の表示だけなら場当たり的だけどTicksを手で書いてしまう手はある

yのLogを取って、TicksオプションでLogScaleを使う：
Plot3D[,,,,Ticks->{Automatic,LogScale[2,5],Automatic}]
とかはどう？
ここに霊がでてるけど：
http://mmi.tudelft.nl/~wouter/publications/model/texture.html

自分でもテストしたらLogScaleなんてなかったから,
本当に霊かと思ったら、
Needs["Graphics`Graphics`"]
がいるようだ。
でもLogScaleの使い方はいまいちよくわからないので
753が調べてまとめて教えておくれ。

>>751
実質的には
ForAll[{a,b,c},A={{a,c},{c,b}};Eigenvalues[A]//FullSimplify]

だけになる。

実際のＭ演算は実質的にはなにもしていないのでは？

>>752ですが誰か分かりませんか？

>>752
微積分の本を読んだ方が良い、、、。

>>760
dy/dx、d^2y/dx^2、d^3y/dx^3、・・・
と求めたいんですよ
1回微分なら自分でもできますが、高階微分となるとちょっと手計算では・・・

>>761
高階微分でもヤコビアンの微分の項が加わるだけでしょ．
式だけ手計算で出して面倒なところだけ
Mathematicaにやらせればいいんじゃない?

>>761
d/dx=1/x' d/dt
だから、例えば
x[t]=t-Sin[t]
y[t]=1-Cos[t]
d[f_]:=D[f,t]/x'[t]
とでも定義すれば
d^n y/dx^n は Nest[d,y[t],n]
で計算できるけど、これじゃダメなの？

>751
イメージはそんな感じです。もっと簡単に整数の行列で、

　1: ∃matrix={a[i,j]}.∀i.a[i,i]=1…対角成分が1の行列は存在する
　2: ∀matrix={a[i,j]}.∃i.a[i,i]=1…どんな行列も対角成分のどれかが1

を考えるとき、751のalistを使って

　１．Resolve[Exists[alist, Element[alist, Integers], ForAll[i, Element[i, Integers], a[i, i] == 1]]]
　２．Resolve[ForAll[alist, Element[alist, Integers], Exists[i, Element[i, Integers], a[i, i] == 1]]]

のようにすると、1はTrue、2はFalseとなってほしいのですが、
論理式が表示されるだけでTrue/Falseを計算してくれないのが
どうしてなのかわからないのです（ちなみに５版です）。

>>763
ありがとうございます

私の手計算では、2回微分のとき、
d^2y/dx^2 = dt/dx d/dt(dy/dx)
= 1/(1-cost) ・ (cost (1- cost) - sint sint)/(1-cost)^2
= (cost - cos^2 t - sin^2 t)/(1-cost)^3
= (cost - 1)/(1-cost)^3
= -1/(1-cost)^2

という感じになるんですが、
mathematicaにsimplifyさせると
-1/4 (Csc(t/2))^4
になっちゃうんですよね・・・

この両者をTrueQで比較してもうまく行かないんですが、
どうしたらいいでしょうか？

>>765
それくらいなら手でも検算できるが（cost=1-2sin(t/2)^2だから）、
mathematicaでaとbがおなじかどうかは、
a-b//Simplify
として運良く0になればわかる。

>>766
なるほどｗｗｗ
昔からmathematicaの三角関数の変形は融通がきかなくて困りますが
こういうチェック方法も面白いですね

そんなことも知らんで使っていたのか。

>>764
具体的に計算機がどう動くかイメージできないいじょう数式処理ソフトもだめだとおもう。
計算機がどうやればできるかだね

ただ宣言すればできるなら、数学はかわるとおもうけど　裏打ちを支えるのが数学だから相互に発展するんだろうとおもいます。

>>763

Nest「ｄ、ｙ、２」だとｙ’がでてきてそれいじょうすすまないんですけど
どうしたらいいでしょうか？

y[t]=じゃなく、y[t_]=...にした？
x[t]=もx[t_]=だが、、、。

かわりません

f[t_] := {t - Sin[t], 1 - Cos[t]};
D[f[t], {t, 4}]

で4回微分出来たけど

orz

>>773
>>763 をよめ！

>>773のようなしろうとはこまるな

たぶんさ〜、微積分の教科書を読んだ方が早いと思う。

Nest[d, y[t], 2] // Simplify:
-(1/4) Csc[t/2]^4
Nest[d, y[t], 3] // Simplify:
1/8 Csc[t/2]^8 Sin[t]
Nest[d, y[t], 4] // Simplify:
-(1/16) (4 + 3 Cos[t]) Csc[t/2]^10
...
Nest[d, y[t], 10] // Simplify
-(1/4096)(3193801 + 4517832 Cos[t] + 1553157 Cos[2 t] +
233370 Cos[3 t] + 11340 Cos[4 t]) Csc[t/2]^28

とできるけどな。あってるかどうか知らんが、、、。

>>770
d[f_] := D[f, t]/D[x[t], t]
って書けばちゃんと動きましたよ

Simplify[Nest[d, y[t], 2]]とします

あ、そういうことか。

こういう計算をMathematicaで済ますのは感心しないな。

計算するためにはアルゴリズムを考えなきゃいけないし、
アルゴリズムが分かれば自力で計算できる。

この場合Mathematicaは答え合わせとして使うべき。

>>780
5回微分ぐらいしてみてください
関数形がひどいことになるんでｗ

Nest[d,y,2] 770
Nest[d,y[t],2] 779

d[y]=0
d[y[t]]=sin(t)/(1-cos(t))
d[y[z]]=0

ｔを意識しないといけない定義なんだね
一般性にかけるのかな？

>>782
陽にy=f(x)と書けないからね

d[f_,z_]:=D[f,z]/D[x[z],z]

にするとd[y[z],z]=sin(z)/x'(z) になっちゃうなあ

>>784
　x[t]=-->x[t_} (>>773 の指摘）

Nest表現はどうなるの？

眠くて勘違いしてたらごめん。
Solve[x == 1 - Cos[t], t]
out=
{{t -> -ArcCos[1 - x]}, {t -> ArcCos[1 - x]}}

後者を採用すると、
y[x_] := t - Sin[t] /. t -> ArcCos[1 - x]

y[x]
out=
Sqrt[1 - (1 - x)^2] + ArcCos[1 - x]

D[y[x], {x, 5}] // Simplify
out=
(3 (-5 + 12 x - 12 x^2 + 8 x^3))/((-2 + x)^4 x^3 Sqrt[-(-2 + x) x])

でどう？

　すごいですね。
　　
　わたしも形はちがっても差は０になりました。　眠いですね

MatrixPowerがあって、何でMatrixLogがないんじゃボケ。

新関数のSendMail使ってみたけど、どうもSSLによるMail送信がうまくいかないですねえ。
非暗号化ならうまくいくけど。誰かgmailのsmtpでSendMail送信できた人いますか。

Mathematica Home Editionはなんで日本じゃ買えないの？

売ってないから

代理店が妨害してるから。
だれか消費者庁に通報しろ！

入試問題検討委員とかいうのにさせられた(´･ω･`)
高校数学なんてすっかり忘れてる(´･ω･`)
Mathematica様（と皆様）だけが頼りです

Sin[x]+Cos[x] = 1/3 のとき、Sin[x]^3 + Cos[x]^3 を求めよ

なんて、Mathematicaで解ける？

因数分解しろよ

sin(x)^2やsin(x)^3にはsinやcosで1次式に直せるでしょ。

三角関数の展開で出来るはずだよ。

>>794-795
すいません、出先なのでお礼だけ
GW明けに試してみます

えええ？いくらなんでもそれくらい手で計算できるのでは？

こういう問題って
(sin(x) * cos(x) = 〜
の式があると簡単になるよね。

hoge

sx=1/3-sy
sx^3+sy^3=1/27-sy/3+sy^2=1/27-sx/3+sx^2=1/2{(2/27)-(sx+sy)/3+1)}=77/153

↑1
3/27

ユークリッド互叙法の関数Gcd(x,y)を再定義をせずに自作関数で作りたいのですが、ご教授いただけないでしょうか。

gcd[x_, y_] := ({x0 = x, y0 = y};
While[y0 > 0, gcd[x0, y0] = gcd[y0, Mod[x0, y0]]]; Return[x0])

こういう雰囲気でいけないかなと試行錯誤してます。

>>802
グローバル変数x０、y０を使い、関数自身をダイナミックに変更してるようで
ちょっとよくわからないのですが、簡単に、
gcd[x_, 0] := x
gcd[x_, y_] := gcd[y, Mod[x, y]]
ではいけないのですか？

ありがとうございます。
そう書きたいのですが、再定義禁止（言葉が悪くかったです、オーバーロード禁止？）という条件で
やりたいです。

>>804
Fortranの様に手続き的に,しかも再帰呼び出しなしということですか？
だったらこんな感じでどうですか？
gcd[x_, y_] :=
Module[{x0 = x, y0 = y, r},
While[y0 > 0, r = Mod[x0, y0]; x0 = y0; y0 = r]; x0]

もう少し簡単にして、これでも動きます：

gcd[x_, y_] :=
Module[{x0 = x, y0 = y}, While[y0 > 0, {x0, y0} = {y0, Mod[x0, y0]}];
x0]

>>805
>>806

ありがとうございます！

１００以下の正整数のうちpq（p,qは相異なる素数）の形のものからなるリスト

を作りたいのですが教えてくれませんか？

（ヒント：nがpqの形になるための必要十分条件は、nの約数の個数が４、
かつnが立方数ではない（すなわち、n^(1/3)が整数でない）ことである）

ヒントがあってもわかりません。お願いします。

計算機に頼るのならば、
ヒントを使うよりも素直に2〜47の相異なる2つの素数の積のうち
100以下のものをリストアップした方がはやい気がする。
一個一個の数字の約数の個数を数えるのは骨が折れる。

>>808
FactorIntegerで素因素分解して、２つの素数の積でなりたち、指数が１、のを選べばできると思うけど、
そんな解き方でいいんだろうか？
a=Table[FactorInteger[i], {i, 100}];
b=Cases[a,x_/;Length[x]==2];
c=DeleteCases[b, {___, {_, x_ /; x > 1}, ___}]

あ、最後cを元に戻さないといけなかった、、、

NDSolveに詳しい人にお聞きしますが、
（2階線形、２変数）偏微分方程式を長方形の境界を持つ領域（0<x<a,0<y<b)で解きたいとき、
境界条件として、向かい合う辺での値を与えて(u(0,y)=f(y),u(a,y)=g(y))解くようなことはできないでしょうか？

正整数nに対し、2^nの先頭の数字を与える関数をf5[n]を書け。これを用いて2^n（１≦n≦１００）の
戦闘の数字のリストを作れ。またそのうち、１が何個あるか調べよ。

わからないのですみませんが誰か教えてくれませんか？お願いします。

ToString,StringTakeでもできるが、頭を使って解くんじゃないのか？

>>814
頭を使ってもわからないから質問してるんじゃないのかな？
僕は初心者だからもちろんわからん。

>>813
力づくでダサくてもいいなら、

ToStringで数を文字列とし、
StringTakeで文字列の一部をとりだす。
文字列としての数を、数に変換したければ、ToExpressionでできる。

それで所望の関数f5[n]ができたら、1から１００までのリストxに
Map[f5,x]とすれば数字リストが得られる。
１が何個あるかは、Countをつかえばよい。

>>815にもだが、
このような力づくの方法ならそれほど難しくはない。
mathematicaの解説本があるからそれを一度読んだ方が良い。
mathematicaのhelpも自分で読むべきだし。

mathematicaブックって持ってた方がいい？

mathematicaブックってmathematicaのhelp内にはいってない？

Mathematicaはバージョンアップが速すぎてテキストが追いついていない。
しかしヘルプが充実しているので、ヘルプを見ればだいたい分かる。
また、デモンストレーション・プロジェクトを眺めていると、自分のやりたいことのヒントが見つかる。

持っていると便利だし
本の方が画面より見やすいから
古い版のやつでも3000円位で売っていたら買いだろうね。

>>816
で、答えは？本読んでもよくわからん。

Clear[f];
f5[n_] := StringTake[ToString[2^n], 1];
Count[Table[f5[n], {n, 1, 100}], "1"]

Clear[f];
f5[n_] := StringTake[ToString[2^n], 1];
Count[Table[f5[n], {n, 1, 100}], "1"]

Clear[f]; →Clear[f5]; ね。

１００ステップの2次元ランダムウオークを作れ。
ただし、東西南北をそれぞれ、｛１、０｝｛ー１、０｝｛０、ー１｝｛０、１｝で
表すことにする。
（ヒント：りすとdirectotions={{1,0},{-1,0},{0,-1},{0,1}}を作り、1、2、3、4の範囲の
乱数kを発生させ、directotionsのk番目の要素をとるとランダムな方角が得られる）

すいませんが練習だと思って解いてくれませんか？課題が出せないんでお願いしますm(_ _)m

ばかたれ
自分でやれ

こんな？

directotions = {{1, 0}, {-1, 0}, {0, -1}, {0, 1}};
NestList[# + Part[directotions, IntegerPart[Random[ ] 4 + 1]] &, 0, 100]

>>827
RandomChoiceというのもあるよ。
NestList[# + RandomChoice[directotions] &, 0, 100]

http://reference.wolfram.com/mathematica/ref/RandomChoice.html
にきれいな答えが出てるね。

>>828
>>829
どうもです。
なるほど、そんな関数もあるんですね。

整数nの各位の数字の合計を出力する関数f4[n]を書け。
例えば、f4[1234]=10.

わかる人います？教えてくれませんか？

>>831
中学生かよ
自分でやれ

ここで質問は許さん!!

>>832
別にいいじゃん。教えてやれよ。
わからないならともかく。

f4[n_] := Plus @@ (Characters[ToString[n]] // ToExpression)

f4[x_] :=Total[IntegerDigits[x]]

831==834以外には勉強になるな、、、。

835です
8の噂とかを尋ねて、たまたま通りかかったばぃ
>>836
上手いな…おまん若いじゃろ!!そんな関数知らん
かった、ありがとうなぁ
>>831
問題投げかけといて、音沙汰なしかい!
ぼてくりこすぞ〜

次のお題 マダー？

ありがと

>>839や他のみなさんへの課題
てか教えて

正整数nのケタ数を与える関数f2を書け。たとえば、f2[40172]＝５

↑
問題間違えました。

０より大きく１より小さい既約分数で分母がn以下のものを
小さい方から並べてできるリストを出力する関数f5[n]を書け。たとえば、

f5[n]={1/4,1/3,1/2,1/3,1/4}

（このような数列をnに対するFarey数列といい、いろいろおもしろい性質が知られている。
たとえば、隣り合う項をa/b,c/dとすると、必ずab-cd=-1となっている。）

すいませんがお願いします。

Length[Characters[ToString[n]] // ToExpression]
かね。（未検証）

ゴメン>>843は>>841のレスね

>>843
f5[4]をしたら１になりました。違うかもしれません。
すいません。

>>843
間違えました。>>843さん、わざわざすいません。
841は解けました。本当にすみません。

正整数nの各位の数字を逆転させた数をnの逆転数といい、r(n)と書く。
たとえば、f[1234]＝4321

r(n)=n/4となるような最小の正の整数を求めよ。

お願いします。意味わかりません。

>>847
死ね

>>847
意味わからんぞぃ！
例えばのfはrでええんかぁ
んでn/4はリバースする前のnが4の倍数なんか？
4の倍数を利用しないブルートフォース版じゃが
n = 1;
While[n/4 != FromDigits[Reverse[IntegerDigits[n]]], n++];
Print[n]

宿題は自分でやれよ

>>849
ありがとうございます！！！！！！

こういう課題って、色々示唆に富んでいてホント勉強になるわ。
演習問題の書込み→回答は
このスレッドの価値を上げますね。だれも損しない。

あほか

解けないなら黙ってろよ

くだらなすぎて

下らなすぎて解けないってか？w

本に載ってる有名なことばかりだろ
知らないバカ

お前他人が楽しそうに遊んでいると怒りを覚える人間だろw
刃傷沙汰起こす前に病院行けよ

おまんMathematica使いなんか？それともポケットのなかで贋金を鳴らす評論家か？
上の問題はたしかに入門編レベルじゃけん、Mathematica使いには簡単じゃろ！
しかしなぁMathematica使いなら、すて台詞ばかりはいてどうするんかぃ？
そや上のf5[n]について、プログラム板でのマルチ投稿(マルチは気に入らんが…)
があって、そのレスが以下ばぃ！
-----------------------------------------------------
292 ：デフォルトの名無しさん：2010/06/16(水) 15:22:37
>>291

あまり美しくないが。

f5[n_] := DeleteDuplicates[Sort[Flatten[Table[Table[Which[0 <= a/b <= 1, a/b, True, {}], {a, 0, n}], {b, 1, n}]]]]
-----------------------------------------------------
これでもええんじゃけど、思考的にうかぶのは上の形を関数３個で表わした形じゃろ？
おぅ！>>857！！
こんなん入門書０章じゃけんど応用がきかんとだめばぃ！
みんなにわかりやすぅ〜く説明したってくれんかぁ！！

いっとくけんど、くだらんとか簡単すぎて書けんとかなしじゃど！!

おりさまか？おりはココの住人でもないし、Mathematica使いでもないばぃ！

ヒマだから誰か簡単めの学校の課題貼り付けてよ
本屋に行っても「例題で解く〜」くらいしかいい本がないんだよね〜

くだらん

>861
つまらんなぁ…まぁがんがれやぁ
つか、おりもあきたけん

>>847
>>849

ANS=879^n12 n=0....infinity

8712 87912 879912 8799912 .......

おぅ!>>863!!
>>849
と問われたけん、仁義のれすばぃ！
>ANS=879^n12 n=0....infinity
が一般論なんか？んなら突然の879
を説明せんかい！！

つか、ぶぅぶぅ〜

あっ　訂正します。
879^(0)12 は　８７１２と読んでください。　n=1..infinity　です。
8712
87912
879912
8799912
87999912
87999912
879999912
.......

証明もお願いします。

ほう！なるほろ…おもろい
>証明もお願いします。
おりか？むっつりスケベの
>861
にまかす!!
こういう板はむずい…

正の整数リストから元の順序を変えずに重複するものを取り除く関数を書け。
{3,1,1,5,2,1,2}み適応し、{3,1,5,2}になることを確かめよ。
（Unionを使うと{1,2,3,5}になるので注意）

なんか方法ある？

unionのヘルプ見たらそのまんまの答え発見しちゃったよ。。

ここ宿題スレじゃねえよ

またケチだけの糞虫が来たな。
文句しか付けられないのかお前は。
能力か人望のどちらか位は備えていないと人生キツいだろう。

８６８ありがとうございます。

タクシーの番号を見てラマヌジャンが
「その数は2通りの仕方で2つの正の立方数（立方数とはn^3）の和で表される
自然数のうち最小のものです」さてその数は何か？まずf[n]を作ること。

答えはたしか９１になればいいと思います。
僕にとってムズイので教えてください。

追記
91 = 6＾3 + (-5)＾3 = 4＾3 + 3＾3

訂正
「まずf[n]を作ること。」は消して。いらんこと書いた。

>タクシーの番号を見て…
>答えはたしか９１…

おぅ！おまんラマヌジャンしっとるんか〜
つか、タクシーの番号が二桁ってなんかい！
逸話も間違っとるし、-5は自然数じゃなか！

問題にはそう書いてあります。

タクシーの番号を見てラマヌジャンが
「その数は2通りの仕方で2つの正の立方数（立方数とはn^3）の和で表される
自然数のうち最小のものです」さてその数は何か？
とね。

正の立方数で考えてくれませんか？

91じゃないかもしれません。すいません。

逸話は、
1729 = 12^3 + 1^3 = 10^3 + 9^3
んで、あとはかんがえれ

わかりません(>_<)

教えて〜

↑このアホに教えてやれよ。まぁ俺様もアホだが（笑）

例えば，
c:\work\TEST.nbで保存したファイルを開いたときに，
c:\workを取得する簡単な方法がないでしょうか？

Export[]を使うときにフルパスを一々指定しないで済ませたいのですが…

相対パスを指定したら？
..\

タクシーの番号を見てラマヌジャンが
「その数は2通りの仕方で2つの正の立方数（立方数とはn^3）の和で表される
自然数のうち最小のものです」さてその数は何か？

くだらん

ラマヌジャンはもっと高いところにいた
そんなエピソードは忘れろ

はいはい、かまってちゃん。

ラマヌジャンの問題早く教えてください

>>882
NotebookDirectory[]
は？

>>883
>>888
レスありがとうございます．

申し訳ありません．
Mathematicaのバージョンを書き忘れてました．


NotebookDirectory[]はまさに欲していた機能でしたが，
ver6以上でした．

私が使っているのはver5.1なのでまだ実装されていなかったようです．
とりあえずこのまま使い続けることにします．

ラマヌジャンの問題早く教えてください

↑このアホに教えてやれよ。まぁ俺様もアホだが（笑）

NotebookDirectory[]で検索したところ，
ver6以前でも使える方法がいくらか紹介されていました．

ttp://www.mofeel.net/1164-comp-soft-sys-math-mathematica/5574.aspx
ttp://blog.hulinks.co.jp/2010/02/mathematica-gis-2.html

これによると，例えば

SetDirectory["FileName" /. NotebookInformation[EvaluationNotebook[]] // First // ToFileName];
list = Table[{x, Sin[x]}, {x, 0, Pi, 0.1}];
Export["sine.txt", list, "Table"]

で sine.txt を現行ノートブックと同じフォルダに生成することができました．

因みにversionは

Version Number: 5.1.1.0
Platform: Windows

です．ありがとうございました．

君が解答案を提示できたら添削してやる。

いちおう、コメントだけ。。。
t=p^3+q^3=r^3+s^3
なら、p^3+q^3-(r^3+s^3)=0
を満たす四組が必要条件じゃろ！

んでブルートフォース ばぃ！！
つか、この手の問題ってそんなに
綺麗に解ける問題なんか？

つか、
Select[PowersRepresentations[#, 2, 3] & /@ Range[20000],
Length[#] == 2 &]
{{{1, 12}, {9, 10}}, {{2, 16}, {9, 15}}, {{2, 24}, {18, 20}}}

Map[Total, r^3, {2}]
{{1729, 1729}, {4104, 4104}, {13832, 13832}}
で、
1^3+12^3=9^3+10^3=1729
なんか？

6より大きな偶数nは二つの素数p,qの和で書けると予想できている。（ゴールドバッハ予想）
6より大きな偶数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はどこで入れるの？最初から違う？

ex.

goldbach[n_] :=
Select[Table[{i, n - i}, {i, 3, n/2}], (PrimeQ[#] == {True, True}) &]

↑
すごいね。
ありがとう。マジで感謝！！

正整数nが対称数とは、nの逆転数とnが等しいものである。
例えば１２３２１は対称数である。
正整数nに対し、それが対称数ならnを出力し、
そうでなければ、nとその逆転数を加え、それが対称数かどうか調べる。
この操作を対称数が得られるまで繰り返し、
得られた対称数を出力するような関数f[n]を書け。
たとえば、
78→78+87=165→165+561=726→726+627=1353→1353+3531=4884
なのでf[78]=4884である。

わかんないです。助けてください。

ex.

f[n_] := Module[{r = FromDigits[Reverse[IntegerDigits[n]]]},
If[r == n, Print[n], f[n + r]]]

すいません。本当にありがとうございます。

De Polignacは「5以上の任意の奇数は(2^m)+p(m>=1,pは素数)と書ける」と予想した。
この予想が正しいか、実験して確かめよ。

ごめんなさい。お願いします。

n = 5;
While[Apply[Or, PrimeQ[n - 2^Range[Floor[Log[2, n]]]]] == True,
n += 2];
Print[n]
127

なめとんのかぁ？ぜんぜんだめじゃん！

問題がわかりにくかったかもしれません。もう1度書きます。

De Polignacは「5以上の任意の奇数は(2^m)+p
(m>=1,pは素数)と書ける」と予想した。
この予想が正しいか、実験して確かめよ。

>>904
だからn=127で既に予想が崩れてるって...

予想がただしけりゃ、
n - 2^Range[Floor[Log[2, n]]] /. n -> 127
のリストには素数が少なくとも一つはあるはず
じゃろ？しかしこれがFalseでWhileが終了...

そっかー、ありがとう

n人のクラスの中に誕生日が同じ人が1組以上いる確率を求める
プログラムbirthday[n]を書け。
birthday[50]を1000回試行して平均確率を求めよ。

確率苦手です。

漸近値:1-(365P50)/365^50=0.97

birthday[n_] :=
Table[Random[Integer, {1, 365}], {n}] // (Length[Union[#]] != Length[#]) &

Table[birthday[50], {1000}] // (Count[#, True]/Length[#] &) // N

プロは凄いな。尊敬します。

確率の授業で必ず習うじゃん

Mathematica programmingを？
確率の授業で必ず習うってことはない

n=3^k+a1*3^k-1+a2*3^k-2+・・・・+ak-1*3+ak

と表される。ここにa1,a2,....,akは-1，0，1のいずれかである。
たとえば25=3~3+0*3^2+(-1)*3+1
nに大志、上のようなリスト{a1,a2,....,ak}を与えるような関数を書け。
25に対してはリスト{0,-1,1}を与えるはずである。
またn=1000に対して実行してみよ。

お願いします。

>>912

n = 25; lis = {};
While[n > 0, m = Mod[n, 3];
lis = Prepend[lis, If[m != 2, m, -1]];
n = If[m != 2, Quotient[n, 3], Quotient[n, 3] + 1]];
Drop[lis, 1]

でどうですか？n=1000でも試してみて下さい

手抜き版
nにm=3^(k+1)+3^k...+1を加算すると、akは{-1，0，1}+1->{0,1,2}で、
一般的な3進数のためIntegerDigits[#,3]で係数を取得する。つぎにmを
引いて(各係数から1を引く結果akの範囲は{-1，0，1})nに戻したあと、
//〜で整形処理ばぃ

ex.
f[n_] := IntegerDigits[n + (3^(Floor[Log[3, n]] + 2) - 1)/2, 3] - 1 //
If[First[#] == 0, Rest[#], #] & // Rest

>>912

n = 10000
Clear[aa, bb, cc];
aa = IntegerDigits[n, 3];
k = Length[aa] + 1; bb = Table[0, {k}]; cc = bb; Do[
bb[[k - i]] = aa[[i]], {i, Length[aa]}]
Do[Which[bb[[i]] == 2, {bb[[i + 1]] += 1; bb[[i]] = -1},
bb[[i]] == 3, {bb[[i + 1]] += 1; bb[[i]] = 0}], {i, k}]
Do[cc[[k - i + 1]] = bb[[i]], {i, k}]

n=10000 のとき
{1, 1, 1, 2, 0, 1, 1, 0, 1}ーー＞{1, -1, -1, -1, -1, 0, 1, 1, 0, 1}

になりますね。
>>914　さんの方法の初等的なやりかたですが、（頭のいい人はちがうなあああ）

これはＡｎｄｒｅ　Ｗｅｉｌ　の整数論初歩の練習問題（証明せよという形式）に出てきますね。

Ａｎｄｒｅ　Ｗｅｉｌ　の問題は
　天秤と｛１，３、。。。、３＾n-1}ｇの錘群をつかえばすべてのＮｇが計量できることを証明せよ

のようなものでした。
シロウト向けの講義いい加減な記憶ですから、責任はとりません。

>>913
n=2345 {0, 1, -1, 0, 0, -1, -1}
>>915
n=2345 {0, 1, 0, 1, -1, 0, 0, -1, -1}

アラ違う？！

>>914
ちょっと一本取られました。綺麗な解答ですね。

>>917
914,915では正解でした。検算
In:=
3^7 + 3^5 - 3^4 - 3 - 1
Out:=
2345

すみません。正解は913,914です。

ありがとうございます

>>931をヒントの再帰版ばぃ

f[n_, l_: {}] :=
If[n == 0, Rest[l],
f[Quotient[n + 1, 3], Prepend[l, Mod[n + 1, 3] - 1]]]

f[25]
{0, -1, 1}

どや？

>>931でなく>>913だったorz

マセマティカによりモンテカルロ法を用いて、何か平面図形の面積を求めてください。
そして、その数値と求積の公式や積分を利用して得られる理論値との違いを調べてください。
すなわち答えることは以下の4点です。
（１）どうゆう図形について考えましたか？
（２）その図形の面積をモンテカルロ法で求めるとどうゆう結果がでましたか？
（３）その図形の面積を数学で求めるとどうなりますか？
（４）あなたがモンテカルロ法で計算した結果は理論値に対して何％のずれがありますか？

すいませんがおねがいします。

>>921

９２１のやりかた　ｆ「２３４５」＝{0, 1, -1, 0, 0, -1, -1}
９１７のやりかた　ｃｃ＝　ｆ「２３４５」＝ {0, 1, 0, 1, -1, 0, 0, -1, -1}

で９２１はあやまり？（けたくだりについて考え不足でした）

９１７は再帰ではなく、Ｃでもかけます。　めもり、計算速度の効率がいい。　でもつまんない＞＜；；

>>923　はべつにＣでもＢａｓｉｃ　でもいいんじゃない。　とりあえず円でも計算したら。
乱数の性質と精度と計算時間とビット長の関係でも報告したら、いいんじゃない。

>で９２１はあやまり？

ん？まちがっとるかぁ？
2345 = {1x3^7+0x3^6+1x3^5-1x3^4+0x3^3+0x3^2-1x3^1-1}
-> {1, 0, 1, -1, 0, 0, -1, -1} == {1, a1, a2, ... ak}

で、{a1, a2, ... ak}の{0, 1, -1, 0, 0, -1, -1}を返す
のとちがうんかぁ

＞{a1, a2, ... ak}の{0, 1, -1, 0, 0, -1, -1}を返す
のとちがうんかぁ

それはずるい
はなしならん

>>925
マセじゃないといけないらしいです。
モンテカルロ法が今ひとつわからないです。お願いします。

円の1/4の面積＝Pi/4 を計算する。

ss[x_, y_] = x^2 + y^2 - 1;
pp := (mx = Random[]; my = Random[]; ss[mx, my]);
tt[n_] := Table[pp, {n}]
Area[n_] := Count[Negative[tt[n]], True]/n // N
Truearea = Pi/4
ErrorDiff[n_] := Area[n] - Truearea

Table[ErrorDiff[n],{m,1000,10000,1000}]

で検討してみてください

>>918　ちょっと一本取られました。綺麗な解答ですね。

はは　畳の上の水練ですな　クリティカルパスの設計で修羅場をくぐっていませんね

>>930

クリティカルパスの設計そのものが全く分かりません。
プログラムは昔BASICを自己流でやっただけなので、ヘルプを見ながら何とか作っている状態です。

でも、このスレはとても勉強になるので感謝しています。

>>929
すげー。ありがとう。Ｓランに改名しろよ。

>>927
あんなぁ、おりが>>914で『手抜き版』と謳ったけん、
921で>>913さんの純正の正面突破ロジックを、
おりなりに圧縮&再帰化した『再帰版』じゃけど！
>>925の疑義がでたけん、簡単な説明をしたばぃ！
んなら…
-------------------------------------------------
＞{a1, a2, ... ak}の{0, 1, -1, 0, 0, -1, -1}を返す
のとちがうんかぁ

それはずるい
はなしならん
-------------------------------------------------
って、どこがずるくて、はなしならんのか説明せんかい！

なんかソロバン塾で、計算もしないでご高説を垂れるお方
が沸いてますが、徒な無意味エントロピー増大お疲れさま。

1段目の辺の長さnの正方形、2段目がn-1の正方形、...、
n段目が辺の長さが１の正方形であるようなピラミッドがある。
各段は長さ1の正方形のレンガが敷き詰めてある。
レンガの総計が平方数であるようなピラミッドは何段か？

お願いします。

しらん

Mathematicaのｽﾚとして趣旨ずれてきてないかｗｗ

いい加減算数スレに行けよ荒らしども

>>934
Mathematicaと関係ない、やぽー知恵袋で聞けば？

934です
でもこの問題をＭでやれと言われました。

>>939
っMatlab

1+2^2+3^2+....+n^2=n(n+1)(2n+1)/6

n=24 のとき　７０ｘ７０になる

Ｓランさんありがとう

ポーカー
トランプのカードを
{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点

お願いします。

は？

できるかボケェ

わかるかたがいらしたらお願いします

しらん

>>943
まずは「ポーカー　確率」で検索して、自分で考えて下さい。

見てもわかんないです

>>950
そうですか。では「教えて君」を検索して、自分を考えて下さい。

>>951
教えてくれないのですか？わかるなら教えてください。

Fランが教えてくれるよ。

>>953
　　すみません。
この問題は非常によい問題です。（教育の観点から）
問題の記述と解決の努力のなかでいろいろ学ぶところの多いとおもいます。　出題者の先生もそれを期待しているのでしょう。
時間がかかっても自分で試行錯誤してとかれることを進めます。
ポーカをあまりよくしませんが、ルールの表現と判定などいろんなやり方があると思います。

時間がかかっても自分で試行錯誤<-大事
とっかかりの参考例…某ソースの流用版！
carddeck =
Flatten[Map[({{H, #}, {S, #}, {D, #}, {C, #}}) &, Range[13]], 1];
f := RandomSample[carddeck, 5]
EvaluatePokerHand[a_List] := With[{
dist = Sort /@ Map[Length, Split /@ Sort /@ Transpose[a], {2}],
straight =
MemberQ[Sort /@ Take[Partition[Range[13], 5, 1, 1], 10],
Sort[Last[Transpose[a]]]]
},
If[straight,
If[ dist[[1]] == {5},
If[Sort[Last[Transpose[a]]] == {1, 10, 11, 12, 13},
100(*"royal flush"*), 50(*"straight flush"*)], 20 (*"straight"*)],
If[ dist[[1]] == {5}, 25(*"flush"*),
Switch[dist[[2]],
{1, 4}, 40(*"four of a kind"*), {2, 3},
30(*"full house"*), {1, 1, 3},
15(*"three of a kind"*), {1, 2, 2},
10(*"two pair"*), {1, 1, 1, 2}, 5(*"pair"*), {1, 1, 1, 1, 1},
0(*"nothing"*)
]]]]
(EvaluatePokerHand /@ Table[f, {1000}]) // Mean // N

ありがとうございます

あーあ、学ぶチャンスを失ってしまいましたね。

n=10a+bに対して、N=n+100とおきます。
Nの倍数ばかりでてくる数列を作成してください。
そして、その数列に登場する数がすべてのNの倍数に
なっていることを証明してください。
ただし、作成する数列としては、ak=NkのようにNの倍数を含むものや、
連続するN数の積をとるといったような、
明らかにそれとわかるものは除外します。

すいませんがわかる方教えてくださいますようお願いします。