〓 Mathematica 3 〓
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
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/
0643１３２人目の素数さん
2007/06/12(火) 21:27:03元の数と入れ替えたのもとの数がなにかもわからないし、まず、説明方法を考えたほうがいいと
おもうよ。
例えば、関数をつくりたいなら、インプット要素と期待するアウトプットを書くとかね。
0644632
2007/06/12(火) 21:34:40任意の自然数において成り立つ法則だと思うのですが、例えば、
元の数 12
入れ替えた数 21
この2数の和 → 33（11の倍数）
この2数の差 → 9 （9の倍数）
これは2桁の数の例ですが、3桁、4桁、n桁と拡大したときに、
規則的な倍数変化をするように思ったのです。
0645632
2007/06/12(火) 21:41:23インプットは任意のn桁の自然数です。
アウトプット二パターンあり、一パターン目は、そのn桁の自然数を
入れ替えた数との差のが何の倍数になっているか調べることです。
（99の倍数など）
二パターン目は、n桁の自然数と入れ替えた数との和が何の倍数
になっているか結果が帰ってくることです。
0646１３２人目の素数さん
2007/06/12(火) 21:42:10あいかわらずわかりづらいです。
インプットは、N桁の整数
アウトプットは、なに？
あと、33(11の倍数)の3ではなく11のほうを選ぶことも必要なの？その場合、その基準は？
0647１３２人目の素数さん
2007/06/12(火) 21:43:22ニアミス。
アウトプットはわかりました。
何の倍数ってのは、素因数分解じゃなくて、特定のものをかえすひつようがあるの？
0648632
2007/06/12(火) 21:47:07何度もすみません。
>あと、33(11の倍数)の3ではなく11のほうを選ぶことも必要なの？その場合、その基準は？
例えば、二桁の数を 10a+b とおくと、入れ替えた数は 10b+a となります。
その和は 11a+11b なので、11の倍数です。
その差は 9a-9b なので、9の倍数です。
このような11の倍数であるとか、9の倍数であるといったアウトプットが欲しいのです。
0650１３２人目の素数さん
2007/06/12(火) 21:57:55都合のいい例ばかりで一般的な規則がわかりません。
たとえば1392の場合、どんな値をどういう理屈で返すのでしょう？
0651１３２人目の素数さん
2007/06/12(火) 22:05:580652１３２人目の素数さん
2007/06/13(水) 18:58:46func[r_] := (
s = func2[r]; (*func2はなんらかの複雑な処理*)
t /. FindRoot[Sin[t] + t == s, {t, 0}]
)
なる関数を定義したとき、逆関数的に
FindRoot[func[r] == 6, {r, 0}];
などとして、返り値tが目的の値になるようなrを求めようとすると
エラーになります。
なにがいけないのでしょうか？
0653１３２人目の素数さん
2007/06/14(木) 00:42:350654652
2007/06/14(木) 10:28:17func[r_] := (
t /. FindRoot[Sin[t] + t == r, {t, 0}]
)
FindRoot[func[r] == 6, {r, 0}];
このコードはエラーになります。
（ただしエラーと共に正しい答えも…？）
ただこのコードだと、Sin[t]+tにt=6入れれば欲しい答えが出ますが、
そうではないという意味でさっきはfunc2を挟みました。
実際にはfuncはrを受け取ってゴニョゴニョしてループとか条件分岐とか
あって、それからFindRootやって、出てきた答えをまたゴニョゴニョして
返す、といったちょっと複雑な関数なため、逆に辿るのは無理です。
つまり、「ちょっと複雑めな関数を外から逆関数的に使いたくて
FindRootしてみたら、その関数の中にあるFindRootでエラーになる。
なんとかならんか」ということです。
0655１３２人目の素数さん
2007/06/14(木) 15:26:060656652
2007/06/15(金) 23:56:42上の例はwell-definedではないのでしょうか
実際のfuncの中身は、とある光学系における幾何光学的な光線追跡で、
入射光角度を引数に最終的な到達点の座標を返すものです。
(mathematicaにやらせる仕事じゃないかもしれませんが)
逆にここに到達させるためにはどんな角度で…ってやろうとしたら
行き詰ったわけで。
0657１３２人目の素数さん
2007/06/16(土) 01:48:34FindRootで逆関数をもとめられてるわけじゃないんだから。
0659１３２人目の素数さん
2007/06/21(木) 22:25:09皆さんはどうですか？
0660１３２人目の素数さん
2007/06/21(木) 23:28:310661１３２人目の素数さん
2007/06/26(火) 21:23:150662１３２人目の素数さん
2007/06/29(金) 02:53:490663１３２人目の素数さん
2007/07/01(日) 00:45:02指数部が -3<#<3 の時はそのまま
それ以外のときは3の倍数(EngineeringForm形式)
をデフォルトに設定したいのですが，
ヘルプを見てもどこに設定していいやら分かりません。
(オプション Number/Format らしいのですが…)
式は
NumberForm[%, ExponentFunction -> (If[-3 < # < 3, Null, 3 Quotient[#, 3]] &)]
で合っていると思うのですが
設定するオプションの位置をご存知の方
どうぞ宜しくお願いします。
(未だにこの独特のオプション設定ウィンドウの扱い方が分からない…)
0664１３２人目の素数さん
2007/07/01(日) 22:11:21漸化式 a_0=0,a_1=1,a[n+1]=a[n]+a[n-1] で与えられた数列の一般公を求めたいのですが・・・
RSolve[{a[n+1]==a[n]+a[n-1],a[0]==0,a[1]==1},a[n],n]
として shift+enter 押してもできない・・・
0665１３２人目の素数さん
2007/07/01(日) 22:12:330666１３２人目の素数さん
2007/07/01(日) 22:17:30フィボナッチ数って知ってる??
0667１３２人目の素数さん
2007/07/01(日) 22:46:050668１３２人目の素数さん
2007/07/02(月) 02:17:30Mapleがそれに続くだろう。
0669１３２人目の素数さん
2007/07/02(月) 02:48:51それ実行する前にClear[a]とか?
0670習い始め
2007/07/02(月) 18:47:22RSolve[{a[n+1]==a[n]+a[n-1],a[0]==0,a[1]==1},a[n],n] //FullSimplify
でOK
0671１３２人目の素数さん
2007/07/03(火) 20:53:332DPlotで横軸に取るプロット数を任意に決める方法はありますか？
たとえば、Plot[Sin[x],{x,-1,1}]において、プロット数を10個に指定して描くことはできますか？
ListPlot以外の方法でどうかお願いします！
0672１３２人目の素数さん
2007/07/03(火) 21:08:54無理。
ListPlotを使って、Plotと同じ引数をもつ関数かけばいいでＦＡ。
0673671
2007/07/03(火) 21:36:29分かりました、有り難うございます。
ところで、こちらがnbファイルをアップして、プログラムの無駄な部分を検証して下さる方はいますか？
計算時間を短くしたいのです。
無理承知で一応。
0674671
2007/07/03(火) 21:38:29用途は大学院の論文作成です。自分の能力からみても営利的な働きはありません。
0675１３２人目の素数さん
2007/07/04(水) 00:56:17とか？曲率を保とうとするので、xが等間隔にはならないけど。
0676１３２人目の素数さん
2007/07/04(水) 11:49:330677１３２人目の素数さん
2007/07/04(水) 11:52:56一目見てできる作業にMathematica使う必要ないでしょ。
0678１３２人目の素数さん
2007/07/04(水) 12:06:170679１３２人目の素数さん
2007/07/04(水) 12:07:500680１３２人目の素数さん
2007/07/04(水) 12:16:10とりあえずやってみます。
ありがとうございました
0681１３２人目の素数さん
2007/07/04(水) 12:21:480682１３２人目の素数さん
2007/07/04(水) 12:32:55試してみます
0683１３２人目の素数さん
2007/07/04(水) 21:08:140684１３２人目の素数さん
2007/07/05(木) 02:11:21やってみますってだけで結果の報告無しとかいうやつが多いね。
0685１３２人目の素数さん
2007/07/05(木) 08:13:110686１３２人目の素数さん
2007/07/10(火) 15:35:150687１３２人目の素数さん
2007/07/10(火) 16:28:300688１３２人目の素数さん
2007/07/10(火) 16:42:01MM[n_] := MatrixForm[Table[1, {i, n}, {j, n}] - IdentityMatrix[n]]
0689１３２人目の素数さん
2007/07/10(火) 22:17:45できました！
ありがとうございました。
Mathematicaプログラミングを勉強したいのですが、
オススメの本はありますか？
0690１３２人目の素数さん
2007/07/11(水) 02:10:481 - IdentityMatrix[n]
0691687
2007/07/11(水) 13:04:05ありがとうございます。
こういう書き方もできるのですね。覚えておきます！
でも、実は全成分１ではなく文字列の行列で、
正方行列ではなく、m×ｎの行列を求めたかったんです。
分かり易くシンプルに質問させて頂きました。
IdentityMatrixで作ったｎ×ｎの単位行列とm×ｎの行列の演算ができないので、
ii成分だけ0にするように書き換えてみました。
ちょっと手間かかり過ぎですかね？
0692１３２人目の素数さん
2007/07/12(木) 12:24:35それで期待通りの動作をしてるんならいいんじゃないでしょうか。
ていうか何か聞きたいならコード出すなりもっと具体的に書かないと
めちゃくちゃ親切でヒマな人がいたとしてもどうしようもない。
0693１３２人目の素数さん
2007/07/12(木) 20:49:270694１３２人目の素数さん
2007/07/13(金) 03:07:11これ使うと高校数学なんかは余裕なの？例えば整数問題苦手なんだけどこれだと問題文
入力するだけで答えが出るとか？
0695１３２人目の素数さん
2007/07/13(金) 09:40:46苦手なものを苦手なまま回答だけ手に入れようと考えているようではダメだ！
ぶっちゃけ、算盤や電卓より高度ではあるが所詮は道具だから
理解をするのは人間の仕事
0696１３２人目の素数さん
2007/07/13(金) 11:02:480697１３２人目の素数さん
2007/07/13(金) 12:36:160698694
2007/07/14(土) 09:15:30参考までにこれ使ってabcd<a+b+c+d a,b,c,d整数ですべて求めよ。みたいな問題解く
ことってできるの？取りあえず2次方程式はだいたい解けるらしいけどこういうのは
どうなのかなって思って。
0699１３２人目の素数さん
2007/07/14(土) 09:19:160700１３２人目の素数さん
2007/07/14(土) 12:16:06感覚的には， g[f[x]_]:=(h[x_]:=f[x])
という感じで・・・ 実際これを使った時
入力：
g[x^2+3];
Print[h[x]]
Print[h[1]]
結果：
x^2+3
4
を得たいということです
0701１３２人目の素数さん
2007/07/14(土) 12:23:59その例だけだと、定義の中で関数定義なんて発想にならないと思うので、ちゃんと考えて書いてみて。
0702１３２人目の素数さん
2007/07/14(土) 12:31:26これだとうまくいかないということです
g[s_] := (h[x_] := s);
g[x^2 + 3]
Print[h[x]]
Print[h[1]]
結果
3+x^2
3+x^2
となって3が代入されたことにならない
定義の中で定義したいわけです
0703１３２人目の素数さん
2007/07/14(土) 12:34:00訂正：1が代入されたことにならない
0704１３２人目の素数さん
2007/07/14(土) 12:48:10ぽかーん。
0705700,702
2007/07/14(土) 19:08:240706１３２人目の素数さん
2007/07/14(土) 19:19:450707１３２人目の素数さん
2007/07/15(日) 00:20:11自己解決した時も、どう解決したのかは書いてくださいよ
ネットで質問する時の基本事項ですよ
>698
それは余裕で。
0708698
2007/07/15(日) 00:49:59で解けるようですが上の問題もそうなんですか？それとも別のコマンド？
この分ならチャート式くらいの問題は全部コマンドひとつで解けそうで興奮し
てます。
0709１３２人目の素数さん
2007/07/15(日) 16:37:200710１３２人目の素数さん
2007/07/16(月) 00:43:210711１３２人目の素数さん
2007/07/17(火) 15:19:56自分も5.2でいい。6はよくわからん。
0712１３２人目の素数さん
2007/07/18(水) 00:09:43の実行結果
Mathematica 5.0 Mac 版：2
Mathematica 5.0 Win 版：2 + Abs'[0]/2
Mathematica 5.1 Win 版：5/2
厳密解：5/2
これは Mathematica 5.0 以前のバグ？仕様？
0713１３２人目の素数さん
2007/07/18(水) 12:27:58これはまたWebで登録しなきゃいけないの？
0714１３２人目の素数さん
2007/07/18(水) 20:39:33アップグレードも馬鹿高いし、数学は基本的に完成された学問というイメージが強くて劇的な機能強化があるようにも思えず、
非数学分野の自分としては十分に満足して使っていました。
こんな自分は間違っているでしょうか？
さっさと最新版にアップグレード汁！という方がおられたらご意見をください。
0715１３２人目の素数さん
2007/07/18(水) 21:14:010716１３２人目の素数さん
2007/07/18(水) 22:15:480717１３２人目の素数さん
2007/07/18(水) 22:33:44あほらし。正規ユーザーなめてる。
0718◇宮城200か12-54＠Linuxザウルス : ◆mzpjBMNq3U
2007/07/18(水) 23:03:20著者： 竹内薫 発行年月日：2007/07/20 サイズ：新書判 ページ数：198 ISBN：978-4-06-257560-7 価格（税込）：1,260円 近日販売
内容紹介 理系学生なら、絶対必要！ 必ず役立つ！数値計算はもちろん、数式計算やグラフィックまでなんでもこなす、理系学生必携の「数式処理ソフト」。
科学ライターの著者が初心者に向けて、使い方と使いこなしのテクニックを楽しく解説します。「ペンローズのツイスター」にもチャレンジしよう。
高価な市販品に匹敵する数式処理ソフト「Maxima」付き。（Windows Vista，Windows XP対応）
複雑な数値計算から、式の展開や方程式を解く数式処理、さらに計算した値をCGで表示するグラフィックスなどの能力を備えた「数式処理ソフト」は、理系の学生なら是非もっていたいソフトウェアです。
しかし、市販品の数式処理ソフトは大変高価で、簡単に使うことはできません。本書は科学ライターの竹内氏が、Maximaを使って数式処理ソフトの使い方と使いこなしを、分かり易く解説します。
なお、掲載した全てのコードに対応するMathematicaのコードも収録してあるので、Mathematicaのユーザーにも楽しんでいただけます。
目次
1．「数式処理」って何なの？
2．Maximaの概要
3．まずは遊んでみよう
4．コーネル大学の算数問題
5．ファインマン先生の予備知識
6．ペンローズのツイスター
7．超ひもが棲む世界？
8．ボームの量子ポテンシャルを「見る」
9．いろいろな数学の問題に挑戦する
0719１３２人目の素数さん
2007/07/19(木) 02:24:21サポに連絡したらすぐ更新できんじゃないの?
マザボ交換でライセンスの更新が必要になったとき
メールしたらすぐに再発行されたよ。
英語版でWolframの直接サポートだけど。
0720１３２人目の素数さん
2007/07/20(金) 16:59:50License Certificateってパスワード書いてなくね？
0721１３２人目の素数さん
2007/07/20(金) 18:28:46例えば
(1,2),(2,4),(4,3),(7,1),(8,2),(9,4),(11,2)
見たいなやつ
0722１３２人目の素数さん
2007/07/20(金) 18:38:14初歩的な質問ですけど
院生です
0723１３２人目の素数さん
2007/07/20(金) 19:25:44理系大学ならソフトがあると思う。
0724１３２人目の素数さん
2007/07/21(土) 00:30:08生協でアカデミック版が安く売っているだろ
0725１３２人目の素数さん
2007/07/21(土) 00:59:54全然関係ないが、何度見ても「（品物の名）が売っている」には慣れる事ができない・・・
0726１３２人目の素数さん
2007/07/23(月) 08:12:32よくy軸が変なところになってしまいます。原点から
表示してもらいたいのですが、どうすればよいでしょうか？
どうかよろしくお願いします。
0727１３２人目の素数さん
2007/07/23(月) 10:09:27どういうオプションがあるかいろいろ書いてある。
0728１３２人目の素数さん
2007/07/23(月) 14:09:15いま見てみました.
PlotRange->{0,100}
などとすると、できそうです。
ありがとうございました。
0729１３２人目の素数さん
2007/07/23(月) 14:12:23AxesOriginのヘルプも見とき。
0730１３２人目の素数さん
2007/07/23(月) 15:02:34パスワードが通じません。
0731１３２人目の素数さん
2007/07/23(月) 18:53:480732１３２人目の素数さん
2007/07/24(火) 03:06:42ありがとうです。
0733ageage
2007/07/24(火) 08:54:480734１３２人目の素数さん
2007/07/24(火) 08:56:420735ageage
2007/07/24(火) 08:57:440736１３２人目の素数さん
2007/07/24(火) 10:46:43とかやっても、２００桁表示されません。
どうやれば２００桁表示してくれるのでしょうか？
0737１３２人目の素数さん
2007/07/24(火) 12:02:180738１３２人目の素数さん
2007/07/24(火) 12:05:360739１３２人目の素数さん
2007/07/24(火) 21:54:50初回限定版特典のドラマCDが最高だった
0740１３２人目の素数さん
2007/07/24(火) 21:57:260741１３２人目の素数さん
2007/07/25(水) 10:16:475/100にしたら表示されました。
ありがとうございました。
0742１３２人目の素数さん
2007/07/27(金) 02:34:47それはスゴイ
■ このスレッドは過去ログ倉庫に格納されています