〓Mathematica〓

[表示 : 全て最新50 1-99 101- 201- 301- 401- 501- 601- 701- 801- 901- 1001- 2chのread.cgiへ]
Update time : 12/19 02:05 / Filesize : 264 KB / Number-of Response : 1002
[このスレッドの書き込みを削除する]
[＋板最近立ったスレ＆熱いスレ一覧 : ＋板最近立ったスレ／記者別一覧] [類似スレッド一覧]


↑キャッシュ検索、類似スレ動作を修正しました、ご迷惑をお掛けしました

〓Mathematica〓

1 名前：１３２人目の素数さん [2001/03/20(火) 02:34]: ないので立てました。
Mathematicaについて語るスレです。

大学で初めて出会った時は感動したなぁ。(NeXT版…)
Plot3Dでよく遊んだものです。
四次方程式の一般解を出力して、壁に貼ったりしてました（ｗ
348 名前：１３２人目の素数さん mailto:sage [02/12/19 22:15]: mathmaticaって高スペックじゃないとダメかな？？
おせーてください
349 名前：１３２人目の素数さん mailto:sage [02/12/20 03:26]: >>348
計算量による。
職場じゃPenII&384MBで使ってるけど、
MALTAB自体が重いことはない。
350 名前：349 [02/12/20 03:32]: MATLABじゃねぇ、mathmaticaだ
鬱
351 名前：１３２人目の素数さん mailto:sage [02/12/20 10:17]: >348
「ダメ」というのがどれ以上がダメなのか？

こういう計算が何秒以内ならＯＫでそれ以上ならダメとかいう指標が
ないと主観の問題だからね

高性能すぎて困ることはないから今だとPentium4 3.06GHz(Hyper threding)＆メモリ１ＧＢ
くらいにしとけば誰も文句は言わない

Mathematicaってマルチスレッドの恩恵に預かる場面ってあったっけ
352 名前：１３２人目の素数さん [02/12/22 10:54]: >>349

職場でmathematica使ってるなんて凄いね。
俺なんか、職場ではexcel（分析ツールなし）とwindows標準の電卓（普通の方）
で事足りる。技術職だけど。
353 名前：１３２人目の素数さん [02/12/22 11:19]: >352
ＭＡＬＴＡＢ（ＭＡＴＬＡＢの綴りまちがい）だって言ってるんだよボクェ！
354 名前：山崎渉 mailto:（＾＾）sage [03/01/11 12:31]: （＾＾）
355 名前：１３２人目の素数さん [03/01/14 22:18]: すいません。質問なのですが研究室で
mathematicaで波動関数のアニメーションを作ろうと思うのですが現在の環境
clerone 1.3G
memory 512MB
VGA G550 (1024x768)
windows 2000 SP3
液晶ディスプレイ
の環境で動かすとコマ数が64にすると画面やけのような現象が起こり、
ついにはハングアップしてしまいます。
それで充分な描画表示能力のあるハードを買おうと思います。
理想としては256コマを裁けるような感じにしたいのです。大体どのくらいのスペックが必要なのでしょうか?
あまりコストがかかるようならば、ある程度は妥協したいと思ってます。
宜しくお願いします。
356 名前：１３２人目の素数さん mailto:sage [03/01/20 04:35]
357 名前：１３２人目の素数さん [03/01/20 10:46]: (･∀･)ｹﾞﾊﾊﾊﾊ
358 名前：１３２人目の素数さん mailto:sage [03/02/02 09:04]
359 名前：教えて君 [03/02/02 18:26]: リストの操作で　mogi={23,43,59,66,33,25,56,89,51,45.66}
からi番目の要素からｊ番目の要素を足すという操作を自動で行いたい時は、どうしたら
いいのですか。
数列ならシグマをつかえばいいのですが、それに似たことです
360 名前：１３２人目の素数さん mailto:sage [03/02/02 18:37]: Sum[mogi[[k]], {k, i, j}]
361 名前：１３２人目の素数さん [03/02/02 23:35]: マスキャドを電卓がわりに使ってます。

後、メイプル（微分方程式の解析解なんか求めるとき）

mathematicaはつかいずらい、印象があります。

本はいっぱいでてるけど・・・
362 名前：教えて君 [03/02/03 18:18]: A*x=0 ここで　Aは、11*11　の行列　でこの線形方程式を解くには、どの関数を
使ったらいいですか
363 名前：１３２人目の素数さん [03/02/04 01:05]: 外国の数学本を読む人へ　５冊（入門書から一般書まで）　

page3.auctions.yahoo.co.jp/jp/auction/c30452348
364 名前：１３２人目の素数さん [03/02/04 14:44]: 関数の定義で、1つの変数をあたえて、２つの戻り値を得ようと思います。
関数の内容は、自分で定義した関数やSolveなど２１個の評価を行い、最終的には
2つ値を求めるものです。1つの式にできれば、問題ないのですが、それはできません
お願いします
365 名前：１３２人目の素数さん mailto:sage [03/02/04 22:30]: >364
単にリストにすればいいような気がするが。。。
366 名前：１３２人目の素数さん [03/02/04 23:49]: 364のものです。
foo[n_] := Module[{}, ]を使って解決しました。
367 名前：１３２人目の素数さん [03/02/05 17:16]: A.x=0 Aは係数行列　　ｘは変数ベクトルをあらわします。
この方程式をとくのに、最適な関数を、教えてください
368 名前：通りがかり [03/02/05 19:14]: >367
NullSpace[A]:
Ax=0を満たす自明でない解となるベクトルのリスト
369 名前：通りがかり [03/02/06 00:03]: >362
>367

>368 の方法で良いけど、ここも参考になるかも。
www53.tok2.com/home/xxxxxxxx/math/026.html
370 名前：１３２人目の素数さん [03/02/06 01:21]: >362（=>367だとは思うが）

det(A)≠0なら
x=(0,0,0,0,0,0,0,0,0,0,0)'

det(A)=0なら
NullSpace[A]で基本解の導出
371 名前：１３２人目の素数さん [03/02/07 18:13]: mathematica でPlotしたグラフをWordで表示する方法を教えてもらえませんか
372 名前：１３２人目の素数さん [03/02/07 19:05]: プロットしたグラフを右クリックして
copy as より bitmap で保存して word に貼り付ければ？
373 名前：１３２人目の素数さん [03/02/07 20:20]: 372さん、ありがとうございます
374 名前：１３２人目の素数さん [03/02/08 01:15]: ボイコット・ザ・Mathematica。
もっとフリーソフトをつかう、作るべし。
MathematicaにせよMapleにせよ囲い込まれるだけだし、
ある日突然サポートがとまったりするリスクを考えよ。
バージョンアップで毎年金を毟り取られるのに快感を覚えるな。
他と著しく違う変なシンタックスに染まっても、システムが
替われば、そのスキルは無用の長物。
　もっと全数式処理ユーザーが結束して、どうでもよいところは
共通仕様化することを要求していくべきだし、なるべくフリーな
ソースのある、内部で採用で採用されているアルゴリズムが
公開されているものをつかうようにしないと、本当に中世の
暗黒時代が再来するだけだよ。
375 名前：１３２人目の素数さん [03/02/08 09:54]: どなたか
Mathematicaの入門書を教えてください
どうしても自力じゃ使いこなせません
よろしくお願いします
あるていど使いこなしてみたいんです
376 名前：１３２人目の素数さん [03/02/08 13:57]: はやわかり　Mathematica　第２版　榊原　進
は、どうですか
377 名前：１３２人目の素数さん [03/02/08 14:09]: メイプル　マスキャド

と比較してどうなんでしょ？
378 名前：１３２人目の素数さん [03/02/08 14:16]: ますます出来るMathmatica　＝標準版＝って
どうですか
わかりやすいですか？
でも14000は高いですよね？
379 名前：１３２人目の素数さん mailto:sage [03/02/08 22:59]: >374 じゃ、 MAXIMA ね。
まあ、正直僕も Mathematica の税金的な面は痛い所。
昔は３万円！でバージョンアップしてくれたのにね。。。
特に日本だと代理店経由しないとだめで、個人の気楽な
用途からどんどん遠ざかります。

>375
何したいのかによりますよ。
あ）とにかく操作入門レベル。コンピューター使用暦ゼロに近い
い）直面した問題に対処したい。コンピューター使用暦数年。
う）情報処理系的側面重視。コンピュータープログラミング要望。
でいろいろありますね。

まあ自分が知っている分野で、 Mathematica という言葉が
表題にあれば、大抵簡単な入門はのっているよ
380 名前：１３２人目の素数さん [03/02/08 23:32]: 理工系の教養としての情報科学―Mathematicaによる入門
ローマンメーダー (著), Roman Maeder (原著), 時田節 (翻訳)
www.amazon.co.jp/exec/obidos/ASIN/4894717751/qid=1044714600/sr=1-1/ref=sr_1_0_1/250-2915690-8097054
381 名前：１３２人目の素数さん [03/02/08 23:33]: >>380
Mathematicaでのプログラミングにはお勧め
382 名前：１３２人目の素数さん [03/02/09 11:42]: Mathematicaって
普通のアプリにあるようなバグフィックスのアップデータって無いの？
バグなんか無いって事なんだろうか。
383 名前：１３２人目の素数さん [03/02/09 12:20]: >>382

バージョンアップ時に小出しにバグフィックス致します。
バージョンアップ料金は１０万円以上かかります。
384 名前：１３２人目の素数さん [03/02/09 12:40]: 俺の持ってるのは3.0.1.1xってなってる
3.0.0.0持ってる人もいるの？
385 名前：１３２人目の素数さん mailto:sage [03/02/09 15:02]: MuPADってだめなの？
386 名前：１３２人目の素数さん [03/02/09 19:36]: まてぃまてぃかで
世界征服したいんですけど？
387 名前：１３２人目の素数さん [03/02/11 18:22]: Plot した図を　Showで表示します。
２つの違う関数を表示するのだけど、その線がどっちかって説明とA4に印刷するた
めの方法を教えてもらえませんか。
簡単に言えば、レポートに提出するためのものです
388 名前：１３２人目の素数さん mailto:sage [03/02/11 22:38]: >>385
フリーなのは学生と教員だけだろ。
389 名前：１３２人目の素数さん mailto:sage [03/02/11 23:24]: >387
help に線の太さやパターンを変える方法がのってるよ。
wolfram のサイトにいけば、日本語で読める可能性有り。
390 名前：１３２人目の素数さん mailto:sage [03/02/15 23:16]: >>388
個人利用はフリーと聞いたが
391 名前：１３２人目の素数さん [03/02/21 15:50]: risa.cs.ehime-u.ac.jp/index.html
392 名前：１３２人目の素数さん [03/03/01 00:06]: mathematicaのstudent版って卒業してしまうと使えなくなってしまうのですか？
買ってから４年間、みたいな制約があるんでしょうか？
393 名前：Lisper mailto:sage [03/03/08 15:57]: >>392 GNU Maximaなら永遠に無料でつ．シミュレーション板Maxima擦れを参照．
どうせMathematicaをそれ程クリティカルな業務に使ってないのでしょうから，Maximaで十分でつよ．
3Dグラフがマウスでぐりぐり回せますぜ．どうでつか．
394 名前：山崎渉 mailto:（＾＾） [03/03/13 13:23]: （＾＾）
395 名前：１３２人目の素数さん mailto:sage [03/03/16 08:28]
396 名前：とある馬鹿　 ◆BAKAB.w.so mailto:ヽ( ･∀･)ﾉ●ｳﾝｺー [03/03/16 12:33]: Mathematica　を早く使いこなせるように私はなりたい　ううう　( つДT)

大学でしかできないからなぁ・・・

やっぱり自腹で買ったほうがいいのかな
それとも研究室の金で買ったほうがいいのかな･･･
397 名前：１３２人目の素数さん mailto:sage [03/03/16 13:00]: 大学でできれば十分ジャン。

たまには家に帰りたいもんだ。
398 名前：とある馬鹿　 ◆BAKAB.w.so mailto:sage [03/03/16 14:16]: >たまには家に帰りたいもんだ

なんか（・∀・）ｶｺｲｲ!!!!!
399 名前：質問！ [03/03/16 14:27]: （関数を定義する）
仮に・・・a_{i}:実数
s[a_{1},a_{2},....,a_{n}]:＝a_{n-1}　：　与えられた数のn-1番目を取り出す
という関数sを定義したいですが，変数をリストとしてではなく，
与えたいと思ってます．どうすれば良いのでしょうか？
ちなみに，一般のｎについてできるようにしたいです．
ちょうど，Max[a_{1},a_{2},....,a_{n}]の入力の仕方で
答えを得るようにしたいということです．
s[{a_{1},a_{2},....,a_{n}}]から，a_{n-1}を得るのは，簡単ですけど，
s[a_{1},a_{2},....,a_{n}]：これのときは？＃ｎ：一般の自然数＃
400 名前：399 [03/03/16 18:21]: 誰かわかる方いませんか？
401 名前：１３２人目の素数さん [03/03/16 20:27]: >>399-400
これでいいだろ。

In[10]:=
s[x__]:=List[x][[-2]]
In[11]:=
s[1,2,3,4,5,6]
Out[11]=
5
402 名前：１３２人目の素数さん [03/03/16 20:40]: これでもいい。

In[16]:=
s[x__]:={x}[[-2]]
In[15]:=
s[1,2,3,4,5,6]
Out[15]=
{1,2,3,4,5,6}
403 名前：402訂正 mailto:sage [03/03/16 20:53]: In[16]:=
s[x__]:={x}[[-2]]

In[17]:=
s[1,2,3,4,5,6]

Out[17]=
5
404 名前：399-400 [03/03/16 21:43]: なるほど！！！
そういうことで回避できるんだー！
401-403 : ありがとうございました!
うーんこれは使える・・・
405 名前：１３２人目の素数さん [03/03/18 22:30]: >>183
先生, 質問です.
StudentからProfessionalにUpgradeすると,
パスワードの発行形態は後者に準じたものに変わるのでしょうか?
406 名前：１３２人目の素数さん mailto:sage [03/03/18 22:39]: >>393
使えるのなら、使っていいんじゃないの？
version3.0 for Students を持っているけど、いまだに使っているよ。
インストールしなおす時には、もちろん、きーじぇ・・・スミマセン。
407 名前：kawabata君の友人 [03/03/19 00:48]: >>405
はい、発行形態がかわります。
ただし、Proffessionalへアップグレードしても
　(１)ホームユース(名前は確かじゃないですが「一本」のライセンスで、
　　いつも使用している研究室や職場だけでなく、自宅でも使用することが
　　できるシステム)(この場合は自宅の端末専用に1年間有効のMathIDベースの
　　パスワードが発行されます。当然ＣＤもMathID仕様のものが追加で交付
されます。１年ごとに追加料金がかかります
　(２)そのほかにも北米ライセンスを、一時的に他の地域で使用する場合、
　　　または、その逆で、インターナショナルライセンスを北米で使う場合
　　　などにも、時限パスワードが追加交付されるはずです
408 名前：１３２人目の素数さん [03/03/19 01:24]: マス　書いてもいいか
409 名前：１３２人目の素数さん mailto:sage [03/03/19 09:47]: >>407
ということは,
For StudentからProffessionalへアップグレードしても"MathIDは常に必要"ということですか?
最初からProffessionalの場合みたく, 実質的には(自分の所有する)複数のPCに入れ放題ということにはならないのか……
410 名前：kawabata君の友人 [03/03/20 00:43]: >>409
プロフェッショナルはMathID依存でない１８桁～２５桁のパスワードがでる
はずです。しかしホームユースや地域外での利用に対しては、別途料金がかかり
その料金徴収、有効期限管理、追加で使用を認めるマシンを限定し特定する必要
があるため、MathID形態のＣＤ並びにパスワードを追加交付しているのでしょう。
411 名前：baggio [03/03/20 02:06]: 質問です。mathematicaを使って、非線形連立方程式を解きたいのですが、どうしたら
いいのでしょうか？初心者の質問で申し訳ないですが、教えてください。
412 名前：１３２人目の素数さん mailto:sage [03/03/20 22:41]: >410
うーん。通常の Professional ではなく、プレミアサポート
（定価３０％、PC版なら約１０万）でホームユースが
出来る、と考えた方がいいでしょう。
413 名前：１３２人目の素数さん [03/03/20 23:44]: >411

代数方程式なら
Solve[方程式のリスト，未知数のリスト]
とします。
それ以外の場合はFindRootを使います。
その場合は解の初期近似を与える必要があります。
FindRoot[方程式のリスト，{未知数,その初期値}のリスト]
414 名前：１３２人目の素数さん [03/03/23 03:14]: Mathematicaで整数計画問題解けますか？
415 名前：baggio [03/03/28 01:10]: >413 132人目の素数さんありがとうございます。
一回やってみます。わざわざ解答ありがとうでした。
416 名前：１３２人目の素数さん [03/04/15 01:53]: Linux版のMathematicaで、スタンドアローンのプロフェッショナル版を
買おうかと思っているのですが、Ｌｉｎｕｘのバージョンをアップデート
したり、あるいはハードを変更したり、マザボードやＣＰＵをとりかえたり
ネットワークカードを変えたり、ＩＰアドレスを変えたりしたら、
使えなくなってしまうのではないかと恐れています。
　何が、保存されていれば（たとえばＭＡＣアドレスとかＩＰアドレス）、
環境を変化させても、ライセンスを取り直さなくても使えるのでしょうか？
417 名前：１３２人目の素数さん [03/04/15 16:13]: MathematicaではLogListPlotやLogLogListPlotのような片対数グラフや両対数グラフを
表示できて便利なんですけど、エラーバー付きのLogListPlotやLogLogListPlotのような
ものはないんでしょうか？

対数表示でない通常のエラーバー付きグラフならErrorListPlotがあるんですが
ErrorLogListPlotやErrorLogLogListPlotのようなものがないのか知りたいです。
418 名前：１３２人目の素数さん mailto:sage [03/04/15 16:25]: >>417
臭うぞ！ Q.manの臭いがプンプンするッ！
419 名前：mathmania ◆uvIGneQQBs [03/04/15 16:45]: >> 417
Exp,Log,ErrorListPlotを組み合わせて自分で関数を作れ。
>> 418
417がQ.manならば、Q.manは変数変換の公式なぞろくに使えないだろう。
420 名前：１３２人目の素数さん [03/04/15 18:01]: Q.manってだれ
421 名前：Q.man [03/04/15 18:23]: ゎたすがQ.manです。
422 名前：１３２人目の素数さん [03/04/16 10:14]: あっそ
423 名前：P.man [03/04/16 23:09]: ゎたすがP.manです。
424 名前：１３２人目の素数さん [03/04/19 18:31]: (･∀･)ｹﾞﾊﾊﾊﾊ
425 名前：山崎渉 mailto:（＾＾）sage [03/04/20 03:59]: 　　 ∧＿∧
　　（　　＾＾）＜ぬるぽ（＾＾）
426 名前：apico [03/04/23 20:52]: Mathematicaで３次元のグラフを表示させる方法を教えて下さい。
以前販売元へ質問状をFAXとメール送ったのですが、その後連絡がありません。
日本人には冷たいのかな？？
閑話休題、
[1]曲面z=Sqrt[x*y]と３平面z=0,x=2,y=3によって囲まれた立体を表示させたい
のと
[2]曲面z=Sqrt[a^2-y^2](y≧a/2)と４平面x=0,x=2a,y=a/2,z=0によって囲まれる
立体を表示させたいのです。
表示させようと思ってすでに半年。自分でもどうやって聞いたらよいのか分からな
いまま経ってしまいました。もしこれが無理なら日本語でこういった質問の出来る
メールアドレスを教えて下さい。御願いします。
427 名前：mathmania ◆uvIGneQQBs [03/04/25 18:11]: >> 426
表示しようと思って半年もかかるくらいなら、手書きの方が正確で早い。
だが、それでもコンピュータで表示したいなら、
? Plot3D
を使うか、
あるいは
<< Graphics`をやって
? ParametricPlot3D
でヘルプを出して関数を理解してくれ。
428 名前：１３２人目の素数さん mailto:sage [03/04/25 18:50]: >>426
誰でも、すぐできることだから。
429 名前：１３２人目の素数さん mailto:sage [03/04/25 19:45]: はじめまして
今日mathematica 4.2 for students買ってきました
今後お世話になるかもしれませんがよろぴく
430 名前：１３２人目の素数さん mailto:sage [03/04/27 20:40]: >>427, 428
「囲まれた立体」を表示させるのは
Plot3Dでは不可能だし
ParametricPlot3Dでも[1], [2]を表示するのはすくなくとも困難だ。
まぁ428には不可能だな。
431 名前：１３２人目の素数さん mailto:sage [03/04/28 22:16]: >>426
? ClipFill
432 名前：１３２人目の素数さん [03/04/29 07:49]: ClipFillねー
これでいいのかな？

Plot3D[Sqrt[x*y],{x,0,2},{y,0,3},
　PlotRange->{0,3},
　PlotPoints->30,
　ClipFill->Automatic,
　BoxRatios->{2,3,3}]
433 名前：１３２人目の素数さん mailto:sage [03/04/29 22:25]: >>432
Clipしてないのだが？
434 名前：１３２人目の素数さん [03/05/01 18:17]: mathematicaって将来性あるの？
435 名前：１３２人目の素数さん [03/05/02 08:41]: MaximaやMuPADなど無料のソフトがもっと
グレードアップしてくると
今の商売は難しくなるね。
我々はその方が有り難いけど。
436 名前：１３２人目の素数さん [03/05/02 08:42]: 優良スレに認定されてます

◇ 「知」の欺瞞
academy2.2ch.net/test/read.cgi/philo/1047993277/
◇ アーベル賞にセール
news2.2ch.net/test/read.cgi/newsplus/1049380704/
◇ 数学の本
science.2ch.net/test/read.cgi/math/1044371030/
◇ 代数学
science.2ch.net/test/read.cgi/math/1045779496/
◇ 数学セミナー
science.2ch.net/test/read.cgi/math/1050563258/
◇ Poincare 予想解決
science.2ch.net/test/read.cgi/math/1049433413/
◇ Lie群・Lie環
science.2ch.net/test/read.cgi/math/997388738/
◇ グロタンディック
science.2ch.net/test/read.cgi/math/1011949239/
◇ 基礎論
science.2ch.net/test/read.cgi/math/1043049229/
◇ アーベル賞にセール
science.2ch.net/test/read.cgi/math/1049442770/
◇ 素数判定は「決定的」多項式時間
science.2ch.net/test/read.cgi/math/1028813059/
◇ スペクトル系列
science.2ch.net/test/read.cgi/math/1031726106/
◇ unix と数学
pc.2ch.net/test/read.cgi/unix/1043629084/
437 名前：１３２人目の素数さん [03/05/03 20:39]: いまいち盛り上がらないので。
Mathematica programmingコンテスト。
お題は、

トランプ(Range[52])をシャッフルする関数を書け。

とか言ってみる。
438 名前：１３２人目の素数さん [03/05/03 20:43]: ＞４３７
ドナルド・トランプですか？
439 名前：437 mailto:sage [03/05/03 20:53]: >>438
普通のトランプです。
ドナルドさんが52人おられるのでしたらそれでも結構です。
440 名前：１３２人目の素数さん mailto:sage [03/05/15 01:19]: >>437
こんなのとか。

Shuffle[x_List] :=
Prepend[Shuffle[Rest[#]], First[#]] &@
RotateLeft[x, Random[Integer, Length[x] - 1]]
Shuffle[{}] = {}

Shuffle[Range[52]]
441 名前：山崎渉 mailto:（＾＾） [03/05/21 22:59]: ━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━―
442 名前：１３２人目の素数さん [03/05/25 05:59]: 29
443 名前：山崎渉 mailto:（＾＾） [03/05/28 14:44]: 　　　　∧＿∧
ﾋﾟｭ.ｰ　(　　＾＾）＜これからも僕を応援して下さいね（＾＾）。
　　＝〔~∪￣￣〕
　　＝ ◎――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉
444 名前：１３２人目の素数さん [03/06/02 16:04]: Mathematicaが原因でMATHEMATICA.EXEにエラーが発生しました。
Mahtmeticaは終了します。
問題が解決しない場合は、コンピュータを再起動してください。

と出たので再起動したのですが、何回再起動しても同じエラーが出て
起動してくれません。

フォント関係の設定をいじった後このエラーが出たので
元に戻したいと思っているのですが
設定ファイルはどこにあるかわかりませんか？
レジストリは探してみたのですが、わかりませんでした。

Windows MeでMathemtica4.2のシングルライセンスのを使っています。
445 名前：１３２人目の素数さん mailto:??? [03/06/02 18:40]: 再インストールしる!
446 名前：１３２人目の素数さん [03/06/02 20:41]: Mathematicaって使いにくくない？ライブラリとか凄いとは思うけど
447 名前：mathmania ◆uvIGneQQBs [03/06/02 20:46]: コマンドが長いのは確かに使いづらい。
もっとも、私は最近Mathematicaをやっていない。
C言語(ForTranも少し)を使ってる。
理由はというと、、、Mathematicaが使えなくなったから。
まぁ、使いやすさの面でいえば、大きなプログラムには、C言語の方がいいだろう。
448 名前：444 mailto:sage [03/06/02 23:50]: >>445
再インストールしても同じエラーが出るんです。
設定ファイルが残ってるんじゃないかと疑ってるのですが
その場所が分かりません。
449 名前：１３２人目の素数さん mailto:sage [03/06/03 22:42]: >>446
原始的なプログラミング言語と不統一な仕様の
mapleよりは使いやすいと思う。
言語としてつかえてなんぼだから
素人にはお薦めできない、って感じはあるかもね。

>>447
補完機能をつかうべし。
450 名前：１３２人目の素数さん [03/06/07 17:45]: 初心者にお勧めのMathemticaの
本は何でしょうか？
451 名前：１３２人目の素数さん mailto:??? [03/06/07 21:09]: >448
OSを再インストールしる!
それでだいじょうぶだ! 安心しろ

>450
白石さんのがいいぞ!
安心して買え
452 名前：１３２人目の素数さん [03/06/07 21:25]: あなたが探してるのってこれだよね？二日間無料だしね♪

alink3.uic.to/user/angeler.html
453 名前：１３２人目の素数さん [03/06/07 22:15]: mathmaniaスレかと思った

　　　　　　　　　　　　と言ってみる
454 名前：１３２人目の素数さん mailto:sage [03/06/08 00:02]: >>450
Mathematicaプログラミング J.ゲイロード･N.カーミン･R.ウエリン
を読んで、その後で
プログラミングMATHEMATICA. Ｒ・メーダー著、時田節
を読むとよくわかると思われ。

特にメーダーの本は(全くの初心者以外には)一読の価値あり。
非常に下手な翻訳だけど。
455 名前：１３２人目の素数さん [03/06/08 00:09]: >>248
Windows Me はやめた方がいいですよ。
2KかXPにしましょう。
（マジレス）
456 名前：１３２人目の素数さん [03/06/08 02:09]: Mathematicaって有効数字が６桁の表示で、例えば
1.23456789とか1234.56789って直接入力すると
1.23456や1234.56と表示されます。

＃　表示だけがそうなるだけで内部的には入力した数のままです

これを以降の入力からすべて指定した有効数字で表示させたいのですが
どうやったら良いのでしょうか？NumberFormではそれで指定した場合だけ
表示されるのでこれを以降の入力すべてで有効にしたいのでよろしく
お願いいたします。
457 名前：１３２人目の素数さん [03/06/08 10:38]: >>451 >>454

どうもありがとう！
458 名前：１３２人目の素数さん [03/06/08 10:49]: >>456
ちょっと邪魔くさいけど、関数を定義すればどう？たとえば
　F[x_]:=N[x, 30];
　F[E^3]
とするとe^3が30桁表示されるはず。
459 名前：mathmania ◆uvIGneQQBs [03/06/08 13:02]: 456の方法でも有効数字は確かに表示される。
だが、SetPrecision[E^3,Infinity]か、SetAccuracy[E^3,Infinity]はどうだろう？
まぁ、分数表示になるが、その数に2^52(?)を掛ければ小数点表示がわかるだろう。

あぁ、それと、私はMathematica復活したわけではないので。
460 名前：444 mailto:sage [03/06/08 21:28]: >>451
結局OS再インストールしましたよ。。。
したくなかったけど。

>>455
大学のPCだから自由にできません。
経費で2000買ってくれなかった。
461 名前：１３２人目の素数さん [03/06/09 17:46]: >>458-459
Mathematicaでの以降の入力すべてで指定した有効数字で表示させたいということ
でしょうから、その方法ではその関数内の引数にだけ有効であって意味がありません。
何かシステムの設定でできそうな気もしますがそれらしいものが見当たりませんね。
どうやるんだろうか？
462 名前：１３２人目の素数さん [03/06/09 18:00]: >>461
SetPrecisionは計算精度を変更させるだけであってその表示は依然として
有効数字６桁です。
463 名前：１３２人目の素数さん [03/06/10 19:55]: >>456
編集→環境設定→書式設定→式の書式→表示設定→PrintPrecision
464 名前：１３２人目の素数さん [03/06/10 23:57]: >>463
書式設定でしたか。見落としてました。大変助かりました。ありがとうございます。
465 名前：mathmania [03/06/17 23:17]: 駄スレ漏れスレage
466 名前：１３２人目の素数さん [03/06/18 00:04]: あなたが探してるのってこれだよね？この中にあったよ♪
endou.kir.jp/betu/linkvp/linkvp.html
467 名前：１３２人目の素数さん [03/06/21 17:16]: 4次元データを視覚化したいのですが、mathematicaに限らず
何かよい方法(ツール)はありませんか?

等高線図を複数枚、平行になるよう並べて
ナナメ上から見下ろすようにすればよいと思うのですが、
よいアプリケーションを探しています。
他の方法でも、可能ならば何でもいいです。
468 名前：１３２人目の素数さん mailto:sage [03/06/22 17:01]: >>467
いってることが良くわかりません。等高線図を斜めから見下ろす？

mathematicaでの実現ならば、三次元にスライスしたときの
データをどうやって表示するかという話になるのだろうけど、
それはデータの形式に依存するからなんともいえないし。
469 名前：467 [03/06/23 13:26]: >>468
間違いがありました。すみません。
まず画面に出したいのは4次元データ中の3次元分なので、
「3次元データを表示したい」と訂正します。
データは格子状データです。
古い人間なのでfortran形式で言うと
do k=1,nz
do j=1,ny
do i=1,nx
write(6,*) x(i),y(j),z(k),data(i,j,k)
end do
end do
end do
の形式の等間隔のデータと思って下さい。
やりたいのは、これまた解かりにくい表現かもしれませんが
gnuplotで
set contour base
unset surface
set view 60,10
くらいの時にプロットした等高線図を複数枚、一つの画像にプロットしたいのです。
470 名前：１３２人目の素数さん mailto:sage [03/06/23 20:03]: >>469
gnuplot云々はわかりませんが、
ListContourPlot3Dかな?
471 名前：１３２人目の素数さん mailto:sage [03/06/24 22:42]: Ver.5 ですね。
それに $1880 5.0 Eng vs \448800 日本語版 4.2only
US ってこれでも Premier Service 付きだって。。。。

潮時かね
472 名前：467 [03/06/26 18:39]: 忙しかったので、見るの遅くなって済みません。

>>470
>ListContourPlot3Dかな?

手持ちのマニュアルには
ListContourPlotとListPlot3Dしか見付かりません。。
(ver5からの新機能ですか?)
それら二つから推測すると、
ListContourPlot3D[{{{},{}},{{},{}}}]
のような3重ネストの配列を入れてプロットするのでしょうか?

値段高い…
473 名前：１３２人目の素数さん mailto:sage [03/06/26 23:19]: >>472
4.2にはあります。
Onlinehelpで出てくるはずです。

または、
<<Graphics`
?ListContourPlot3D
474 名前：467 [03/06/27 14:28]: みなさま、御親切に
本当にありがとうございます！
ちょっといろいろ試してみます。
475 名前：１３２人目の素数さん [03/06/27 15:22]: 足してnになる1以上の整数の組み合わせを並べる時、並べ方が最大になるような
整数の組み合わせをMathematica（あるいはC言語でも）で求める方法を教えてください。

わかりにくいので例を挙げますと　n=5　のとき
足して5になる1以上の整数の組み合わせは
5, 41, 32, 311, 221, 2111, 11111　であり
並べ方が何通りあるかはそれぞれ　1, 2, 2, 3, 3, 4, 1
で最大値は　2111　の場合で　4　に、
n=6の場合は　321　と　2211　の場合で最大値は　6　になります。
一般に、数字の個数（311なら3個、11111なら5個）がm個、1がk_1個、2がk_2個、…nがk_n個の
場合、並べ方はm!/(k_1!・k_2!・・・k_n!)になります。
476 名前：mathmania ◆uvIGneQQBs [03/06/27 15:36]: Re:>475
最大値をどうやって計算するかは分かるだろう。
あとは、k_1,…,k_nの与え方だが、
これはk_n,k_{n-1},…,k_2,k_1の順に再帰的に与えればいいのではないか？
n=5のときを例に挙げると、
まず、k_5=1ととり、1*k_1+2*k_2+3*k_3+4*k_4=0になるk_1,k_2,k_3,k_4の組合わせは0,0,0,0しかない。
次にk_5=0ととるとき、1*k_1+2*k_2+3*k_3+4*k_4=5となるk_1,k_2,k_3,k_4を選ぶことを考える。(a)
(a)の中でk_4=1ととり、1*k_1+2*k_2+3*k_3=1となるようにk_1,k_2,k_3を選ぶようにする。(k_1=0,k_2=k_3=0しかない。)
(a)の中で、今度はk_4=0ととり、1*k_1+2*k_2+3*k_3=5となるようにk_1,k_2,k_3を選ぶ。(b)
(b)の中で、k_3=1ととると、(k_2,k_1)=(1,0)と(k_2,k_1)=(0,2)となる。
(b)の中で、k_3=0ととると、k_2=2,1,0の3通りにわかれる。
k_2=2ととるとき、k_1=1となり、k_2=1ととるとき、k_1=3となり、k_2=0ととると、k_1=5となる。

この文では分かりにくいかもしれないが、とにかくやってみてくれ。
477 名前：474 mailto:sage [03/06/27 16:10]: >>476
その考え方はわかりました。
初心者のものでプログラムの書き方を勉強してきます。
478 名前：１３２人目の素数さん mailto:sage [03/06/27 16:11]: 477=475です
479 名前：１３２人目の素数さん mailto:sage [03/06/27 17:06]: >>475
どうぞ。

<< DiscreteMath`Combinatorica`
<< Statistics`DataManipulation`
Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@
Partitions[5]
480 名前：１３２人目の素数さん mailto:??? [03/06/27 20:12]: Mathematica Ver.5の情報はどこにあるんですか?

何が変わるのでしようか
おしえてね!
481 名前：１３２人目の素数さん [03/06/27 20:24]: 過去にレスがついていなかったのですが、
ますます出来るMATHEMATICA-標準偏-
www.catkk.co.jp/mathcd/index.html
どうでしょうか？
482 名前：１３２人目の素数さん mailto:sage [03/06/27 20:25]: >>480
www.wolfram.com
483 名前：475 mailto:sage [03/06/27 23:37]: >>479
よくわからないのですが、Partitions[5]の　5　を　n　に変えれば
足してnになる1以上の整数の組み合わせの並べ方が最大に
なる整数の組み合わせが求められるんでしょうか？
484 名前：483 mailto:sage [03/06/27 23:39]: すいません。”　5　を　ｎ　に”　は　”　5　を任意の整数に”という意味です
485 名前：１３２人目の素数さん mailto:sage [03/06/27 23:43]: >>483
f[n_Integer] := Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@
Partitions[n]
ってこと?

もしかして、数学っぽいnについての式を欲しがってるんですか?
それならそれは数学の問題。
プログラムで求めることはできません。
Cでなんて論外。
486 名前：479 == 485 mailto:sage [03/06/27 23:45]: >>484
See f in >>485.
487 名前：475 mailto:sage [03/06/27 23:54]: >>486
ありがとうございます。手元にMathematicaがないので
少し勉強して、学校ので試してみたいと思います。
数学的に（nを使って）解くことが出来たらそれが一番いいのですが、
それはここで聞くことではないので他のスレで聞いてみます。
488 名前：１３２人目の素数さん [03/07/04 13:36]: 再セットアップしたらアクティベーションできなくなった
489 名前：１３２人目の素数さん mailto:??? [03/07/04 19:39]: www.wolfram.com/

うわーん!
Ver.4.2の情報しかないじゃないか。
しかも日本語だー
490 名前：１３２人目の素数さん mailto:sage [03/07/05 01:36]: >489
どこかに English ってあるからそこから英語ページへ。
トップにあるよ。
まあ大規模行列と数値計算の速度改良が主ですね。
試行錯誤でアルゴリズム作るときは、全作業時間において
コンパイル系とためはるつもりなのかな

なお日本では日本語版のみしか売れない。
誰だよな日本語ヘルプキット強制バンドル計画したやつは。
491 名前：１３２人目の素数さん mailto:??? [03/07/05 21:11]: ttp://www.wolfram.com/products/mathematica/newin5/
に五版の概要がありますね。

英語版でいいから新しいやつを動かしてみたいねえ。
ほんと。
492 名前：１３２人目の素数さん [03/07/05 23:02]: 私が作りました♪
nuts.free-city.net/index.html
493 名前：488 mailto:sage [03/07/07 15:52]: 誰か助けてくれ
494 名前：１３２人目の素数さん mailto:sage [03/07/07 19:47]: >>493
サポートに連絡しろ。
495 名前：１３２人目の素数さん mailto:sage [03/07/08 01:13]: >>494
わかったやってみる
㌧くす
496 名前：475 mailto:sage [03/07/08 16:22]: >>479
Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@
Partitions[5]　　を実行したところ

Part::"partd": "パート指定5[[1]]の長さはオブジェクトの深さを超えています．"

　 1/5[[1]]!　

となったのですが、直すべきところがわからないので教えてください。
497 名前：１３２人目の素数さん mailto:sage [03/07/08 21:15]: >>496
ちゃんと上の二行を入力してください。

# そのエラーはPartitions[5]が評価されていないため。
498 名前：１３２人目の素数さん [03/07/08 22:23]: 最大尤度法で推定をしていて、結構複雑なモデルなんですが、推定したい変数が５個
あるんです。それで、尤度関数を作って、偏微分してFindRootつかって求めようした
んですが、出来ないんです。

そこで、正規分布使って、普通に正規分布を推定してみようって思って、同じような
方法を使ったんですが、これでもうまくいきません。エラーとしては

{{初期値の尤度}},{{初期値の尤度}}は、｛a,b}={1,2}では長さ２のリストになって
いません

というような感じです。最大尤度法をやってるプログラムのサンプルどっかないかなあ
って探したんですが、見つからないし・・・・。泣きそうなんですけど、どなたかお助け
ください。
499 名前：１３２人目の素数さん mailto:sage [03/07/08 22:43]: >>498
その説明で他人がなにかわかると思うか？

> んですが、出来ないんです。

ふーん出来ないんだ。(どこがどう出来ないのか全くわからん)

> {{初期値の尤度}},{{初期値の尤度}}は、｛a,b}={1,2}では長さ２のリストになって

文法間違いだな。 (syntax errorなんだからソース見ないとどうしようもない)

> って探したんですが、見つからないし・・・・。泣きそうなんですけど、どなたかお助け

まぁ冷静になれ。
500 名前：１３２人目の素数さん mailto:sage [03/07/08 23:08]: >>497
わかりました。
501 名前：498 [03/07/08 23:11]: そうですか・・・・すいません。どっかサンプルプログラム探すしかない
ですかね・・・
正規分布を推定する場合、乱数を出して、その後、尤度関数を出せたとして
K1[a_,b_]=D[L[a,b],a]
K2[a_,b_]=D[L[a,b],b]
として、
FindRoot[{K1=0,K2=0},{a,-1},{b,2}]
としてるんですけど。。この部分では、間違いがないとは思ってるんですが
尤度関数が間違ってるのかな・・・・。はあー泣きそう。。。
502 名前：１３２人目の素数さん [03/07/08 23:12]: ☆頑張ってまーす！！☆女の子が作ったサイトです☆
　　　　　　　☆見て見て！！
yahooo.s2.x-beat.com/linkvp/linkvp.html
503 名前：１３２人目の素数さん mailto:sage [03/07/08 23:20]: >>501
Mathematicaの基礎からやりなおしってかんじ。
とりあえずダメなところは指摘しとくが、
マニュアルか教科書(上でも挙がってる)読み直せ。

> K1[a_,b_]=D[L[a,b],a]
= と := の違いが良くわかってるのなら何もいわないが。

> FindRoot[{K1=0,K2=0},{a,-1},{b,2}]
K1 = 0 は K1 という変数を 0にする。
書くなら {K1[a, b] == 0, K2[a,b] == 0} 。
504 名前：498 [03/07/08 23:20]: あー勿論
FindRoot[{K1[a,b]=0,K2[a,b]=0},{a,-1},{b,2}]
としてます。。。。。

サンプルプログラム載ってるサイトってないですかねえ・・・・。。。
505 名前：503 mailto:sage [03/07/08 23:22]: >>504
その2行目をもういちど読み直せ。
506 名前：498 [03/07/08 23:23]: :=と=の違いってえーと、すみません、良く分ってないです。。。
確か:=は、なんか固定するみたいな意味だった気が・・・・
507 名前：498 [03/07/08 23:25]: >>505
==ですね・・・すみません、それは、やってます・・・・。
508 名前：503 mailto:sage [03/07/08 23:28]: >>506
で、聞いてる方としてはあいかわらずどこがうまく行かないのか
わからないんだが。

で、どううまくいかいの？
エラーメッセージは？
509 名前：498 [03/07/08 23:37]: エラーメッセージは、FindRootのところなんですよー。
>>498に書いたようなエラーメッセージが出ます。
FindRootの文法は、合ってると思うんですが。。。
510 名前：503 mailto:sage [03/07/08 23:45]: >>509
498を読んだところでは第一パラグラフと第二パラグラフは
別の手段を試してるようにしか読めない(FindRootで失敗したから別の手段)。
どの関数でエラーが出力されたかもわからない。
エラーがあなたの言語に翻訳されているがそれが適切かどうかわからない。
すでに書いたがsyntax errorはソースをみないと何もわからない。

技術的なことを質問する方法をしらない、
シンタックスエラーが自分で解決できない、
レベルだったらソース張るべきだよ。

変にプライドがあるのかそういう「初心者」は多いけどね。
情報を出さずに、自分のやってることは間違ってない!、みたいな。

たんに誰かが作ったライブラリが欲しいのなら、
まず、googleなりwolframのサイトなりを検索し、
もっと適切な質問のしかたをすべき。

おれは今日はもう寝ます。
511 名前：503 mailto:sage [03/07/09 00:02]: >>509
寝ようと思ってたが気づいたので。
多分 L がリストを返してるからだ。
512 名前：１３２人目の素数さん [03/07/09 14:24]: ndist2=NormalDistribution[0,4];
p1=Table[ Random[ndist2],{100}];

L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}]

K1[a_,b_]=D[L[a,b],a]
K2[a_,b_]=D[L[a,b],b]

FindRoot[{R1[a,b]==0,R2[a,b]==0},{a,{1,2}},{b,{-2,4}}]

ソースは、これです。
>>511
あーそれは、どうすれば、解消できますかね。リストで返さずに関数として認識
してもらわないとつらいっすね・・・。
513 名前：498 [03/07/09 15:12]: あっ最後の式は、
FindRoot[{K1[a,b]==0,K2[a,b]==0},{a,{1,2}},{b,{-2,4}}]
です。
514 名前：１３２人目の素数さん [03/07/09 16:05]: 今課題をやっています。

誰か偏微分係数の出し方を教えてください。

ちなみにMathematika 4.1を使用しています。
ヘルプの使い方がわからなくて・・・
515 名前：514 [03/07/09 16:19]: 解決すますた。
失礼。
516 名前：503 mailto:sage [03/07/09 17:37]: >>498
ソースはコピペしてそのまま張れ。
タイプミスもあるしな。

πってPiだよな？
<<Statistics`NormalDistribution`
してるんだよな？

原因はp1のまえにある'{'だ。'{'はリスト構築子。
L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}]

↓におきかえろ。
L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-((p1[[i]] - b)^2/(2 *a^2)))],{i,1,100}]
517 名前：１３２人目の素数さん [03/07/09 21:13]: このスレッドの人から見てＳ言語はどうですか？
社会学系だとmatlabやmathematicaより使われていますが…
518 名前：498 [03/07/10 16:44]: >>516
ありがとうございました。それで正規分布の推定は、できました。ただ、本当に
求めたいのは、
L[a1_, a2_, b1_, b2_, b3_] = Sum[y1[[i]]*Log[F[a1 - (b1*
x1[[i]] + b2*x2[[i]] + b3* x3[[i]])]] +
y2[[i]]*
Log[F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])] - F[
a1 - (b1*x1[[i]] + b2*x2[[i]] +
b3*x3[[i]])]] +
y3[[i]]*Log[
1 - F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])]], {i, 1,
100}];
という結構複雑なもので、これも｛を（にしてみたんですが、やはりリストで返してる
みたいです。Sumの扱い方が悪いんでしょうか・・・・
てか、見にくくてすみません。
519 名前：１３２人目の素数さん mailto:sage [03/07/10 19:28]: >>518
a1,a2,b1,b2,b3,F,y1,... のどれかにリストが入ってんだろ…
同じことを何度も書くが、部分だけ張るなよ…

プログラミング自体初めてなのかも知れないが、レベル低すぎる。
一冊Mathematicaの本読め。
520 名前：498 [03/07/10 21:10]: すみません、一遍に貼れないんで２回に分けますけど、これが全部です。。。

<<Statistics`ContinuousDistributions`

ndist=NormalDistribution[0,1];

x1=Table[ Random[ndist],{100}];

x2=Table[ Random[ndist],{100}];

x3=Table[ Random[ndist],{100}];

y=Table[2*x1[[i]]+3*x2[[i]]-x3[[i]],{i,100}];

y1=Table[If[y[[i]]<=-5,1,0],{i,100}];
y2=Table[If[-5<y[[i]]<=5,1,0],{i,100}];

y3=Table[If[5<y[[i]],1,0],{i,100}];
521 名前：498 [03/07/10 21:11]: F[x_] = Integrate[Exp[z]/{1 + Exp[z]}^2, {z, -Infinity, x}];

L[a1_,a2_,b1_,b2_,b3_]=Sum[y1[[i]]*Log[F[a1-(
b1*x1[[i]]+b2*x2[[i]]+b3*
x3[[i]])]]+
y2[[i]]*Log[F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]-
F[a1-(b1*x1[[i]]+b2*x2[[i]]+
b3*x3[[i]])]]+y3[[i]]*Log[1-
F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]],{i,1,100}];

K1[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a1];

K2[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a2];

K3[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b1];

K4[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b2];

K5[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b3];

FindRoot[
{K1[a1,a2,b1,b2,b3]\[Equal]0,
K2[a1,a2,b1,b2,b3]\[Equal]0,
K3[a1,a2,b1,b2,b3]\[Equal]0,
K4[a1,a2,b1,b2,b3]\[Equal]0,
K5[a1,a2,b1,b2,b3]\[Equal]0},
{a1,0},{a2,1},{b1,1},{b2,3},{b3,1}
]
522 名前：503 mailto:sage [03/07/10 21:20]: >>521
F[x_]の定義を良く読め。

それからMathematicaでプログラミングするときは
出力を確認しながらやれ。
F[x_] の行をセミコロン無しで打つか、
F[x]とか試しに打ってみながらやってれば間違えないだろ？

あとは自分でやれ。
523 名前：１３２人目の素数さん mailto:sage [03/07/11 14:51]: Solve[]や、NSolve[]で、求めた複数の解のうち、
ひとつだけを使いたいのですが、
たとえばt = %[[-1]], t = Last[%]などとしても、
{x ->0.234hogehoge}と　t　に入ってしまいます。

この x->を取り除くにはどうすればいいのでしょうか？
524 名前：523 mailto:sage [03/07/11 15:00]: x /. %
でできました。トンクス
525 名前：１３２人目の素数さん mailto:sage [03/07/12 00:43]: >521
えー、 Sum[<<>>,{i,1,100}] ってあるけど
数値やるなら NSum じゃ駄目なん？

どうせFindRoot なんだし
526 名前：１３２人目の素数さん [03/07/14 23:27]: あのー四捨五入するような関数ってないですか？切り捨てでもいいんですけど。
マニュアルで調べてみたんですが見つからなくて・・・・
527 名前：１３２人目の素数さん [03/07/14 23:28]: Roundじゃないの？
528 名前：１３２人目の素数さん [03/07/14 23:36]: 0.5足して切り捨てろ。
529 名前：山崎渉 mailto:（＾＾） [03/07/15 12:38]: 　__∧＿∧_
　|（　　＾＾）|　＜寝るぽ（＾＾）
　|＼⌒⌒⌒＼
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉
　　 ~￣￣￣￣
530 名前：１３２人目の素数さん [03/07/16 23:06]: shop.ponycanyon.co.jp/kids/p14796.html
531 名前：１３２人目の素数さん [03/07/17 07:26]: age
532 名前：１３２人目の素数さん [03/07/17 14:11]: >>530
これなら3,800円か。安いな。
533 名前：１３２人目の素数さん [03/07/17 16:26]: 質問なんですが、例えば
x=k/10
y=k+1
For[k=1,k<11,k++,
　　r=x+y]
みたいに繰り返し計算をさせる場合に、このrは、毎回更新されてしまうんで
その結果を格納してリストの形にしたいとするとどーすればいいですか?
勿論、これぐらい単純だとTable使えばいいんですが、やりたいのはもっと
ややこしいんで、結果をリストにする格納する方法を教えてください。
よろしくお願いします。
534 名前：１３２人目の素数さん mailto:sage [03/07/17 19:13]: >>533
単純であろうがなかろうがTableでいいはずだ。

どうしてもというなら
l ={}
For[k=1,k<11,k++,r=x+y;l = Prepend[r,l]]
等。

しかしmathematicaではFor文を使う必要はない(その上
遅いので使うべきでもない)ので再考することを勧める。
535 名前：１３２人目の素数さん mailto:sage [03/07/18 14:39]: とぇ
536 名前：１３２人目の素数さん [03/07/22 12:25]: FindMinimumを使って最小値とそのときの変数の値を求めると、
{234,{a -> 3,b -> 2}}
みたいな感じだと思うんですが、これをいっぱい出して、a,bの平均値を出すには
どうすればいいでしょうか？
537 名前：１３２人目の素数さん mailto:sage [03/07/23 01:12]: >>546
漠然としすぎ。

こういうことか？

<< Statistics`DescriptiveStatistics`

data = {x, y} /. (#[[2]] &) /@
Flatten[Table[
FindMinimum[Sin[x + y], {x, i}, {y, j}], {i, 0, 3}, {j, 0, 3}],
1];

Mean /@ Transpose @ data
538 名前：１３２人目の素数さん mailto:sage [03/07/23 01:12]: >>537
s/546/537
539 名前：537=538 mailto:sage [03/07/23 01:13]: 536駄・・・まちがえすぎだ。
540 名前：536 [03/07/24 13:07]: すみません、537のプログラム見て解読してたんですけど、ちょっと良くわかんなくて。
たとえば、
FindMinimum[f[a,b],{a,1},{b,2}]
とした場合、{最小値,{a -> 0.89898},{b->0.1}]
みたいにした場合、a->の部分を消すのはどうすればいいですか？
541 名前：536 [03/07/24 14:36]: なんか読みにくいですね・・えーと置換なんかでうまくa->の部分を消せない
かなあと・・・無理なんですかねえ。。。
542 名前：１３２人目の素数さん [03/07/24 14:37]: 開店セール見に来てーーーーーーーーーー！DVDを見ようー！GO！
▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼
開店セール1枚　500円より→→→www.dvd-exp.com　
アダルトDVDが激安販売→→→www.dvd-exp.com　
アダルトDVDの価格はなーんと1枚900円→→→www.dvd-exp.com　
最短でご注文いただいた翌日にお届け！→→→www.dvd-exp.com　
セルビデオ店の人気商品が24時間お買い物が出来ます。
見やすい画面・注文簡単・時間をかけてゆっくり選んでください！
　　　　　　　　　　www.dvd-exp.com←へGO！
　　　　　　　　　　　DVDエクスプレス
543 名前：１３２人目の素数さん mailto:sage [03/07/24 15:18]: >>541
簡単に消せるように"->"つきで出力されている。

s = {a -> 1, b -> 2}
{a,b} /. s
544 名前：１３２人目の素数さん [03/07/24 22:08]: s はリストみたいになっているから
s[[1]] や s[[2]] とかで取り出しても良いんじゃない？
545 名前：１３２人目の素数さん mailto:sage [03/07/25 04:10]: >>544
543で解決してます。”a->”みたいなのを取り除きたいってのが541の質問だろ。
ってか、房は黙ってろ（ｗ
546 名前：536 [03/07/25 12:02]: >>543
ありがとうございました。出来ました。
>>544
それだと、a -> 0.89898などと出ます。
547 名前：１３２人目の素数さん mailto:??? [03/07/25 13:33]: s[[2,2]]
とやるんじゃ
548 名前：１３２人目の素数さん [03/07/25 14:11]: 【クレジットカードのショッピング枠を現金化！】

クレジットカードで買い物をすると、商品代金の８０～９０％のキャッシュバック！！
１０万円の買い物をすれば８～９万円の現金がお手元に！！

■■朝日カード■■
www2.pekori.to/~ryu/asahi/
549 名前：１３２人目の素数さん mailto:sage [03/07/25 21:29]: >>547
自分で一回やってみれば？違うのよくわかるから。
ほとんど使ったことないでしょ？分りもしないのにアドバイスしようとしない
方がいいじゃないのか？
550 名前：１３２人目の素数さん mailto:sage [03/07/25 22:37]: >549

残念ですが{a->1,b->2}[[2,2]] で 2 が出てくるよ。
[[]] は別に行列以外でも使える、数式構造の場所指定でした。

まあこのやり方は、-> の値を取り出すのにふさわしくないとは
僕も思いますがね。
551 名前：１３２人目の素数さん mailto:sage [03/07/25 22:40]: In[4]:=
{a -> 1, b -> 2} [[{1,2},2]]
Out[4]=
{1,2}

List[Rule[a, 1], Rule[b, 2]]

In[6]:=
FullForm[{{1,2},{a,b}}]
Out[6]//FullForm=
List[List[1,2],List[a,b]]
552 名前：１３２人目の素数さん mailto:sage [03/07/25 22:52]: >>550
質問の意図と違うから。その意味での違うなんだけど？なにがでしただよｗ
553 名前：550 mailto:sage [03/07/25 23:02]: >552

確かに s[[2,2]] だと、本当に欲しい {1,2} は
出てきませんね。失礼スマソ。

ただ -> に対して、通常考える /. 以外の [[]] 指定方法に
目を付けた点は悪くはないのでは、とちょっと擁護
してみたかったのでした。
554 名前：１３２人目の素数さん mailto:sage [03/07/26 10:30]: >>553
> 目を付けた点は悪くはないのでは、とちょっと擁護

何がやりたいのかわかりにくい、バグに弱い、
一般的な状況(a^2 + aを求めるとか)では使えない…
センス悪すぎるよ。
555 名前：１３２人目の素数さん [03/07/26 20:47]: リストから第一行目、第十一行目、第二十一行目・・・というふうに取り出したい
んですけど、どーすればいいですか？
556 名前：１３２人目の素数さん mailto:??? [03/07/26 22:46]: >554

Mathematicaのルールどおりだからバグに弱くない
おまえがセンス悪、
それとも正確ﾜﾙか
557 名前：１３２人目の素数さん [03/07/26 23:07]: 　　　　　　　　　　　　　　　 _,.. ----　.._
　　　　　　　　　　　　　　,. '"　　　　　　　｀丶、
　　　　　　　　　　　／　　　　　　　　　　　　` ､
　　　　　　　 ,..-‐/　　　 ...:　 ,ｨ　,.ｉ .∧　,　　ヽ.
.　　　　　　　,:'　　.l　.::;',. :::;/..://:: /,':/　 ', l､ .i　　ヽ
. 　　　　　 ,'　　..::| .::;',' :;:','フ'７フ''7/　　 ',.ﾄ',_|,　, ',.',
　　　　　　　,' 　 .::::::!'''l/!:;'／　/'ﾞ　 /　　　'! ﾞ;:|:､.|､| 'l
. 　　　　 ,'. 　.:::::::{　ｌ'.l/　　､_　　_,.　　　　　　'l/',|.';|
　　　　　　ｌ　 :::::::::::';、ヾ　　　￣　　　 `‐-‐'/! ';. '
. 　　　　 !　:::::::::::/ `‐､　　　　　　　　ゝ　　　|'ﾞ　|
　　　　　　 | ::::::::／　　　＼　　　　､＿,　_.,.,_　ﾉ:::　!　
　　　　　　 |::::／. 　　　　_rl｀': ､_　　　///;ﾄ,ﾞ;:::::./ 　　　マターリと良スレ。
..　　　　　　`´　　　　　 /＼＼　 `i;┬:////ﾞlﾞl ヾ/　　　　　　　　　　数ヲタがうなぎを食べられますように・・・
　　　　　　　　 ,.:く::::::::`:､＼　〉lﾞ:l 　/　!.|
.　　　　　　　　／:.:.:.:＼:.:.:.:.`:､ｿ/:.:| 　　| |
　　　　　　　　 /.:.:.:.:.:.:.:.:.:＼:.:.:.:У:.:;l　　 /./
.　　　　　　　 /:.:.:.:.:.:.:.ｒ'´`‐,`､:/.,.:‐{　　 |　!`:､
　　　　　　 ,'.:.:.:.:.:.:.:.:.';_,ﾟ.,ﾉ.:./,:':.:.:.:',　　|　|`､:|
　　　　　　 !:.:.:.:.:.:.:.:.:.:.ﾞ､:.::／:.:.:.:.:.:.ヽ, /　,!:.:`､
558 名前：１３２人目の素数さん mailto:sage [03/07/27 00:52]: 鰻食いてーー
559 名前：１３２人目の素数さん mailto:sage [03/07/27 00:53]: でも貧乏だからとても買えねーよ(;´Д`)
560 名前：１３２人目の素数さん mailto:sage [03/07/27 01:07]: >>556
> Mathematicaのルールどおりだからバグに弱くない
あきれますな。悪いこと言わんからもっと謙虚になれ。
プログラムくんだ経験自体無いだろ？

たとえばFindRootで解が見付からなくてそのまま式自体が
かえったらどうなる？

[[]]ってのはかなりError proneだぞ。
561 名前：１３２人目の素数さん mailto:sage [03/07/27 09:37]: >555
tensor_a[..., Table[1+10 i,{i,0,5}], ....]
てな感じ。ポイントは添字指定で数字ではなく
リストを使うこと。
562 名前：550 mailto:sage [03/07/27 09:53]: >560
/. でも FindRoot そのままの式が返ってくると
あまりうれしくない気がするけど。。。

/. を使うと a^2+a /. みたいな純関数的使い方が出来て
見通しが良くなるのは確かに同意。
つまり /. と -> はペアで使うのが基本なのは当然だけど
Lisp 的な感覚として式の木構造を考えることが出来る [[]]も
無下に NG とするのもどうかも。。。。

あきらめ悪いっすかね
563 名前：１３２人目の素数さん [03/07/27 10:59]: Mathmatica 4を使っています。
Module[{i, j},
　ans = {};
　Do[
　　Do[
　　　ans = Append[ans, {i, j}], {j, i + 1, 3}
　　], {i, 2}
　]
　ans
]
を実行すると
{{Null, 2 Null}, {Null, 3 Null}, {2 Null, 3 Null}}
と変なNullという文字が引っ付いてくるのですが、
どうすれば、まともな答えが戻ってくるのでしょうか？
564 名前：１３２人目の素数さん mailto:sage [03/07/27 13:15]: >>562
エラーメッセージが出て、文字がそのままかえってくる
のと、
「何も言わずに、間違った答えがかえってくる」
のでは大違いだぞ。

全ては式である、ということの一例として[[]]がつかえることを
理解していることは必要だが、無闇に使うのは否定すべきだよ。
565 名前：１３２人目の素数さん mailto:sage [03/07/27 13:19]: >>563
　]
　ans
->
　];
　ans
566 名前：１３２人目の素数さん mailto:sage [03/07/27 14:01]: >>563
おれは>>534だから一応Tableによる解決を書いておく。
このほうがすくなくとも数倍速い。

UpperInds[n_Integer] := Flatten[Table[Table[{i,j},{j, i + 1, n}],{i,n - 1}],1]
567 名前：563 [03/07/27 14:45]: >>565
ありがとうございます。解決ししました。

>>566
ありがとうございます。
じつは、「内側のループのなかで、i,jを使った条件判定をやった結果、条件
を満たした時のみ、i,jに依存する要素をansに加えていく」というのがやりたい
ことなので、ちょっとそのままでは教えていただいたやりかたは使えないよう
です。
568 名前：１３２人目の素数さん mailto:sage [03/07/27 15:15]: >>567
それなら
UpperInd[n_Integer] := Module[
{ConvertFunc = Sin[#[[1]] + #[[2]]] &, CondFunc = EvenQ[#[[2]]] &},
ConvertFunc /@ Flatten[Table[
Select[Table[{i, j}, {j, i + 1, n}], CondFunc], {i, n - 1}], 1]]
あるいは、
Table[If[..., ..., Gomi],...]としてあとでGomiを除くとか。

綺麗なのは、

ConvertFunc /@ Select[Flatten[Table[{i,j},{i,n},{j,n}],1], CondFunc]

(CondFuncのなかで#[[1]] < #[[2]]も判定)
みたいなやつかな。
569 名前：１３２人目の素数さん [03/07/29 01:55]: すみません、ものすごく初心者なんですけど教えてください。
非線形の最適化ってできるのですか？
お願いします。
570 名前：１３２人目の素数さん [03/07/29 11:57]: ウルフラムさん
セルオ－トマトンの原理で宇宙の起源を解明しようとしてるらしい
一体どうしちゃったんだろうか
571 名前：１３２人目の素数さん mailto:Sage[] [03/07/29 11:58]: なんかMathematica5.0ての海外では売ってるらしいんですけど。
572 名前：１３２人目の素数さん mailto:sage [03/07/29 22:13]: a[0]=0
a[1]=1
a[n+2]=a[a+1]+a[n]
みたいな漸化式を使ってa[100]とかを求めるにはどうすればいいですか？
573 名前：１３２人目の素数さん [03/07/29 22:39]: age
574 名前：１３２人目の素数さん mailto:sage [03/07/29 22:58]: >>572
a[0] = 0;
a[1] = 1;
a[n_] := a[n - 1] + a[n - 2]
575 名前：574 mailto:sage [03/07/29 23:00]: >>572
100は計算回数(2^100)が多すぎる求められんと思うぞ。
計算結果を覚えながらやれば別だが:
b[0] = 0;
b[1] = 1;
b[n_] := (b[n] = b[n - 1] + b[n - 2])
みたいに。
576 名前：１３２人目の素数さん [03/07/29 23:04]: レスありがとうございます！
早速やってみます
577 名前：１３２人目の素数さん [03/07/29 23:06]: 出来ました！すごく感謝！
578 名前：１３２人目の素数さん [03/07/30 02:25]: x+y=1
xy=-1
というとても簡単な連立方程式を解かせたいのですが、
Solve[{x + y == 1, x - y == -1}, {x, y}]
とやっても上手くいきません。どうすればよいでしょうか？
579 名前：１３２人目の素数さん [03/07/30 02:26]: >>578
ミスったｗ
Solve[{x + y == 1, x y == -1}, {x, y}]
これでやって出来なかったです。お願いします。
580 名前：１３２人目の素数さん mailto:sage [03/07/30 10:56]: >>579
できるけど。
581 名前：１３２人目の素数さん [03/07/30 11:47]: 買ったけど使い方がさっぱりわからない。友達に売ろ。
582 名前：１３２人目の素数さん [03/07/30 11:53]: ムカついたんでnyに流そっかな・・・
583 名前：１３２人目の素数さん [03/07/30 14:31]: >>575
b[1244]まではそれで出るんですけど、
b[1245]以降はHoldなんたらとか出て計算できないんですが
どうすればもっと大きい数まで表示出来るようになるでしょうか？
584 名前：math.1st ◆M9pCfogc9g [03/07/30 14:51]: Re:>583 再起呼び出しの関数は、ネストできる数に制限があることから、使わない方がよい。
使用メモリを増やして計算すると大きい数も表示できて、速い。
585 名前：１３２人目の素数さん [03/07/30 15:01]: >>584
ありがとうございます。
もうひとつ質問させてください。
n≦10の時
a[n_]:=a[n-1]+2
n＞10の時
a[n_]:=a[n-1]+1
のように、範囲を指定したい時はどうすれば良いでしょうか？
586 名前：１３２人目の素数さん mailto:sage [03/07/30 21:34]: >>585
そんな問題Cでやった方が速い。2^32進数の作り方もついでに勉強。
まったく典型的「ユーザさん」だね。
あなた甲南大学とやらの学生ですか？
そういう丸投げ的利用法は、あとあとためにならないですよ。
587 名前：585 [03/07/30 23:42]: >>586
いえ、高校生です。
588 名前：１３２人目の素数さん mailto:sage [03/07/30 23:55]: 高校生ならいいってもんでもないな。
589 名前：１３２人目の素数さん mailto:sage [03/07/31 01:49]: >>583
末尾再帰にすれば良い(演習問題)。
関数型言語では、関数は全て再帰で書く。

末尾再帰に書き換えないのであれば、
$RecursionLimit=3000
などと再帰の上限階数を増やす。

末尾再帰については自分で調べろ。

>>584
くだらんレスをするな。

>>585
例えば:
a[n_] /; n <= 10 := a[n - 1] + 2
a[n_] /; n > 10 := a[n - 1] + 1

>>586
丸なげ反対はわかるが、
> のように、範囲を指定したい時はどうすれば良いでしょうか？
これは言語仕様について聞いているんだろう。
590 名前：１３２人目の素数さん mailto:sage [03/07/31 23:36]: >>585
その後 FullForm で、出てくる式の構造を暇な時に見ると
将来、ひょっとしたら役に立つ可能性があるかもしれない。
591 名前：585 [03/08/01 02:39]: >>589
ありがとうございます。
末尾再帰ってのは調べてみたんですけどよく分かりませんでした…。
n = 2000; a = 0; b = 1; While[n >= 1,
If[n ≠ 1, a = a + b; n--, Print[b]; Break[]];
If[n ≠ 1, b = a + b; n--, Print[a]; Break[]]]; If[n == 0, Print[0]]
結局上のを書いてみたんですが改良点が有ったらご指摘お願いします。

f[n_]:=a = 0; b = 1; While[n >= 1,
If[n ≠ 1, a = a + b; n--, Print[b]; Break[]];
If[n ≠ 1, b = a + b; n--, Print[a]; Break[]]]; If[n == 0, Print[0]]
あと、上のようにやってf[2000]とかやっても0とか出ちゃうんですけど、
何が間違ってるのでしょうか…。
592 名前：１３２人目の素数さん mailto:sage [03/08/01 11:58]: 高校生でこんなことをやってるのか。
凄いな
593 名前： mailto:sage [03/08/01 12:39]: え～今までMathematicaなんて使ったことないのに、卒研で使えと言われまして・・・
さっぱりわからないんで教えていただけないでしょうか？

ワイブル分布において
平均Mean[dist] == 1と標準偏差StandardDeviation[dist] ==1/3の連立方程式を
αとβについてSolveで解こうとしたんですが、
『逆関数が使用されています．値は多価逆関数のため失われている可能性があります．』
と出て解けないんですよ。

これってどうすればいいのでしょうか？
594 名前：１３２人目の素数さん [03/08/01 16:39]: Play[Sin[700 t + 25 t Sin[350 t]], {t, 0, 4}]
595 名前：１３２人目の素数さん mailto:sage [03/08/01 16:53]: >>593
解けないんじゃなくて、そのメッセージの後に代表値は出てるよな？（ｗ
596 名前：１３２人目の素数さん mailto:sage [03/08/01 17:07]: Mathematicaとは関係ない話なんだけどwolfram reserchの日本語の
トップページ(www.wolfram.co.jp/)ってさ、リロードするたびに違う
種類の画像とキャッチコピーが表示されるじゃん。キャッチコピーは
「世界の計算方法」とか「数学、科学、工学、教育における革新的言語」
とか色々あるんだけど。あれの中の一つに「技術計算の未来を定義」
ってのがあったのね。それを読んだときに「未来を予測する最良の方法
は、それを発明してしまうことだ」っていうAlan C. Kayの言葉を思い出し
て、「ああそうか、数学で何か発明するってはまさしく何か定義すること
だな」と思って少し感慨を覚えたんだよ。どうでもいい話だけど。

ってか英語の方のトップページ見たらMathematica5って書いてあるじゃん！
6/23に出たのか。全く知らなかった。
597 名前：１３２人目の素数さん [03/08/01 17:31]: 例えば、
p1=Plot[Sin[x],{x,-10,10}];
というように最後にセミコロン「；」を付けてもグラフが表示されてしまうのですが
これを表示させないようにするにはどうすればよいのですか？
598 名前：ビッグバン宇宙論は完全に間違いでした！！！！！ [03/08/01 18:39]: 科学者よ、恥を知れ！
ビッグバン宇宙論は完全に間違いだった！
科学の原則を無視した、デタラメのインチキ理論だったのだ。
そして、そのビッグバン宇宙論の世界的な浸透は
アメリカ、ユダヤ・キリスト教勢力による世界支配のための思想的な戦略なのだ！
また、ビッグバン宇宙論の思想によって戦争が起こり、
貧富の差がひらき、終末的な絶望感が世界に蔓延しているのだ。
ビッグバン宇宙論は世界の平和を揺るがす、悪の元凶となっているのだ。
ビッグバン宇宙論とは、
「宇宙は『無』からビッグバン（大爆発）によって誕生した」という理論である。
この理論は、ユダヤ・キリスト教の創造神話（神が天地を創造した）そのものである。
ビッグバン宇宙論の実態は、科学理論ではなく宗教思想なのである。
『無』は科学的に証明できるものではなく、
そして、『無からの誕生』も科学では証明できるものではないのだ。
ビッグバン宇宙論が科学の正統であるという思想を、世界中の人々に
浸透させる戦略が成功したことにより、ユダヤ・キリスト教勢力の
世界における優位性が確立されていったのだ。（２０世紀に）
そして、その思想的支配の最大の例が、アメリカやイギリスによる
イラク戦争なのだ。
ビッグバン宇宙論の浸透により、世界中に終末思想（世界の終わり）が蔓延してしまっている。
そのことにより、自己中心的、せつな的、短絡的な考え方が社会に広がっている。
科学的に間違っているビッグバン宇宙論から脱却しなければならない。
そして、宇宙は無限だということを理解しなければならない。
人間は本当の宇宙観、世界観を構築し、
新しい時代に進んでいかなければならないのだ。
ビッグバン宇宙論が世界を支配している限り、平和な世界にはならないのだ。
そのことを科学者は重く受けとめるべきである。
平和な世界を！！！！！！！！！！！！！！！！！！！！！！！
599 名前：１３２人目の素数さん [03/08/01 18:49]: 　　　　　　　　　　　,-rヽ、　　　　　　 __,.- 、
　　　　　　　　　　 /　）ｧﾍ　　　　　 / 　　 } 　
　　　　　　　／ ,ｨ　,.ﾆｨ　　　　 /´ヽo‐´
　　　　　　　／＿〈ﾚ'´∠-‐‐- ､__　/　　ｌ
　　　　ﾄ<,__,＞　 _ ,,,,, _　　　＜/　　　|
　　　　（＼_ ､__∠-､‐､‐ヽ　　､ .l　ｌ　 /
　　　　　〉／ _.ィ‐'´ヽﾄ､Alヽﾄ､ヽｌヽl　　　|
　　　イァﾉ.ｲ'"_　　　 ;==､ヽ　ｌ　　 L._
　　　　,ｲ |　　ｌ /lﾊ　　　}Viゝ } /|　　　ｌ:.:.ヽ
　　　／.:L.ヽﾊﾄﾞ上!　　　ﾞ‐'　ﾉﾄ､　　レ' {
　　〈-‐'´.-'´ >ゝ" ` _　....、´ ´7/､｀ー'´　.ｲ
　　　ヽ‐'´,-／ィヘ、　ヽ _ﾉ　　,彳　｀ｰ-‐'´ /
　　　　｀´l　　ﾉィﾍｨ＞､__,.. ‐'´ {　 l.　　　/
　　　　 l　 ! /〈__.L..､ﾉ __.-‐'´| 　! 　　〈
. 　　　　l　! ＞'´-､イ⌒'ｰ-‐'´　　、　　/
　　　　　　 | /＿,ィ´i　l. . : : : : . . ,' 　｀ｰ./
　　　　　〉､: : :.l _,イ: : : : : : : : !　　 /　
　　　　　　/. :! : 　´　　 : : : : : :/　　 / 　
　　　　　l: : : 　　　　　　　　　 /　
　　　　　 /. :　_. -‐- ._　　　. : : : :.〈
　　　　／_,.ィ´　　　　　｀ｰ-､;_: : : :|
　　　＼::ﾉ　　 j　__.　-‐ ‐- エ＞
　　　　　 /‐ '' "__,... -‐‐ｒｰ､-､_|　　
　　　　／-‐ｧ' ´/　 |　　l　ヽ　 ` 、
　　 .／　／　 /　　l 　　! 　`、　　ヽ、
600 名前：１３２人目の素数さん mailto:sage [03/08/01 19:29]: おい、そのAAはロリスレに貼れ
601 名前：ぼるじょあ ◆yBEncckFOU mailto:（＾＾） [03/08/02 03:01]: 　　　　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　・３・） (　　＾＾）＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕
　　＝ ◎――――――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ
602 名前：１３２人目の素数さん mailto:sage [03/08/02 03:41]: >>591
> n = 2000; a = 0; b = 1; While[n >= 1,
> ...
論外なのでコメントできん。一冊以上Mathematicaの本を読め。

末尾再帰を用いた解答
f[n_] := f[n, 0, 1];
f[0, r1_ , r2_ ] := r1;
f[n_, r1_, r2_] := f[n - 1, r2, r1 + r2];

> f[n_]:=a = 0; b = 1; While[n >= 1,
a = 0;のセミコロンで文が切れている。
f[n_] : = a
b = 1;
と入力しているのと同じだ。
603 名前：１３２人目の素数さん mailto:sage [03/08/04 00:48]: >>602
論外ですかｗ
そうゆう本高そうで…。今度学校の図書館で詠んでみます。
ありがとうございました。
604 名前：１３２人目の素数さん [03/08/04 08:23]: 597に応えてあげて
605 名前：１３２人目の素数さん mailto:sage [03/08/04 20:18]: >>604=>>597
606 名前：１３２人目の素数さん mailto:sage [03/08/04 23:22]: >595
日本国内向けに日本語文書作業中。
もっとも買う人がいないから強制バンドル
価格Up 中。

日本は非英語で高度な技術内容が学べ得るけど
その分仲介業者が入る余地が多くなってしまった。。。。
607 名前：１３２人目の素数さん mailto:??? [03/08/05 23:48]: Ver.5使ってます。
起動直後にガイドが出るくらいで
見た目はそんなに変わってない。
608 名前：１３２人目の素数さん [03/08/05 23:49]: ★おまんこナビ登場！！★日本の美しきオマンコ★
endou.kir.jp/marimo/link.html
609 名前：１３２人目の素数さん [03/08/05 23:50]: あなたの近所のあの娘が出てるかも
いつもすまして歩いてるあの娘が
乱れまくり！！
www.j-girlmovie.com/
610 名前：１３２人目の素数さん mailto:sage [03/08/06 15:32]: Sin[a°]=1/2なるaを求めるにはどうすればよいでしょうか？
Solve[Sin[a°]==1/2]では求められませんでした。
611 名前：１３２人目の素数さん mailto:sage [03/08/07 00:11]: Solve[Sin[a] == 1/2, a]
612 名前：１３２人目の素数さん mailto:sage [03/08/07 14:32]: >>611
違う。
613 名前：１３２人目の素数さん mailto:sage [03/08/07 14:59]: >>612
Solve[Sin[a Degree] == 1/2, a]
で
Solve[Sin[a] == 1/2, a]
と同じ程度に、求めることができているので、
質問者が何を求めているのかよくわからない。

Solve[Sin[a Degree] == 1/2, a] // N
ってのが欲しいのか。
614 名前：１３２人目の素数さん mailto:sage [03/08/07 15:25]: >>613
スペースが必要だったのか。ありがとう。
615 名前：１３２人目の素数さん mailto:sage [03/08/09 12:52]: >>614
そういう誤解か。
a\[Degree]はそれ自体がシンボルだから、
Solve[Sin[a\[Degree]] == 1/2, a\[Degree]]
はできるが、aは求められない。

スペース(かけ算)が見た目に嫌なら、
a\[InvisibleTimes]\[Degree]とすればよし。
616 名前：１３２人目の素数さん mailto:sage [03/08/09 15:56]: >>597
だれも書かないな。
$DisplayFunction = Identity
617 名前：１３２人目の素数さん [03/08/12 22:43]: >616
正確にいきましょう。

p1 = ListPlot[Table[{i, i*Random[]}, {i, 1, 8}],
PlotStyle -> RGBColor[1, 0, 0], DisplayFunction -> Identity];
p2 = Plot[Exp[x], {x, -2, 2},
PlotStyle -> RGBColor[0, 0, 1], DisplayFunction -> Identity];
Show[{p1, p2}, DisplayFunction -> $DisplayFunction];
618 名前：１３２人目の素数さん mailto:sage [03/08/13 00:54]: >>617
それのどこが正確なんだ？
PlotStyleなんて何の関係もないし。
# 意味がわかってないままコピペしているようにしかみえないのだが。

で、これのどこが正確でないんだ？
$DisplayFunction = Identity;
Plot[Sin[x],{x,0,1}];
# オプションの方がいい場合「も」あるがな。
619 名前：１３２人目の素数さん mailto:sage [03/08/13 01:56]: >>617-618
要はオプションを関数の中で指定しているかグローバルに指定しているかの違いだろ。
620 名前：１３２人目の素数さん mailto:sage [03/08/13 03:57]: >>619
違う。
$DisplayFunction は DisplayFunctionというオプションのデフォルト設定
を与えるものであって、オプションそのものではない。
オプションを設定することと、
(オプションが省略された時の)デフォルトを変更することは、厳密には別のことだ。

オプションを関数の外で設定するというのは
SetOptions[Plot, DisplayFunction -> Identity]
だ。
621 名前：山崎渉 mailto:（＾＾） [03/08/15 18:32]: 　　　 (⌒V⌒)
　　　│ ＾＾ │＜これからも僕を応援して下さいね（＾＾）。
　　⊂|　　　　|つ
　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン
622 名前：１３２人目の素数さん mailto:??? [03/08/18 18:59]: おーい、Ver5持ってる奴
x^100-1==0
を解いて、根の絶対値を調べてくれ。
１にならないから。
623 名前：tebe [03/08/18 19:06]: 半年間もかけてMathematicaの使い方を学ぶという
非常にエグイ授業を大学でやったよ。
Mathematicaって機能ありすぎて、僕みたいな素人
には宝の持ち腐れって感じだな
624 名前：１３２人目の素数さん [03/08/18 19:11]: ｢今私のオッパイ見てたでしょう。｣ってそりゃ見ますって。
これだけの立派なオッパイが目の前にあったら。
淫乱女教師がフェラ、パイズリで同僚の先生を責めまくります。
挿入部のアップも豊富なのでオマンコもモロ見えですよ。
無料動画をどうぞ
www.exciteroom.com/
625 名前：１３２人目の素数さん [03/08/18 19:19]: Mathematicaの有名なバグ

Integrate[1/x^2, {x, -1, 1}]

古いバージョンだと２と評価してくれるんですが
新しいバージョンだと∞になってしまいます。
Ｖｅｒ５だとどうなんでしょうか？
626 名前：１３２人目の素数さん mailto:sage [03/08/18 21:50]: >>625
はぁ？
バグってるのはあなたの頭では？
627 名前：１３２人目の素数さん mailto:sage [03/08/20 01:54]: 　1　　　　　　　　1
∫(1/x^2)dx＝-1/x｜＝-1/1-(-1/(-1))＝－２
　-1　　　　　　　-1

じゃないの？
628 名前：１３２人目の素数さん [03/08/20 01:56]: >>627
ずれた
　1　　　　　　　　　　1
∫(1/x^2)dx＝-1/x｜＝-1/1-(-1/(-1))＝－２
　-1　　　　　　　　　-1

じゃないの？
629 名前：１３２人目の素数さん [03/08/20 01:59]: >>628
1/x^2>0だから積分値は正だろ？
630 名前：_ mailto:sage [03/08/20 02:03]: homepage.mac.com/yamazaki8
631 名前：１３２人目の素数さん [03/08/20 02:06]: >>625-629
ｘ＝０で発散しているから∞だろ
632 名前：１３２人目の素数さん mailto:??? [03/08/20 08:37]: 625がまぬけだったということでよろしいでつか
633 名前：１３２人目の素数さん [03/08/20 10:12]: >>632
否、これは確かに古いバージョンでは２ではないが－２と評価しているぞ。
いまのバージョンだと∞に修正されてはいるが。
634 名前：supermathmania ◆ViEu89Okng [03/08/20 10:50]: 発散積分の有限部分は確かに-2になる。
Integrate関数は超関数積分なども一部サポートしている。
これとは関係ないが、主値積分もオプション指定でできる。

ついでに、Integrate[Exp[1/(1-x^2)],{x,-1,1}]と
NIntegrate[Exp[1/(1-x^2)],{x,-1,1}]の結果を報告していただきたい。
期待しないで待っているぞ。
635 名前：１３２人目の素数さん [03/08/20 18:05]: >>625は古いバージョンと新しいバージョンの挙動の違いを言ってるだけだろ
積分値が正しい間違ってるの話をしているわけではないと思われ
636 名前：132人目の素数さん [03/08/20 19:54]: >>634
NIntegrate[Exp[1/(1-x^2)],{x,-0.998,0.998}]=3.88057*10^103でつた。
Integrate[Exp[1/(1-x^2)],{x,-1,1}]は、答えが帰ってきません（；＿；）
mathematica 4.2 for Linux-PPC
637 名前：１３２人目の素数さん mailto:??? [03/08/20 20:50]: ちなみにMapleの結果

int(exp(1/(1-x^2)),x);

/infinity \
| ----- 2 |
| \ hypergeom([1/2, _k1], [3/2], x )|
x | ) --------------------------------|
| / GAMMA(1 + _k1) |
| ----- |
\_k1 = 0 /
638 名前：１３２人目の素数さん [03/08/21 01:41]: オクで割れマセマティカを5000円で売ってるやつがいる。
タイーホされてほしい。
639 名前：１３２人目の素数さん mailto:sage [03/08/21 01:55]: >>638
通報すれ
640 名前：１３２人目の素数さん [03/08/21 02:06]: >>637
そんなの超幾何関数使って元の積分を表現しただけじゃん。
641 名前：１３２人目の素数さん mailto:sage [03/08/21 03:33]: 0.11111111111111111111111111
って入力して実行すると
0.11111111111111111111111111
って出て、それをクリックしてエンター押すと
0.1111111111111111111111111099999999999964597
って出るんですがどうゆう意味があるのでしょうか？かなり気になるんですが。
642 名前：１３２人目の素数さん [03/08/21 03:40]: 意味はない。Mathematicaが糞なだけ。（ｗ
643 名前：１３２人目の素数さん mailto:sage [03/08/21 05:45]: >>641
ネタで言ってるのかどうか知らんが・・・

機械精度の演算がどういう仕組みで行われているかを一度考えて見れ。
644 名前：１３２人目の素数さん mailto:sage [03/08/21 07:45]: それにしてもクソですな。そのくらい改善できそうなものだが。
645 名前：supermathmania ◆ViEu89Okng [03/08/21 15:58]: Integrate[Exp[1/(x^2-1),{x,-1,1}]
NIntegrate[Exp[1/(x^2-1),{x,-1,1}]
と云おうとしたら、間違えてしまった。
これなら計算できるのかな？
646 名前：１３２人目の素数さん mailto:sage [03/08/21 19:01]: >>644
あのさ、Mathematica では小数点を付けた数値は機械精度
で計算されるとマニュアルに書いてあるから、
0.11111111111111111111111111
と書いた時点でこれは
11111111111111111111111111/100000000000000000000000000
とは全然違う値になるわけだ。だから、もし
0.11111111111111111111111110000000000000000000
とかになったら、Mathematicaは機械精度の計算を唄っておきながら
実際にはサポートされていないことになる。

この前も同じような質問があったけど、今の大学って数値計算の理論
とか教えないんだろうか。もし習ったことが無いなら、そういう技術者や
科学者にとって必要最小限のことも教えない大学なんて行く意味がないよ。
>>641や>>642、>>644が通っている大学の学長は、こういう不祥事（と
言っても良いよね）に対する記者会見を開いても良いぐらいだと思う。鶴を
燃やした事件みたいにさ。高校生だったら、まあ仕方ないけど。

いや、今日は社会人になっても、２次方程式の解を計算するルーチンで
解の公式をそのまま使うような馬鹿に遭遇したんで、イライラしてるのかも
しれん。桁落ちが起きるだろ、アホが！お前の先輩の文型出身の奴でも
ちゃんと場合わけして回避してたぞ。
647 名前：１３２人目の素数さん mailto:??? [03/08/21 20:13]: NIntegrate[Exp[1/(x^2-1)],{x,-1,1}]
=0.443994

ですな。
648 名前：１３２人目の素数さん mailto:??? [03/08/21 20:14]: 数値解析の講義なんか、あったって選択だろ。
そんなもんやってる時間があったら、
純粋数学やるわいな
数学専攻ならな
649 名前：１３２人目の素数さん mailto:sage [03/08/21 20:35]: >>646
俺は高校生です。仕方ないです。
650 名前：１３２人目の素数さん mailto:sage [03/08/21 22:29]: >>646
> あのさ、Mathematica では小数点を付けた数値は機械精度
> で計算されるとマニュアルに書いてあるから、

こんなことは書いてない。
おまえの書いた
0.11111111111111111111111111
はArbitrary-Precision Numberだ。
つまり、多倍長数の内部表現が二進であるか十進であるかという実装の問題だ。
偉そうなこと書くんだったらもっと勉強しろ。
一目見たらdoubleで書けないことぐらいわかれよ。

64bit以上のCPUつかってんだったら機械精度かも知れんがどっちにしろ
嘘だしな。
651 名前：650 mailto:sage [03/08/21 22:34]: 一応書いとくが、二進で結構だと考えてるよ。
# 上の馬鹿といっしょにされてはかなわん
652 名前：１３２人目の素数さん [03/08/22 05:18]: このスレ、人を罵倒するやつが多いな。
653 名前：644 mailto:sage [03/08/22 06:20]: >>646
おいおいパターソン&ヘネシーくらい読んだことあるぜ。
0.1111111111111111111111111 != 111111111111111111111/1000000000000000000
(桁は適当ね)になっている実装がクソだと言ってるの。
機械精度なんて本来は一切要らないんだから。
そんなものこの時代のマシン特有の都合だろ？
もちろんMathematicaブックを読めば$MachinePrecisionや
N[]や12.3`なんかについてあれこれ載っている(はず)。
しかし、精度という基本的なことがらに付いてあれこれ頭を悩ませる実装および
言語体系が「クソですな。改善できそうなもんだが」というわけ。

掲示板の書き込みで憂さ晴らすのは不健康だね。
654 名前：１３２人目の素数さん [03/08/22 06:49]: よくわからんがマスマティカはくそってことでつか？
655 名前：１３２人目の素数さん [03/08/22 06:50]: >>653
６４６は書きすぎだと思うが、

> 機械精度なんて本来は一切要らないんだから。
> そんなものこの時代のマシン特有の都合だろ？

それは無いだろー
656 名前：１３２人目の素数さん mailto:sage [03/08/22 11:57]: >>646
知りもしないものを「クソ」あつかいか…
それに少し上のレスぐらい読め。
機械精度なんて関係ないんだよ。

計算機の上では有限の精度しか実現できないことをしっているよな？
0.1111111111111111111111111 == 111111111111111111111/1000000000000000000
にしたところで何か計算したとたんにほぼ確実に10進数表示での
数学的な等号は成立しなくなるんだぞ。
何か意味があるか？
精度の保証だけしとけば十分だろ。

Precision, Accuracyについてマニュアルを読め。
657 名前：１３２人目の素数さん mailto:sage [03/08/22 12:09]: >>646
この時代のマシンのことを無視してプログラム書いてどうすんだよ。

おまえが書くプログラムってどれも実用には耐えないんだろうな。
658 名前：１３２人目の素数さん mailto:sage [03/08/22 20:32]: 私は高校生なので優しくしてください。
0.1111111111111111111111111を任意精度として扱ってもらう簡単な方法はないですか？
いちいち111111111111111111111/1000000000000000000のように入力しなおさいと正しい結果が得られないのは不便なんですけど。
659 名前：644 mailto:sage [03/08/22 21:53]: さあ香ばしくなってまいりました。これも5発売フィーバーか？

>>658 に優しいアシストをすると、そういう小数でどんな計算したいのかを言えば、
より有益なレスがつくのではないかな？ていうか小数なの？有理数なの？
ちなみに手もとのMathematicaブックには「0.142が近似値でなく
厳密な値であるならば、142/1000と有理数として入力しなければならない」
という記述をハケーン。
あと0.11111111111111111111`200 で200桁の有効数字が、
0.1111111111111111111111''200で小数点以下200桁までの正確さが保証されるんだと。
それで十分かな？どうよ？

>>656
> にしたところで何か計算したとたんにほぼ確実に10進数表示での
> 数学的な等号は成立しなくなるんだぞ。
おいおいそれこそ「計算」次第だろ。アタマ疲れてない？
>>657
> おまえが書くプログラムってどれも実用には耐えないんだろうな。
CやF90で書くのとMmaで書くのとには違いが出て当然だろ。言語3～4種しか知らん人？
660 名前：１３２人目の素数さん mailto:sage [03/08/22 22:35]: >>659
ありがとう。
661 名前：１３２人目の素数さん [03/08/22 23:36]: こんなところでしかストレス発散できないやつ。
662 名前：１３２人目の素数さん mailto:sage [03/08/22 23:55]: >>659
つっこみどころ満載だけどおまえにはどんな反論しても
無駄みたいだからやめとく
663 名前：１３２人目の素数さん mailto:sage [03/08/23 00:51]: ＞つっこみどころ満載だけどおまえにはどんな反論しても
というと
＞さあ香ばしくなってまいりました。これも5発売フィーバーか？
の辺りか？(w >659のカキコそのものはべつに間違った内容はないんじゃないの？

ていうか5出るのか。洩れのマスマティカブック4.2版どうなるんだろう。
5対応版買った人に～とか＋'の項目だけ写させてもらおうかな。
みなさんどうしますか？ちなみに洩れはサイトライセンスで使ってる人ね。
664 名前：１３２人目の素数さん mailto:sage [03/08/23 00:53]: なんでこんなにｲﾀｲんだろう・・・
自覚ないのか？
665 名前：１３２人目の素数さん mailto:sage [03/08/23 01:17]: >>663
大分前に出たかと
666 名前：１３２人目の素数さん [03/08/23 02:48]: Mathematicaを持ってない貧乏人がひがんでるスレってここでつか？（ｗ
667 名前：１３２人目の素数さん mailto:sage [03/08/23 14:08]: >>659
>> にしたところで何か計算したとたんにほぼ確実に10進数表示での
>> 数学的な等号は成立しなくなるんだぞ。
> おいおいそれこそ「計算」次第だろ。アタマ疲れてない？

だからほぼ確実に、ってかいてんだよ。日本語ぐらい読めるようになろうね。
割算もできないしSinやCosも計算できないし、かけ算したらどんどん
桁が増えてくし、どうすんの？

> CやF90で書くのとMmaで書くのとには違いが出て当然だろ。言語3～4種しか知らん人？

Mathematicaの設計の話をしてんだよ。

>>660
0.1111111111111111111111111は>>656で書いてあるように自動的に任意精度と
してあつかわれている。
例えば、
x = 0.1111111111111111111111111
y = 1000 x - 111
Precision[x]
Precision[Sin[x]]
Precision[y]
とか入力すると、Mathematicaが精度保証して計算してるのがわかるよ。

高校生だからわからないのだろうが、>>644はアホだぞ。
668 名前：１３２人目の素数さん mailto:sage [03/08/23 15:31]: なんかこのへんのフレームは、任意精度と厳密数値(Schemeでいうexactnum)との
混同に基づいているような気がする。
ということで、高校生君に
> 0.11111111111111111111111111
> って出て、それをクリックしてエンター押すと
> 0.1111111111111111111111111099999999999964597
> って出るんですがどうゆう意味があるのでしょうか？かなり気になるんですが。
の理由を分かりやすく教えてあげれ > 誰か
多倍長整数や高精度近似数は2^16または2^32進数(整数要素の配列)らしい。
あとMathematicaBookの「3.1数」には
123.4 機械精度
123.4567890123456789 システムによっては任意精度 <- なんじゃそりゃ
123.4567890123456789` 機械精度
123.456`200 200桁の任意精度 <- precision
123.456``200 200桁の任意確度 <- accuracy
1.234*^6 機械精度
1.234`200*^6 200桁の精度 <- 任意精度ととっていいのか？
とあるね。しかしこの例が微妙に分かり辛い...
リファレンス「入力の文法」にはdigits.digitsは「近似数」としか書かれてない。
これは任意精度なの？機械精度？だれか教えてくれますか？
わたしは普段整数しか扱わないので知らないままだった。
669 名前：１３２人目の素数さん mailto:sage [03/08/23 16:23]: 0.111... を 1/2^n の和で作ってみると一番近い数が
そのゴミがついた数字になる、というのが一番簡単な
説明だよね。
まあいくつ項を取るかというのは気になる点だろうけどね。
670 名前：１３２人目の素数さん mailto:sage [03/08/23 16:54]: どうも高校生です。
1-0.999999999999=9.99978*10^-13
ってなるのが嫌なんです。これは分数で書くしかないですよね。
999999999999/1000000000000って書いたりするには桁数数えなければなりません。
かなり苦痛なので分数に書き換えるプログラムを書いてみました。
For[i = 0, x ≠ IntegerPart[x], x *= 10; i++]; Rationalize[x]/10^i
xに厳密な数値に直したい小数を入れて、実行すればOKです。
671 名前： mailto:sage [03/08/23 19:37]: ある演算をさせて結果がたとえば
If[x>0,・・・,～～～]
のようになったとき、・・・の部分だけをある変数に代入するにはどうすればいいのでしょうか？
672 名前：１３２人目の素数さん mailto:sage [03/08/23 20:33]: >>668
> 123.4567890123456789 システムによっては任意精度 <- なんじゃそりゃ
わかってて書いてるのかも知れないが、システムによって機械精度の桁が違うから。
計算速度のために機械精度で扱える時はあつかおうとするという
Mathematicaの設計の結果。

> リファレンス「入力の文法」にはdigits.digitsは「近似数」としか書かれてない。
digits.digitsって？

>>670
そのプログラムバグってるぞ。
x = 0.9999999999999999999999999
で試してみ。
# というかRationalize[x,0]の不完全版を回りくどくやっているだけ。
欲しいのは
FromDigits[RealDigits[x], 10]
あたりか？

>>671
状況が良くわからないが、
If[x > 0, ..., ...][[2]]
とか、もしかしたら
Simplify[If[x > 0, ..., ...], x > 0]
とか。
673 名前：１３２人目の素数さん mailto:sage [03/08/23 20:36]: >>668
わかった > digits.digits
機械精度にしろ任意精度にしろ近似数、ということ。
# 内部表現が2進なら必然的にそうなるな。
674 名前：671 mailto:sage [03/08/23 21:46]: >>672
Simplify[If[x > 0, ..., ...], x > 0]
でできました！
ありがとうございました
675 名前：１３２人目の素数さん mailto:sage [03/08/23 22:40]: >>672
あれー、なんぜバグっちゃうんだろー。
Rationalize[0.42857142857142857142857142857]とかの場合3/7になっちゃうから
そうならない版を作ったつもりだったんだけど…。
FromDigits[RealDigits[x], 10]かあ。すごいです。つーか、こんな関数しらなかった…。
676 名前：１３２人目の素数さん [03/08/24 16:33]: Mathematicaで平方完成する方法を教えてください
677 名前：１３２人目の素数さん [03/08/24 17:41]: 降ベキの順に表示させる方法ってないの？
678 名前：１３２人目の素数さん mailto:sage [03/08/24 20:18]: >>676
H[f_, x_] /; PolynomialQ[f, x] :=
Module[{r, b, c},
r = (x - b)^2 + c; r /. Flatten[SolveAlways[f == r, x]]]

Example: H[x^2 + 3 x + 2, x]

>>677
TraditionalFormか、別の方法としては、
HoldFormをもちいた並び替え関数を書く。
679 名前：１３２人目の素数さん mailto:??? [03/08/24 20:56]: 二次の係数は?
680 名前：１３２人目の素数さん mailto:sage [03/08/24 21:09]: >>679
それぐらいは自分で直せよ…。
681 名前：１３２人目の素数さん mailto:sage [03/08/24 21:11]: 日本語ﾍﾙﾌﾟ欲しーーーーーーーーーーー！。
いちいちmathematicaの公式サイトみるの大変…。
つーか、日本語版欲しい。
682 名前：１３２人目の素数さん mailto:sage [03/08/24 22:12]: >681
日本語キット(Help 日本語) が \50k ぐらいで売ってるよ。
というか、もう日本語ヘルプ強制バンドル品しか買えない。

Mathematica のヘルプの英語ぐらいは読めるように
しておいた方が何かと便利だと思うんですけど。。。
やっぱり日本語ヘルプのニーズは多いのかなぁ。
683 名前：１３２人目の素数さん mailto:sage [03/08/25 00:02]: >>682
高すぎ…。
684 名前：１３２人目の素数さん [03/08/25 00:36]: >>678
Mathematicaの関数としては、ないわけですね。自分で作るしかないってことか。

ありがとうございました。
685 名前：１３２人目の素数さん mailto:sage [03/08/25 01:37]: このソフトの中の人って頭良いけどちょっと融通が利かないというか…。
686 名前：１３２人目の素数さん [03/08/25 02:29]: >>685
人が入ってたのかよｗ
687 名前：中の人 mailto:??? [03/08/25 12:20]: ユーザにもっと勉強してもらいたいものだ
688 名前：１３２人目の素数さん mailto:sage [03/08/25 16:52]: >>687
sin[x]とかやった時に
「Sin[x]じゃないですか？」
とか言ってる暇があったらそれで計算しろカス
689 名前：１３２人目の素数さん mailto:sage [03/08/25 17:29]: >>688みたいなカス(例えば>>644)が多いなこのスレ。

実際はMathematicaの中の人は他の同系統のソフトの人と比べると、
「頭(知識)は普通(以下？)だが凄く融通が効く」
人だよ。
690 名前：１３２人目の素数さん mailto:sage [03/08/25 18:22]: >>689
人の平均から比べたら融通が利かない方だと思われ
691 名前：１３２人目の素数さん mailto:sage [03/08/25 18:46]: >>689は日本語を覚えたウルフラム
692 名前：１３２人目の素数さん mailto:sage [03/08/25 22:32]: まあ日本国内在住人にとっては腹の立つほど値段が高いこと、
それ以外は悪くないと思いますけどね。

もっともあの英語マニュアルを通して読めたのは
暇だった学生だったからなので、普通の生活では
ちと厳しいかな。

国産の RISA/ASIR も１０万くらい払っていいから
他の商用並の使い勝手にならんかな。
693 名前：１３２人目の素数さん mailto:sage [03/08/25 22:36]: >>690
それはそうだな。でもプログラミング言語ってのは厳密さがいいんだよ。
「sin?どうせSinと同じだろ」とか「+? ここは普通-だろ」とかやられたら
実際困る。
694 名前：１３２人目の素数さん [03/08/25 22:48]: >>692
要はMathematicaの買えない貧乏人のひがみだろｗ
695 名前：１３２人目の素数さん mailto:sage [03/08/25 23:23]: >>693
もちろんそんなのは分かってる。所謂ジョークだよ。
696 名前：１３２人目の素数さん mailto:??? [03/08/26 14:19]: list={{Aaa[a],Aaa[b],Aaa[c]},{Aaa[d],Aaa[e],Aaa[f]},{Aaa[g],Aaa[h],Aaa[i]}}
というリストがあって(ホントはもっと長い)、こいつを
{{a,b,c},{d,e,f},{g,h,i}}
にしたいのですが、エレガントな方法を教えてつかあさい。
Tableは使わないでね。
697 名前：１３２人目の素数さん mailto:sage [03/08/26 15:32]: >>696
(#[[1]] & /@ # &) /@ list
698 名前：697 mailto:sage [03/08/26 15:35]: >>696
Map[First, list, {2}]
のほうがいいな。
699 名前：１３２人目の素数さん mailto:??? [03/08/26 20:09]: >697
>698
できました。Thanks
さすがですねぇ
700 名前：１３２人目の素数さん [03/08/26 20:23]: 教えてつかあさい。
701 名前：１３２人目の素数さん mailto:sage [03/08/27 01:40]: mathematicaの初心者講座やって下さい。
702 名前：１３２人目の素数さん [03/08/27 02:54]: yada
703 名前：１３２人目の素数さん mailto:sage [03/08/27 07:25]: >>701
Shift+Enter
704 名前：１３２人目の素数さん mailto:??? [03/08/27 10:06]: MacではEnterだけで実行!
705 名前：supermathmania ◆ViEu89Okng [03/08/27 11:35]: Re:>701 MathematicaにはHelpが付いているはずだ。それを読むといいだろう。
まずは、Shift+EnterとExit[]と、+-*/^という演算子と、関数の引数は四角括弧で採ることを覚えよう。
706 名前：１３２人目の素数さん mailto:sage [03/08/27 19:21]: >>705
Exit[]知らなかった！！
707 名前：706 mailto:sage [03/08/27 19:22]: つーか、何の意味があるの？DOSみたいに終了するのかと思いきやしないし。
708 名前：１３２人目の素数さん mailto:??? [03/08/27 20:22]: Quit[]
とすると、カーネルが終了するのだ
わかったかね
709 名前：１３２人目の素数さん mailto:sage [03/08/27 20:41]: >>708
カーネルってなんですか？あとそれ以前にExitはどんなときに役に立つのですか？
710 名前：abc [03/08/27 21:37]: おかねもうけ情報見つけました

pocket.muvc.net/job/index.html
711 名前：１３２人目の素数さん [03/08/28 02:03]: Windowsです
[スタート]の[プログラム]メニューにあるMathematicaでカーネルを２まで同時に立ち上げることができる。
（プロフェッショナル版）
712 名前：１３２人目の素数さん [03/08/28 07:48]: xyではなくてx*yを使うの？
713 名前：１３２人目の素数さん mailto:??? [03/08/28 09:37]: カーネルってのは、
フライドチキンのおじさんのことだ
覚えておけよ
714 名前：１３２人目の素数さん mailto:sage [03/08/28 22:04]: 普通乗算はスペースだよな。
715 名前：１３２人目の素数さん mailto:sage [03/08/31 21:04]: Random[]で出る結果は何通りくらいある？
716 名前：１３２人目の素数さん mailto:sage [03/09/01 06:52]: >>715
意味わかんね
717 名前：１３２人目の素数さん mailto:sage [03/09/01 18:06]: >>716
Random[]で出る数は何種類から選ばれるかということです。
多ければ多いほど優秀、と。
718 名前：１３２人目の素数さん mailto:sage [03/09/01 23:07]: >715
Help に書いてないか？

まあ書いてないとして高々二桁 bit なんで十進数で
33 桁程度ではないか？基本的には剰余だろうから。。。
719 名前：１３２人目の素数さん mailto:sage [03/09/03 03:14]: >>718
ありがとう。
ピッタリ1とかピッタリ0ってありますか？
720 名前：１３２人目の素数さん mailto:sage [03/09/03 17:31]: SeedRandom[Random[]]ってやるとかなり増える？
721 名前：１３２人目の素数さん mailto:sage [03/09/03 22:21]: >719
ぴったりですか。。。原理的にはあるんだけど
実数乱数なら無いとしたほうがよく
整数乱数なら有りそうね。でもこれも Help に書いて無いか？

>720
基本的には Random[] は余りを周期的に利用していると考えて
そう外れがない。
SeedRandom は、この周期のどこから始めるかを指定する
だけなので、増えんよ。
722 名前：１３２人目の素数さん mailto:??? [03/09/05 15:45]: Mathematica Ver.5日本語版リリース! あげ
ttp://www.wolfram.com

でも代理店はどこもVer.4のままれす
723 名前：１３２人目の素数さん mailto:??? [03/09/11 18:05]: Ver.5
クラスルームってのは特別定価４万６千円だとよ
724 名前：１３２人目の素数さん mailto:sage [03/09/17 09:30]: mathematicaで計算実行したら、
”Internal counter overflow.
Some expression evaluations may not run to their
final fixed points. Try redoing your computation
on a 64-bit enhanced version of Mathematica.”
こんなエラー出てしまって計算やってくれないのですが、これってどういうことですか？
ちなみにVer.4です
725 名前：EXCITE 翻訳 [03/09/18 03:00]: 内部カウンタのオーバーフロー。いくつかの表現評価はそれらの最終固定
小数点で実行不可能かもしれません。Mathematicaの64ビット増強された
バージョン中のあなたの計算をやり直してみてください。
726 名前：１３２人目の素数さん mailto:sage [03/09/18 14:35]: k1=1/12; k2=1/6;

g[t_] = 2.30743*Exp[t^3.30352]*t^2.30352

h[t_] = (t*(k2+Sqrt[k2^2-t*k1^3]))^(1/3)

k[t_] = h[g[t]] + k1*g[t]/h[g[t]]

のとき、

NSolve[Integrate[h[t], {t, 0, tk}] == 1, {tk}]

を解きたいのですが、式が複雑すぎてエラーの嵐でした。
FindRootでもできませんでした。
先生によれば工夫すればできるらしいのですが、さっぱりわかりません・・・
どなたか教えていただけないでしょうか？

あと、g[t]の値はある計算結果なのですがそのままだと
2.30743*Exp[0.t^3.30352]*t^2.30352
となっているんですが、
0.t^3.30352　= t^3.30352
でいいんですよね？見たことない表現なので戸惑ったのですが・・・
なんでこんな表現になるんでしょうか？
727 名前：１３２人目の素数さん mailto:sage [03/09/18 22:18]: >>726
g も k も関係ない。

解けないのは式が複雑だからじゃなくて、使いかたが間違っているから。
NSolveはpolynomialを解くものであって、一般の式は解けない。

数値解は
FindRoot[NIntegrate[h[t], {t, 0, tk}] == 1, {tk, {1, 2}}]
等。(ただし、Sqrtの中身が負にならないtの範囲で)

0.t に見えたのは0. tで、0. * tだろう。その解釈は間違い。
728 名前：726 mailto:sage [03/09/18 23:20]: >>727
レスどうもです。
726にも書きましたがFindRootでもできませんでした・・・

また、0. * tって0*tとは違うんですよね・・・？
0. ＝桁数の小さい小数ってことですか？
729 名前：726 mailto:sage [03/09/18 23:31]: >>727

書き込んでくださった
FindRoot[NIntegrate[h[t], {t, 0, tk}] == 1, {tk, {1, 2}}]
でやってみたところできました。
FindRoot[Integrate[h[t], {t, 0, tk}] == 1, {tk, 1}]
でやったらできなかったのですが、なぜこれだとできなかったのでしょうか・・・？
730 名前：１３２人目の素数さん mailto:sage [03/09/18 23:36]: >>729
0.が常に0と区別できないRealsを表すのかどうかは知らない。
FullFormすればその0.の値はわかるけれど。

HelpのFindRootの項を読みなさい。
731 名前：726 mailto:sage [03/09/20 02:02]: >>730
FindRootについて勉強してみます
ありがとうございました
732 名前：１３２人目の素数さん mailto:sage [03/09/20 12:51]: >>731
ヘルプ読むのに10分もかからんだろ。そして読んだら、
> FindRoot[Integrate[h[t], {t, 0, tk}] == 1, {tk, 1}]
でもできないとおかしいと気づくはず。
答えが大体わかってたから、数値計算だけでもとめられる{tk,{1,2}}の
指定をしただけだ。

{tk,1}と指定すると、ニュートン法を行おうとする。
Integrate[h[t], {t,0,tk}]の微分は計算できる(h[t])から、
ニュートン法が適用できるはずだ。

実際遅いだけで、(Ver. 4.2では)解ける。
733 名前：１３２人目の素数さん [03/09/22 17:36]: どなたか教えて頂けないでしょうか。ほとほと困り果てております。
「1カラム→実数、　2カラム→複素数」
のテキストファイルをリストとして読み込ませるためには、どのよう
にしたらよいのでしょうか？
734 名前：supermathmania ◆ViEu89Okng [03/09/22 17:44]: Re:>733 C言語なら楽勝だよ。
735 名前：１３２人目の素数さん [03/09/22 17:59]: >733
何のリストに読み込ませたいんだい？
736 名前：１３２人目の素数さん [03/09/22 18:03]: >733
Excelに貼り付けて文字列関数でCONCATENATEを使って結合すれば？
括弧の部分とかは、固定の文字列として入れればいい。
737 名前：１３２人目の素数さん [03/09/22 18:43]: >>735
Mathematicaで読み込ませたいです。
738 名前：１３２人目の素数さん [03/09/22 18:45]: リストはこんな感じです。
　　　　　↓
実数　複素数
1　　3+I2
2 　　4+I5
3 　　8+I9
739 名前：１３２人目の素数さん mailto:sage [03/09/22 19:26]: >>739
data.datの中身が、
1　　3+2I
2 　　4+5I
3 　　8+9I
なら、
ReadList["data.dat", {Number, Expression}]

I2, I5, I9などとなっていたら、Symbolと理解されてしまうので、
あらかじめ加工しておくか、Stringで読み込んでMathematica内でパースする。
740 名前：１３２人目の素数さん [03/09/23 10:48]: >>739
本当にありがとうございます。これから試してみます！
741 名前：１３２人目の素数さん [03/09/23 17:38]: Readlist[#,Record]& で読み込み
Characters でばらす。
742 名前：１３２人目の素数さん mailto:sage [03/09/25 00:45]: 灯台に再びサイトライセンスで導入されるようだな。
1年の解析Bを取る連中は、課題を解く手間がずいぶん少なくなるだろう。
743 名前：１３２人目の素数さん mailto:sage [03/09/27 13:12]: 今入ってないのか

教官が授業で使っててずっけーって思ってた
744 名前：１３２人目の素数さん [03/10/05 21:23]: f[x_] :=
Print[b];
Return[x]

という関数を定義したとき
f[6]と入れても
bとプリントされるだけで、6が返ってきません。
どこが問題なのか分かる方教えてください。
745 名前：744 [03/10/05 21:37]: f[x_] :=
Module[{},
Print[b];
Return[x]
]

こう定義し直したら
f[6]を実行した時
bがプリントされて、６が返ってくるようになりました。

どうしてよくなったかは分からないですが…
746 名前：１３２人目の素数さん mailto:sage [03/10/05 22:00]: >>745
>>602を見ろ。
747 名前：744 [03/10/05 23:02]: >>746
ありがとうございます。
こう定義し直して、うまく行くことが分かりました。

f[x_] :=
(
Print[b];
Return[x]
)
748 名前：１３２人目の素数さん mailto:sage [03/10/07 23:53]: y^2=x^3-2 (-2<=x<=2)このグラフが書きたいのですが･･･
Plot[Sqrt[x^3 - 2], {x, -2, 2}]こう書くと上手くいきません。
どなたか教えてくださいませm(__)m
749 名前：１３２人目の素数さん [03/10/08 02:02]: >>748
そりゃあ、x^3-2<0の部分は書けないに決まってる。
750 名前：１３２人目の素数さん mailto:sage [03/10/08 15:25]: >>748
<<Graphics`ImplicitPlot`
ImplicitPlot[y^2 == x^3 - 2, {x, -2, 2}, AspectRatio -> 1]
751 名前：748 mailto:sage [03/10/08 21:53]: >>749さん
>>750さん
（ﾟ´Д｀ﾟ）どうもありがとうございます。助かりました。
752 名前：１３２人目の素数さん [03/10/12 18:49]: だれかWebMathematicaを使ってWeb上で数式処理が無料でできるサイトを構築してください！
俺にはMathematicaを個人的に買える予算などないんです！学生じゃないし！
753 名前：１３２人目の素数さん mailto:sage [03/10/13 02:23]: おとなしくみゅーぱっどでも使っとれ
754 名前：１３２人目の素数さん [03/10/15 23:37]: >>753

プロ版じゃないと見栄えが悪いと言ってみる。
755 名前：１３２人目の素数さん mailto:??? [03/10/16 09:39]: じゃ、MAXIMA使えばぁ
756 名前：１３２人目の素数さん mailto:sage [03/10/16 12:08]: フリーウェアでほとんどの数式をグラフにしてくれるソフトあるよね。
どんな極方程式もグラフにしてくれるやつ。すげー楽しい
757 名前：１３２人目の素数さん mailto:sage [03/10/19 14:52]: Functionのコマンドについて質問です
ｘ＾２＋ｙ＾２が０になるときのｘを求めるもので
条件としてｙはすでにリストになっています（例えば、y[[1]]=1,y[[2]]=2）
そこで、functionを使ってｘについての方程式をつくり求めようとしました
が、
A=Table[0,2];
Do[
A[[i]]=Function[{x,y},x^2+y^2][x,y];,{i,1,2}
];
として、A[[1]][1]と出力した場合にx^2+1[1]となってしまい、ｘに１を代入してくれません
理想としては１＋１（=2）となってほしいのです
どのようにすれば、求めることが出来るのでしょうか？
よろしくおねがいします
758 名前：１３２人目の素数さん mailto:sage [03/10/19 16:28]: >>757
貼っても動かない上に全く意味不明だ。
Table[0,2]とはなんだ。
(x^2 + y^2)[1]でなにがやりたいのだろうか。
759 名前：１３２人目の素数さん mailto:sage [03/10/19 16:38]: x^2 + y^2 が0になるような整数は(0,0)しかない。
760 名前：７５７ mailto:sage [03/10/19 18:15]: すいませんでした。実際ｘはもっと複雑な関数でして、わかりやすくするために
ｘにしてしまいました。
つまりやりたいことは、ある２つの変数（ｘ、ｙ）をもった式があります
条件としてｙはすでに決まっていてリスト化されています（例：ｙの１番目＝１、２番目＝２、ｎ番目・・・）
そこで、ｙがｎ番目の時に式＝０となるときのｘを求めたいのです
761 名前：１３２人目の素数さん mailto:sage [03/10/19 20:33]: y = Range[10];
A := Table[Function[x, x^2 + y[[i]]^2 // Evaluate], {i, 1, 10}];

A[[1]][1]

なんども書き込んでるんだけど反映されない…
762 名前：１３２人目の素数さん [03/10/19 21:17]: Plotを使ってグラフを描いたときに目盛りに書かれる数字を大きくする方法ってあるんでしょうか。
グラフを大きくしたいときに表示されている部分の角をマウスでドラッグして大きさを変えているんですが、
目盛りの数字が大きくならないんです。
最初に描かれるグラフの大きさも変えたいんですがこちらも分かりません。
全然数学な内容でなくて申し訳ありませんが教えてもらえないでしょうか。
763 名前：757 mailto:sage [03/10/19 22:29]: ありがとうございます。そのやり方でできました。
あの、そこで例えばx＾2+y^2=100と置いた時に
x^2+1=100, x^2+4=100 ・・・ x^2+100=100と１０個のリスト化された式がでます
そこで、xを求めるためにFindRootを使ったのですがリスト化されたものに対して
答えを出そうとするとうまく行きません
答えをリストとして出すためにはどのようにすればいいのでしょうか
よろしくおねがいします
764 名前：１３２人目の素数さん mailto:sage [03/10/19 23:09]: 少しは自分で考えろ。
765 名前：１３２人目の素数さん mailto:sage [03/10/19 23:25]: >>762
DefaultFont -> {"Bold", 10}
766 名前：１３２人目の素数さん mailto:sage [03/10/19 23:50]: >>765
ありがとうございます。明日早速試してみます。
>>764
すみません、Plot周りのオプションで探してはみていたのですが
「それじゃ見つかるわけがない」と猛反省しています。
767 名前：772 mailto:sage [03/10/20 14:05]: 自己解決いたしました。ありがとうございました
768 名前：１３２人目の素数さん [03/10/22 20:19]: mathematica5のStudent版買った人いる？
769 名前：１３２人目の素数さん [03/10/26 20:25]: 買ったよ。パスワード申請が相変わらず(｀×´)丿ウザイ!!
770 名前：１３２人目の素数さん [03/10/26 23:26]: 最近のMathematicaも
Windows版、Mac版、Linux版とかは一枚のCDで提供されているんですか？
少なくとも3.0の頃はNeXTとかも一緒に入っていて
プラットフォームの乗り換えの時、便利でしたが。
771 名前：１３２人目の素数さん mailto:age [03/10/27 00:45]: Winnyで落としたんだけど、なんかパスワードだか必要で使えん
意味ねー
772 名前：１３２人目の素数さん [03/10/27 01:08]: >>771
keyGeneratorがついてるだろ
773 名前：１３２人目の素数さん [03/10/27 01:23]: >>772
そんなアドバイスしてどうする…
774 名前：１３２人目の素数さん [03/10/27 02:52]: >>769
ネットでできるようになったんじゃなかったけ？
775 名前：１３２人目の素数さん [03/10/27 03:04]: ImportやExportのデフォルトのフォルダを指定する方法ありませんか？
初期設定は、Mathematica.exeと同じフォルダになってるけど
そこだと毎回面倒なので、デスクトップとかマイドキュメントに変えたいです。

今はImport["C:\Documents and Settings\galois\My Documents\sound.wav"]
とか書いてます。
776 名前：１３２人目の素数さん [03/10/27 08:51]: 今日アクティベーションが切れた。
来年から社会人だから糞高い値段かわんといけんのかと思うと鬱
777 名前：１３２人目の素数さん [03/10/28 17:17]: 会社で買ってもらえ（ｗ
778 名前：１３２人目の素数さん [03/10/28 20:59]: Yahooオークションで安いの出てるよ
でもたぶん違法だと思うが
779 名前：１３２人目の素数さん [03/10/30 00:57]: 高いねぇ、教員やっていても無茶苦茶高い。しかもバージョンアップが
毎年あって、その価格もとても高い。こんなのを学校で教育に使うという
ことは、ＭＡＴＨＥＭＡＴＩＣＡの会社に貢物をしているようなもので
正気じゃない。卒業生が馬鹿高いソフトを買うように仕込んで協力している
ようなものだよ。そんな金があるのなら、優等な学生に奨学金を出すべきだ。
780 名前：眠い mailto:sage [03/10/30 01:29]: 残念なことに、知ってる範囲のMathematicaユーザは優等生を兼ねてるね
優等生は金持の息子が多い、生活に余裕がないと自由に勉強する余裕がないからそうなる
781 名前：１３２人目の素数さん [03/10/30 20:19]: 貧乏人はフリーのMuPad,maximaを使えってことでつか
Mathematicaを持っているおいらはバージョンアップしないで
古いやつを大事にしまつ
782 名前：１３２人目の素数さん mailto:sage [03/10/30 22:06]: 忘れるな、MuPadは無料だが中身はブラックボックスだ。その点Mathematicaと一緒。

それと古いMathematicaのバグや劣ったアルゴリズムは懐で温めても直らんよ。
783 名前：１３２人目の素数さん [03/10/30 23:08]: 激しく既出だが、mupadはmathematicaよりmapleの方に近い。
ところで、mathematica4.0を持っているのだが、バージョンアップできるのは
5.0までか？それとも、5.1、5.2まで待てるか？
784 名前：１３２人目の素数さん [03/10/30 23:33]: 4.0のupは5.0までだとDMにあったような
785 名前：１３２人目の素数さん [03/10/31 23:29]: オープンソースの数式処理ソフトって何?
786 名前：１３２人目の素数さん [03/11/05 20:16]: 東西南北に行けるという2次元ランダムウォークをなんとか自力で組み立てた
んですが、格子のない（今いるところから半径１の円周上に動く）ランダムウォークが
どうしてもつくれないんですけど、どうしたら作れますか？アドバイス下さい。
当方まだmathematicaさわってまだ1ヶ月なもんでさっぱり知識ありません(^^;
自分で調べようにもいきなり学校でやれっていわれたもんなんで参考書すらないのです･･･
参考(?)までに格子のあるランダムウォークの肝心な所を載せます↓
Walk2D[n_] :=
FoldList[Plus, {0, 0}, {{1, 0}, {0, 1}, {0, -1}, {-1, 0}}
[[Table[Random[Integer, {1, 4}], {n}]]]]

いきなり長文＆質問すいませんです
787 名前：１３２人目の素数さん mailto:sage [03/11/05 21:33]: >>768
Walk2D[n_Integer] :=
FoldList[Plus, {0,0},
{Cos[#], Sin[#]}& /@ Table[Random[Real,{0,2 Pi}], {n}]]

ListPlot[Walk2D[60], PlotJoined->True]
788 名前：786 [03/11/07 01:22]: >>787
おぉ、すごい！！本当にランダムな角度で動いてますね！！！！
皆さんこんなの朝飯前ぐらいな問題ですか？

どーもありがとうございましたm(_ _)mﾍﾟｺ
789 名前：１３２人目の素数さん [03/11/07 21:57]: 教えてつかあさい。

ttp://phong.informatik.uni-leipzig.de/~kuska/mathgl3dv3/
ここにあるMathGL3d for MacOS X
というのを使いたいのですが、
インストールできません

MacOS XのMathematicaはパッケージなんで
サブフォルダに必要なファイルを入れることができないんです
??
790 名前：１３２人目の素数さん [03/11/08 00:04]: >789
拡張子
.app
を取ると、フォルダになります。

それでいいのかわかりませんが
791 名前：１３２人目の素数さん [03/11/08 13:26]: SolveやNSolveなどで複素数の答えがでてしまう式があります。

答えをCases[%, _Real]などで実数の範囲内を出すことができますが、
答えを編集するのではなく
もともとの式の変数の範囲を固定して
計算の高速化なり、入力する式達の簡略化はできるのでしょうか？

速度は不満ではないのですが、少し気になったので。
792 名前：１３２人目の素数さん [03/11/08 18:48]: >789
コントロールキーを押しながら
Mathematicaのアイコンをクリック
793 名前：１３２人目の素数さん [03/11/13 19:23]: Mathematica でバイナリファイルを読み込もうと

In[1]:= Needs["Utilities`BinaryFiles`"]
In[2]:= ReadListBinary["filename",Int16]

としているんですが処理が全然終わらないんです。
webで検索して
library.wolfram.com/infocenter/MathSource/354/binary.zip
をみつけたのですが、mccがうまく動いていないのか

% mcc binary.tm -o binary
"binary.tm", line 189.15: 1506-195 (S) Integral constant expression with a value greater than zero is required.
(このあともいっぱい似たようなのが表示されます)

となりビルドできません。
実は読もうとしているデータはエンディアンが違い、
読めた後も一工夫しなければいけなさそうでそちらを相談したかったんですが、
それ以前に読めなくて行き詰まってしまっています。
どなたかこの原因をご存知の方、対処法を教えていただけないでしょうか。
Mathematicaのバージョンは Mathematica 4.1 for IBM AIX です。
794 名前：大学生 [03/11/13 22:57]: 大学で課題がでたんですがみなさんの案で私を助けてください！

問題（例えば高校の参考書の問題や例題、大学の問題でも可）をあげてそれの解答の一部分をマスマティカを使って
視覚的に解説せよ。（アニメーションを使うと理想的）

ってのがでたんですが、イイ案ありませんか？人とはちがったのをやりたいんですが・・・
私には知識がありません。おねがいします。

アド載せとくんでマスマティカのソースファイルなんかを添付してくれれば一生感謝します。
どうかおねがいします！もちろん書き込みしてくれても感謝ですｍ（＿＿）ｍ
795 名前：大学生 mailto:mrhnr721@ybb.ne.jp [03/11/13 23:03]: >>794です。アド載ってませんでした。
796 名前：１３２人目の素数さん [03/11/14 02:54]: y=-x/(x^2+2)のグラフとその接点と接線を書きたいのですが、グラフと接点をtableとshowを使って
動かすのまではうまくいくのですが、接線が動きません、てか接線だけすら表示できません。教えてください。
797 名前：朝から呑んでる大学生 mailto:sage [03/11/14 09:41]: >>796
暇つぶしになんとなく
Do[{
bmp=Graphics[{RGBColor[1,0,0],PointSize[0.02],
Point[{x,x/(x^2+2)}],{RGBColor[0,1,0],
Line[Table[{k,k/(k^2+2)},{k,-10,10,.2}]]},{RGBColor[0,0,1],
Line[Table[{p,(2-x^2) p/(x^2+2)^2+2 x^3/(x^2+2)^2},{p,-10,
10,.2}]]}}];
Show[bmp]},{x,-10,10}]
798 名前：１３２人目の素数さん mailto:sage [03/11/14 10:10]: >>797
何呑んでるの？
799 名前：１３２人目の素数さん [03/11/14 19:45]: ｘ軸とｙ軸を描いてy=-x/(x^2+2)のグラフを固定し接線を動かすにはAxes->Trueとなにを使えばいいんですか？
800 名前：１３２人目の素数さん [03/11/14 20:22]: >>799
PlotRange
801 名前：797 [03/11/14 21:00]: >>798
マサイの戦士
802 名前：１３２人目の素数さん mailto:sage [03/11/14 23:21]: >>801
「呑んでる」言うから酒かと思ったらソフトドリンクかい
803 名前：１３２人目の素数さん [03/11/17 17:40]: >>802
マサイの戦士は「ソフト」ドリンクですか？？
804 名前：マティカの戦士 mailto:sage [03/11/18 11:05]: >>803
こんな飲み物
www.otsuka.co.jp/mss/
805 名前：狼 [03/11/18 22:27]: あのーサイクロイドをマスマティカで書きたいんですがどうすればいいですか？？
てかサイクロイドって円の一点を固定（？）して円を転がしてその一点の動いた軌跡じゃないですか？
私は円を転がしてその軌跡が描かれていく様子をアニメーションを使って表そうとしてるんですが、
うまくいきません。円と一点と軌跡をShowをつかってあらわしたいんです。
どうかだれか教えてください。（日本語下手ですいません）
806 名前：１３２人目の素数さん [03/11/18 22:52]: >>805
Mathematicaのヘルプにサンプルがあるだろ
807 名前：狼 [03/11/18 23:30]: help英語でわかりません。コピペしていただけたらありがたいです
808 名前：１３２人目の素数さん [03/11/18 23:59]: XpでMathematica5を使ってると，ときどきカーネルが落ちるんだけど，
同じような人いる？
809 名前：１３２人目の素数さん mailto:sage [03/11/19 00:16]: >>808
重い計算させてりゃＸＰでなくてもカーネルは落ちる
810 名前：１３２人目の素数さん mailto:sage [03/11/19 10:42]: >>808
カーネルの中の人も大変だからな
811 名前：supermathmania ◆ViEu89Okng [03/11/19 12:51]: いいこと教えてあげよう。
Mathematicaを使うときは、なるべく他のアプリケーションは切る。
これを守れば今日からジンセイバラ色(?)
812 名前：supermathmania ◆ViEu89Okng [03/11/19 12:54]: 私がいつかやってしまった事だが、
TeXのdviファイルのプレビューと、Mathmaticaを同時に動かしたとき、フリーズした。
気をつけよう。
813 名前：１３２人目の素数さん mailto:sage [03/11/19 14:41]: >>810
そういえば中の人などいないって聞かないね。
これも中の人がいるということが広く認知されたということか。
814 名前：狼 [03/11/19 15:41]: a
815 名前：１３２人目の素数さん [03/11/19 17:59]: 811 名前：supermathmania ◆ViEu89Okng ：03/11/19 12:51
いいこと教えてあげよう。
Mathematicaを使うときは、なるべく他のアプリケーションは切る。
これを守れば今日からジンセイバラ色(?)

　　　　　　　　　　　　　　ジンセイバラ色(?)

ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ
　　　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　ﾊﾟｼｬ　　ﾊﾟｼｬ　ﾊﾟｼｬ　
　∧_∧　　　　　∧_∧　　　　 ∧_∧　　∧_∧　　　　∧_∧　　　　 ∧_∧
　(　　 )】　　　(　　 )】　　　 (　　 )】【(　　 )　　　【(　　 )　　　【(　　 )
　/　 /┘　.　　/　 /┘.　　　/　 /┘ └＼＼　　　└＼＼　　　└＼＼
ノ￣ゝ　　　　ノ￣ゝ　　　　ノ￣ゝ　　　　ノ￣ゝ　　　ノ￣ゝ　　　ノ￣ゝ
816 名前：１３２人目の素数さん mailto:sage [03/11/19 19:48]: >>811
おまいの環境がわかるわｗ
817 名前：１３２人目の素数さん [03/11/22 02:47]: 行列を貼り付けたのですが
Dimensions[a]で実行したところ
Out[数字]＝｛100,100｝とならずに
｛1,100｝
Out[数字]＝｛100｝
と成ってしまいます。｛1,100｝ならわかるんですが、、

問題解決できるひといませんか？（これだけではわからないかな・・・）
818 名前：１３２人目の素数さん mailto:sage [03/11/22 11:08]: >817
行列ってのは
{ {a1 , a2 , ... } , { b1 , b2 , ... } , ... }

てな形なので、２つの階層のどっちかで
, コンマが抜けているのでは。。。
819 名前：１３２人目の素数さん [03/11/22 17:09]: supermathmaniaって実社会でも干されてんの？
820 名前：１３２人目の素数さん [03/11/22 23:01]: マセマティカでべき級数展開ってできますか？
821 名前：１３２人目の素数さん [03/11/22 23:05]: すいません、ちょっと質問なんですけど、
教科書で使われてるような図をMathematicaで作成したいのですけど、
ある領域を斜線で塗りつぶすような事ってできますか？
822 名前：１３２人目の素数さん [03/11/22 23:55]: >>821
液晶汚してもいいならできるよ
823 名前：１３２人目の素数さん mailto:sage [03/11/23 00:17]: >>821
色塗はできるが斜線は自分で引くしかないんじゃないか。
824 名前：821 mailto:sage [03/11/23 00:58]: そうですか、わかりました～ありがとうございます。
自分で作ってみますね。
825 名前：１３２人目の素数さん [03/11/23 22:21]: >>807
www.synapse.ne.jp/~dozono/math/anime/cycloid.htm
826 名前：１３２人目の素数さん [03/11/25 13:17]: >>818
ありがとうございます。違うところでの入力ミスでした。
でも、数百×数百の逆行列を計算しようとして、結果メモリー不足になり
計算できませんでした。普通のPCでこんな計算するのが間違ってましたｗ
827 名前：１３２人目の素数さん [03/12/04 03:09]: EmacsのMathematica-modeなんてあるのを知った。
フロントエンドよりマシな環境ができたらうれしいな～。

とりあえずコマンドの補完とか括弧の補完とか
それくらいやってほしいなぁ…
828 名前：＿ [03/12/05 03:31]: windows でmaximaをtexmacsで起動できた人っている？
両方インストールした(texmacsのみcygwin)けどtexmacsでmaximaを起動しようとすると,
bad url in "complete"
というerrorが表示され終了されます｡できた人がいれば教えてください
829 名前：１３２人目の素数さん [03/12/05 05:12]: MATHEMATICAって値段が高過ぎ。生協価格が20万円を超えている。
ふざけるなって言いたい。しかも毎年バ-ジョンアップが出てて、
アップデ-トに半額以上払わされる。なんだかぼられているのじゃ
ないのかな？アメリカだと10分の1の値段だったりしない？
830 名前：１３２人目の素数さん mailto:sage [03/12/05 05:51]: >>829
アカデミック価格だと、まあまあ買える値段だよ。
831 名前：１３２人目の素数さん mailto:sage [03/12/05 08:38]: >828
環境が違うのでなんともいえないが
locate コマンドの db がまだ更新されていないのでは？
僕はそうでした。

>829
1/10 ということはない。まあ実質半分かな
向こうでは $2000 を超えていないはず。
日本は、日本語環境 5 万円に１年保守費３０％
強制付きだからね。
もう個人は買うなってこと。。。
１０年前は海外から３万円でアップデートして
くれたのにな
832 名前：＿ [03/12/06 06:13]: >>831
レスthx。更新しても駄目でした｡
xmaximaは動くのですがmaximaをcygwinのxtermで動かす事ができなかったので,
その辺が原因だと考えてます｡誰か分かる人いたら教えてください｡
激しくすれ違いスマソ
833 名前：１３２人目の素数さん mailto:sage [03/12/06 23:45]: you should not ever use Japanese.
834 名前：１３２人目の素数さん [03/12/07 10:07]: {{a,b},{c,d},{e,f}}のようなリストを
{a,b,c,d,e,f}にする簡単な方法はありませんか？
835 名前：１３２人目の素数さん mailto:sage [03/12/07 10:43]: Flatten[%]でいけました。よごし失礼。
836 名前：１３２人目の素数さん [03/12/07 18:02]: >829
生協に商用版がおいてあるの?
837 名前：１３２人目の素数さん mailto:sage [03/12/07 23:53]: この関数を

fbbs = Flatten[Cases[
ReplaceAll[{u, t},
NSolve[{fbb[#1, #2, #3, #4], fbb[#5, #6, #7, #8]}, {u, t}]],
{_Real, _Real}]]&;

こんな↓感じにする、いい方法はないかな。
目的は可読性の向上。

NSolve[{fbb[#1, #2, #3, #4], fbb[#5, #6, #7, #8]}, {u, t}]
ReplaceAll[{u, t}, %]
Cases[%, {_Real, _Real}]
fbbst2 = Flatten[%]&

↑の方法だとうまくいかないのです。
Flatten[%]がそのままはいってしまいます。
838 名前：１３２人目の素数さん [03/12/08 16:31]: >>836
Student版が29000で売ってる
software.univcoop.or.jp/products/mathematica/student.html
839 名前：１３２人目の素数さん [03/12/10 20:33]: MATHEMATICA4.0をUNIXにインストールしようと思ったのですが、
パスワードだけ書き残していてライセンスIDをなくしてしまっていました。
さてここで問題です。
いったい私はどうしたら良いのでしょうか？
840 名前：１３２人目の素数さん [03/12/10 20:43]: ノーヒントですか？
相談なしですか？
841 名前：１３２人目の素数さん mailto:sage [03/12/10 22:05]: >>840
俺のライセンスＩＤなら・・・
842 名前：１３２人目の素数さん [03/12/11 11:35]: ライセンスIDってCD-ROMの袋に書いてなかったっけ？
ライセンスIDって無くしたら大変な気がするけど。。。
843 名前：１３２人目の素数さん [03/12/11 12:27]: 誰かただでください！
844 名前：１３２人目の素数さん mailto:sage [03/12/11 15:33]: C言語から評価式と、入力したパラメータをMathematicaに送り，Mathematicaで計算，
結果をCプログラムに返して表示される、というものを作ろうとしています．
で，関数を書いてみたのですが，いくつかわからない点がありましたので質問させてください．

①MathKernel.exeのパスの指定はどのようにするのでしょうか？
CプログラムのほうからMathKernel.exeのパスを指定し，char変数PATHに格納してあるとすると，
link = MLOpenString(env, "-linklunch -linkname 'PATH'", &err);
でいいのでしょうか？

②パッケージの読み込み時に使うのは以下のようにMLPutStringで良いのでしょうか？
MLPutFunction(link, "Get", 1);
MLPutString(link, "Package");

③一時変数はC言語側で定義するのでしょうか？
たとえば，以下の評価式で変数y,zはCプログラムのほうから入力するとき，y,zはC言語側で定義しますが，一時変数xもC言語側で定義するのでしょうか？(用いる変数はMathmatica内部で使用するのみであってもすべてCで定義するのでしょうか？)
temp = y + z
x = temp^temp
845 名前：844 mailto:sage [03/12/11 15:34]: ④定義関数f[t_]=～を送信するのはどうすれば良いのでしょうか？
f[Pattern[t, Blank[]]から，以下のようにしてみたのですが・・・

MLPutFunction(link, f, 1); //←これでいいのか？
MLPutFunction(link, "Pattern", 2);
MLPutDouble(link, t);
MLPutFunction(link, "Blank");

⑤(x^y + z)のように関数Plusの引数に関数Powerが入っている場合，送信するのは
MLPutFanction(link, "Plus", 2);
MLPutFanction(link, "Power", 2); //←ちゃんとPlusの引数と見なされるのか？
MLPutInteger(link, x);
MLPutInteger(link, y);
MLPutInteger(link, z);
でいいのでしょうか？

長くなりましたが，よろしくお願いします m(_ _)m
846 名前：844 mailto:sage [03/12/11 15:36]: ③に関して・・・
一時変数xも→一時変数tempと結果を代入する変数xも
の間違いです
847 名前：こてはん [03/12/11 19:35]: >>842
それほんと？
探してみたら俺のも行方不明でした・・・。
ピンチ？
848 名前：１３２人目の素数さん mailto:sage [03/12/11 22:06]: >>843
ヤフオクに出てる違法コピー買って摘発されてくださいｗ
849 名前：１３２人目の素数さん mailto:sage [03/12/11 22:29]: 今M2で今年で卒業なんだけど
卒業までにstudent版を買ってら、来年も問題なく使えますか？
それとも自動的にライセンスが切れる仕組みになってるの？
卒業までに残りの貯金でソフトを購入する覚悟はできたのですが。
850 名前：837 mailto:sage [03/12/12 00:30]: >>844

E:\Program Files\Wolfram Research\Mathematica\4.0\AddOns\
MathLink\DevelopersKits\Windows\UnsupportedGoodies\VCFE\Source\
と
ヘルプのアドオンMathLink当たりにすべての答えがあると思います。

自作のprogramはwrapper classを作って
頭にMLの文字とMLINKの引数が無いのでそれを前提として読んで見てください。

１、これを渡していました。

"-linkname \'\"E:\\Program Files\\Wolfram Research\\Mathematica\\
4.0\\MathKernel.exe\" -mathlink\' -linkmode Launch";
err checkはしていませんでした。

> link = MLOpenString(env, "-linklunch -linkname 'PATH'", &err);
ちゅうかsprintfしないと……。

２、たぶん、EndPacket();が必要です。
851 名前：837 mailto:sage [03/12/12 00:30]: ３、y + zとtemp ^ temp を複数のPutFunctionでやる場合は結果を受け取る時に必要です。
それを一つの関数にしてしまいひとつのPutFunctionでやる場合は必要ありません。
F.E.
y + z
temp^temp
それぞれ送信する場合、それぞれ答えが帰ってくる。
x = (y + z) ^ (y + z)
これなら一つの結果しか返ってこない。
temp = y + z
x = temp^temp
これを一気に送信したら、２つの答えが帰ってくる。
temp = y + z;
x = temp^temp
たぶんこうして一気に送信したら一つだけだと思う。

４、初期化時に
char *型の配列を用意して関数を登録(VCFEのサンプルにあるEvaluateを呼び出し)
してそれをPutFunctionしている。
つまり
fx = f[Pattern[#1, Blank[]]&;
を登録して
PutFunction(fx, 1);
PutDouble(t);
EndPacket();
852 名前：837 mailto:sage [03/12/12 00:31]: ５、それはよくわからないからtry&errorで。
僕の場合は４の通り、使う関数をすべてあらかじめ登録してから使っています。

おまけ、
EndPacket後に、答え以外のパケットが帰ってくる場合の処理が面倒です。
それは試行錯誤してみてください。

たぶん例に出したような簡単な処理をするのではないと思いますが、
もし桁数を考えないで計算したいのであれば、多倍長ライブラリで検索して見てください。
GMP,NTLがおすすめ。

むむむ、質問に答える事で、自分の質問の答えがわかった。thank you。
853 名前：１３２人目の素数さん [03/12/12 12:03]: >>847
ホントだよ。
今使えてるならば

C:\Program Files\Wolfram Research\Mathematica\4.2\Configuration\Licensing\mathpass
ってファイルにライセンス情報が書かれていて
そこにちゃんとライセンス番号もあるから大丈夫
今の内に控えておこう。

ちなみにこのmathpassファイルはOSを再インストールしなかったら使い回しできる。
いちいちWolframにパスワードを聞かなくて済むので便利。
OSを入れ替えたり、別のパソコンにインストールしようとした場合は
MathIDが変わるからWolframのお世話にならないとダメ。

ライセンスIDが完全に分からない場合はどうするのかねえ。
Wolframに聞いてみるとかか？
>>839はパスワードはあるみたいだし、是非聞いて欲しい。
報告待ってま～す。
854 名前：１３２人目の素数さん mailto:sage [03/12/12 12:09]: >>849
問題なく使えるけどライセンス違反だぞ。
次回パスワードをもらうときに、困ったりするかもしれなかったりするかもしれない。。。

たぶん困らないけど。。。

昔はWolframに英語のメール出したりしてパスワード再発行しないといけなくて大変だった。
英語苦手だから、訳分からない文章だったろうなぁ。。。
Hey!とかから始めたし。
855 名前：１３２人目の素数さん mailto:sage [03/12/12 12:30]: >>849
俺もstudent版買ってパスワードもらう時、卒業時期を
思いっきりごまかしといたからだいぶ使えそう。
アップグレードもしないと思うし。
856 名前：844 mailto:sage [03/12/12 18:45]: >>850
とても詳しいお答えありがとうございました．
提示していただいたサンプルプログラムやヘルプも参考にしてみたのですが，下記の点についていまいちわかりませんでした．

①
>> link = MLOpenString(env, "-linklunch -linkname 'PATH'", &err);
>ちゅうかsprintfしないと……。

wsprintf(launchstring, " -linklunch -linkname '%s'", PATH);
link = MLOpenString(env, launchstring, &err);
としてみたのですが，提示していただいたサンプルプログラムを見てみると，LinkProtocolの指定を行っていますが，これは必要あるのでしょうか？

②
>たぶん、EndPacket();が必要です
パッケージ読み込みの送信をしてMLEndPacket(link);を送信，その後に評価式の送信を行う場合はMLNextPacket(link);を送信する必要があるのでしょうか？

③
>y + zとtemp ^ temp を複数のPutFunctionでやる場合は結果を受け取る時に必要です。
>それを一つの関数にしてしまいひとつのPutFunctionでやる場合は必要ありません。

y + z の結果tempの値をCプログラムの方に返さない、つまり
temp = y + z;
を送ってやれば，変数定義する必要はないということでしょうか？
またtemp = y + z;の送信をする際，
CompoundExpression[Set[temp, Plus[y, z]], NULL]
の引数NULLを送信するときの関数は何を用いるのでしょうか？
857 名前：844 mailto:sage [03/12/12 18:47]: ④
大変恐縮ではありますが，関数の登録について詳しく教えていただけないでしょうか？
VCFEのEvaluate関数のところを見てみたのですが，正直なところほとんどわかりませんでした・・・
もちろん自分で理解しようと勉強してはいるのですが，苦戦している状況です・・・

⑤
VCFEをビルドしてみたところ
LINK : fatal error LNK1104: ﾌｧｲﾙ "..\..\..\..\Program Files\Microsoft Visual

Studio\VC98\Lib\ml32i2m.lib" を開けません。
とエラーが出てビルドできません．ml32i2m.libは
C:\Program Files\Microsoft Visual Studio\VC98\Lib
にコピーしたのですが，できませんでした．他に何が必要なのでしょうか？

当方MathematicaもC言語も学び始めて間もないので理解が浅く，お手数とは存じますがよろしくお願いしますm(_ _)m
858 名前：839 [03/12/12 21:01]: >>853
むぅ～、mathpassに身に覚えの無い数字とアルファベットが
タブ区切りで入っているのですがこれをどう理解すればよいの？
859 名前：837 mailto:sage [03/12/13 00:27]: >>844,856-857
＊LNK1104を調べて見て。
ようはlibはcopyをするだけで使えるようになるわけではなく、
リンカにこのlibを使うと教えてあげなければ、使えるようにはなりません。
VisualStadioのプロジェクトのプロパティでリンカの部分で指定してあげてください。

また、libはcopyするのではなく、libがあるフォルダのパスを通してあげてください。
menu,ツール->オプション,プロジェクトでライブラリファイルが指定できるはずです。
（VS.Net 2003であれば）

＊僕もMathLinkをバリバリ使っているというわけではないので、間違えが含まれるかも知れません。それを前提として読んでください。
僕の前の発言は少し間違いがふくまれているかもしれません。
こちらでEvaluate関数を改造してしまっていたので、そのまま式を渡しても自作関数を
登録（定義）したことにはならないかもしれませんし、いけるかもしれません。try&errorで。

＊MathLinkは二種類の式の呼び出し方をサポートしています。

ひとつは関数を指定して、値を一つずつ渡す方法。
PutFunction("Plus", 2);
PutReal(hoge);
PutReal(hage);
EndPacket();

もうひとつが、式の文字列をそのまま渡す方法。（TextPacketかEnterTextPacket
（違いはみていない）の関数の引数として式の文字列を渡す。）

PutFunction("TextPacke"/*不安*/, 1);
PutByteString((const BYTE *)"hoge + hage");
EndPacket();
860 名前：837 mailto:sage [03/12/13 00:27]: 上の例だとGetReal();で答えを取れるはずですが、
下の例だとhelpファイルを見れば分かりますが、答えがそのまま文字列として
帰ってくるので、GetString(&OutputString);でとればOKです。

VCFEのEvaluate関数は下の例でやっています。

＊関数の定義は別に難しい事をやっているわけではなく、
Mathematica上で関数の定義をするのとまったく一緒です。

// ベジェ曲線を媒介変数ｔを指定しxyposを取得
// #1 t #2 q0 #3 q1 #4 q2 #5 q3
char szInit_fb[] = "fb = (1 - #1)^3 #2 + 3(1 - #1)^2 #1 #3 + 3(1 - #1)#1^2 #4 + #1^3 #5&;";
// #1 t #2 q0x #3 q0y #4 q1x #5 q1y #6 q2x #7 q2y #8 q3x #9 q3y
char szInit_fb2[] = "fb2 = List[fb[#1, #2, #4, #6, #8], fb[#1, #3, #5, #7, #9]]&;";

この二つの式をMathmatica上で定義しておけば、
fb2[(* 略 *)]
と打ち込むだけで答えが返ってきます。
同じように自作プログラムでは、あらかじめEvaluate関数で上の二つを定義しておき、使う時に
PutFunction("fb2", 9);
/* 略 */
EndPacket();
で答えが返ってきます。

あとNullは固定シンボルであるようです。helpファイルを。
861 名前：１３２人目の素数さん [03/12/13 09:30]: >>858
どう理解って…
そんなの他人に入れてもらったとしか理解できんよ。。。
割れてるんじゃないの？
862 名前：858 [03/12/13 09:43]: >>861
自分で入れたんだけどね、
あそこにはパスなどの情報がそっくりそのままの状態で入ってるの？
識別できたのは自分で任意にいれた名前と所属ぐらいだったよ。
863 名前：１３２人目の素数さん [03/12/13 15:52]: >>862
そうだよ。
ライセンスIDとかMathIDとかパスワードがそのまま入ってる。
だから同じパソコンに同じOSをインストールして、
そこに同じMathematicaを入れたなら
このmathpassファイルはそのまま使える。

Wolframにパスワードを聞いても、全く同じものが返ってくるしね。
これは自分で確かめた。
864 名前：厨房君 [03/12/13 18:02]: わしは某大学の研究室でwindowsとlinux両方にインストールしているのだが
linuxのほうはmathpassにそのままデータが入っていたが
確かにwindowsのほうはデータがそのまま入っているようには見えないなぁ
ちなみに4.0でーす。
865 名前：１３２人目の素数さん mailto:sage [03/12/13 20:33]: >>864
そう？
俺はWindows2000でMathematica4.2だけど
全部そのまま入ってるよ。

まずパソコン名
その次がMathID
そしてライセンスID
パスワードで
最後がユーザー名

そのまま入ってるよ。
866 名前：858=862 [03/12/13 23:47]: >>865さん
もし問題が無ければ自分の手元にあるパスとmathpassに入っている
（タブで）4つめにあるデータを書き込んで一体どういうことなのか、
更に（タブで）3つめのデータからライセンスIDを推測するには
どうしたらよいかご教授いただきたいのですが、
問題ありですかねぇ。
867 名前：１３２人目の素数さん [03/12/14 03:13]: >>866
パスワードの方は公開しても何も問題ないと思う。
個々のMathematicaごとに違うから、バレても意味無いしね。

ライセンスIDの方はWolframとやりとりする時に必要そうだから
公開するのはやめた方がいいんじゃないかな？
そのmathpassに書いてある、訳分からない文字列の４番目の方も
そこからライセンスIDが分かったら困るしやめた方がいいかもね。
868 名前：１３２人目の素数さん [03/12/14 17:02]: ですよねぇ。とほほ・・・
869 名前：１３２人目の素数さん [03/12/14 17:43]: >>868
パスワードの方だけでも見せてくれない？
そのままじゃないってのがどういう状態なのか見てみたいです。
４番目の所だけ伏せて見せてくれたらベスト。
870 名前：868 [03/12/14 22:07]: >>868さん
一応部分的に*（中身は数字です）でふせておきます。
mathpassに入っていたデータをタブの順番通りにコピペすると以下の通りです。

aaaa <-自分で任意に入力しました。
**90F64AED623BCB72B3168A
F76EA93EB4840346902DCBE4A72B8F47E89417B223**
28495937**
aaaa <-自分で任意に入力しました。
aaaa <-自分で任意に入力しました。

です。
ちなみに手元のパスワードは
「49874-10775-47982-49844-53548-***」
です。
>>865の順番通りだとすれば
どう理解すればいいのやら・・・。
871 名前：１３２人目の素数さん mailto:sage [03/12/15 00:08]: 自分で任意に入れたって部分がよく分からないけど
（ユーザー名とパソコン名？じゃもう一つはなんだろ…）
それじゃあライセンスIDは分かりそうもないね。。。
872 名前：844 mailto:sage [03/12/15 22:21]: >>859,860さん
丁寧にありがとうございました。
教えてくださったこととVCFEを参考にしてソースを書いてみたいと思います
またなにかありましたらよろしくお願いしますm(_ _)m
873 名前：大学生 [03/12/17 20:22]: ParametricPlot3Dでメビウスの輪の描き方教えてください。
よろしくお願いします。
874 名前：１３２人目の素数さん mailto:sage [03/12/17 22:16]: ヘルプに載ってるよ
875 名前：８７３ [03/12/17 23:33]: ヘルプ見たんですがトーラスしかのってませんでした。どうか教えてくださいませんか？
お願いします。
876 名前：837 mailto:sage [03/12/18 16:20]: mathworld.wolfram.com/MoebiusStrip.html
877 名前：837 [03/12/19 02:31]: これすげぇ
メビウスの歯車だ。
mathworld.wolfram.com/gifs/moebgear.gif
878 名前：１３２人目の素数さん mailto:sage [03/12/19 03:13]: notebookも落とせるところがよいね
879 名前：844 mailto:sage [03/12/20 17:47]: CプログラムのほうからMathKernel.exeを立ち上げることができましたので，簡単な例をやってみようと思い，
和を計算させるソースを書いてみました.以下はその一部です

char expr[] = "2^4 + 3^2";
const char *OutputString;

case IDC_EVALUATE:
MLPutFunction(link, "EvaluatePacket", 1);
MLPutFunction(link, "ToExpression", 1);
MLPutString(link, expr);
MLEndPacket(link);
MLFlush(link);

while (MLNextPacket(link) != RETURNPKT)
MLNewPacket(link);

MLGetString(link, &OutputString);
wsprintf(szBuf, "OutputString = %s", OutputString);
MessageBox(hDlg, szBuf, "Debug", MB_OK);
MLDisownString(link, OutputString);

MLNewPacket(link);

return TRUE;

この例だとうまく行くのですが，
結果を取得する部分を，この場合は答えは25であり明らかに整数ですので，
MLGetIntegerを用いてみるとアプリケーションエラーで終了してしまいました.
MLGetIntegerなどはどのような場合に使用するのでしょうか？
880 名前：844 mailto:sage [03/12/20 17:48]: また，ほかの評価式の場合，例えば
char expr[] = " Solve[{x + y == 5, 3x + 5y ==13}, {x,y}]";
といった連立方程式の評価を試みてみたところ，同様にアプリケーションエラーが出て終了してしまいました.
結果の取得に用いる関数がMLGetStringでは正しくないためと思われるのですが，どの関数を用いて結果の取得を行えばいいのでしょうか？
MLGetRealList等試してみたのですが，うまくいきませんでした.

よろしくお願いします
881 名前：837 mailto:sage [03/12/21 17:07]: EvaluatePacket & ToExpression でもいけるんですね。
違いはわかりませんが、TextPacke、EnterTextPacketと同じく
式を文字列でおくり、結果を文字列で返す関数だと思います。

結果をGetIntegerでもらう場合はこうです。
PutFunction("Plus", 2);
PutInteger(1);
PutInteger(1);
EndPacket(link);

下の例は、GetStringでいいはずです。
VCFEで投げるとうまく行く＆エラーメッセージが
返ってくるような式ではないので
EvaluatePacket & ToExpressionがまずいか、
EndPacketから、GetStringの間の処理がまずいかどちらかです。
882 名前：837 mailto:sage [03/12/21 17:12]: ToExpressionですね。

In[373]:=
EvaluatePacket[ToExpression[{x + y == 5, 3x + 5y == 13}]]

\:6B21\:3088\:308A\:FF1AIn[373]:=
ToExpression::"notstrbox": "\!$\(\(x + y$\) == \
5\)h文字列あるいはボックスではありません.ToExpressionは文字列かボ\
ックスだけをMathematica入力として解釈することができます．"

\:6B21\:3088\:308A\:FF1AIn[373]:=
ToExpression::"notstrbox": "\!$\(\(\(\(3\\ x$\) + $\(5\\ y$\)\)\) == \
13\)h文字列あるいはボックスではありません.ToExpressionは文字列かボ\
ックスだけをMathematica入力として解釈することができます．"

Out[373]=
EvaluatePacket[{$Failed, $Failed}]

In[374]:=
EvaluatePacket[ToExpression[2^4 + 3^2]]

Out[374]=
EvaluatePacket[25]
883 名前：１３２人目の素数さん [03/12/22 10:54]: a,b,tをパラメータにしてグラフにしたいんですが,うまくいきません。
ｂのところを数字にすればできました。どうしたらいいでしょう？
sol = Table[NDSolve[{x''[t] == aCos[((b - a)) t^2] - bSin[((b - a))t^2],
　　　　　　　　　　　ｙ''[t] == aSin[((b - a) t^2] + bCos [((b - a)) t^2],
　　　　　　　　　x[0] == 0, y[0] == 0, x'[0] == 0, y'[0] == 0},
　　　　　　　　　{x, y}, {t, 0, 1}],
　　　　　　{a, (-1), 1, 0.2}, {b, (-1), 1, 0.2}];
ParametricPlot[Evaluate[{x[t], y[t]} /. sol], {t, 0, 1},
PlotRange -> {{0, 3}, {0, 3}}]です。
884 名前：１３２人目の素数さん [03/12/22 12:16]: Mathematica で、整数や多項式の行列の単因子(不変因子）を求めるには
どうしたらいいんでしょうか？
Maple だと、smith や ismith という関数があるんですが・・・
885 名前：１３２人目の素数さん mailto:sage [03/12/24 00:40]: >>883
このスレの一番最初にある関数は違うかな？
886 名前：１３２人目の素数さん [03/12/26 11:03]: FrameLabelのフォントサイズを指定するにはどうすれば良いのでしょうか？
FrameStyle->{FontSize->36}とか設定してもFrameLabelのフォントサイズ
は変わりません。
887 名前：844 mailto:sage [03/12/27 15:22]: >>837さん
遅レスになりましたが、レスありがとうございました
参考にさせていただき、がんばってみます
888 名前：１３２人目の素数さん [04/01/04 12:08]: age
889 名前：844 mailto:sage [04/01/07 20:43]: 879について再度教えていただけないでしょうか？

EvaluatePacket[ToExpression["Solve[{x + y == 5, 3x + 5y ==13}, {x,y}]"]]
としてMathematicaのほうで計算するとできたのですが，
char expr[] = "\"Solve[{x + y == 5, 3x + 5y ==13}, {x,y}]\"";
MLPutFunction(link, "EvaluatePacket", 1);
MLPutFunction(link, "ToExpression", 1);
MLPutString(link, expr);
とCのほうから送信してやってMLGetStringで受信しても，答えのほうにSolve[{x + y == 5, 3x + 5y ==13}, {x,y}]が出てきて不正な処理落ちしてしまいます・・・

また、VCFEのように，
char expr[] = "Solve[{x + y == 5, 3x + 5y ==13}, {x,y}]";
MLPutFunction(link, "EnterTextPacket", 1);
MLPutByteString(link, (const unsigned char*)expr, strlen(expr));
MLEndPacket(link);
で送信
MLGetString(link, &OutputString);
で受信というようにしてみたのですが，こちらもうまく行かず処理落ちしてしまいました・・・

どうすればよいのでしょうか？正直お手上げ状態です・・・
890 名前：837 mailto:sage [04/01/09 11:33]: VCFEと同じようにして落ちるって言うのはありえないから、
同じように出来ていないって事では。
Debugger立ち上げて一行ずつどう処理しているかみてみるといいよ。

たぶん、必要なFlash関数の呼び出しとかをしていないんだと思うよ。
891 名前：１３２人目の素数さん [04/01/09 14:37]: どなたか
Fourierseriesexpansion[1/(Pi+Cos[x])]
をしてください。
892 名前：１３２人目の素数さん mailto:sage [04/01/09 23:37]: In[1]:=
Fourierseriesexpansion[1/(Pi + Cos[x])]

Out[1]=
\!$Fourierseriesexpansion[1\/\(π + Cos[x]$]\)

ウェブ全体から検索日本語のページを検索
ウェブイメージグループディレクトリ

Fourierseriesexpansionに該当するページが見つかりませんでした。
"fourierseriesexpansion"を含むページは見つかりませんでした｡
893 名前：１３２人目の素数さん [04/01/10 02:37]: 初心者ですが、質問させてください。
Version5で新機能となったMaximizeって機能を使いたいのですが
私のVersionが古くて使えません。

そこで関数に領域を与え最大値を求める関数を実装したいと思いますが
どのようにしたら良いのでしょうか。
教えてください、お願いします。
894 名前：837 mailto:age [04/01/11 20:14]: ベジェ曲線の長さを求める式の不定積分を無理やり計算しようとしているのですが、
僕のマシンが非力なのか計算できません。
どなたかできませんか？
もし計算を途中で止める場合はAlt+.かAlt+,です。

(* www.interq.or.jp/www1/anzawa/doc/doc3.htm *)
(* #1 == t; #2,#3,#4,#5 == 1234 ;*)
fb = (1 - #1)^3 #2 + 3(1 - #1)^2 #1 #3 + 3(1 - #1)#1^2 #4 + #1^3 #5 &;

(* 媒介変数表示で表されたた曲線の長さを求める式
NIntegrateにして、x1などに値を入れ、２番目の引数を{t, 0, 1}にすれば求められる *)
Integrate[Sqrt[(fb'[t, x1, x2, x3, x4])^2 + (fb'[t, y1, y2, y3, y4])^2], t];
895 名前：837 mailto:age [04/01/11 20:22]: >> 893
やったことないけど
FindRoot, FindMinimum, ConstrainedMin, ConstrainedMax,
あたりかな？
素直にヴァージョン上げるのが一番早いと思うよ。
896 名前：844 mailto:sage [04/01/12 23:20]: >>837さん
レスありがとうございます
おっしゃるとおりFlash関数呼び出していませんでしたので，やってみます

グラフィック描画はどのようにすればよいのでしょうか？
(VCFEで簡単なプロットをさせてみようとしてもエラーが出て描画できないのですが，実装されていないのでしょうか？)

VCFEのソースを見たのですが、いまいちよくわかりません
WMFファイルのオープン？とかさっぱりです・・・
プロット画像情報のWMFファイルが作成されていて，それをOpenWMFFileで読み込み，DisplayGraphicで表示させていると解釈していいのでしょうか？
とすればどこにWMFファイルが作成されているのか等々問題点は尽きないのですが・・・

また，DISPLAYENDPKTが来ると描画させているようですが，DISPLAYPKTのときはどのような処理を行っているのでしょうか

？グラフィックウインドウに空文字列または出力文字列を表示させようとしているみたいですが・・・

質問ばかりで申し訳ありませんがよろしくお願いします
897 名前：837 mailto:sage [04/01/13 02:27]: >>894を無理やり計算して見ました。結果は押して知るべし。

FullSimplify[Sqrt[(fb'[t, x1, x2, x3, x4])^2 + (fb'[t, y1, y2, y3, y4])^2]]
(* 出てきた答えのルートの中の式をhogeとしてベキ級数展開 *)
Series[hoge, {t, 0, 4}]
(* そんでこんな形になると仮定。本当は係数の英文字の順が逆がいいな *)
f[t_] := 3 Sqrt[AA + BB t + CC t^2 DD t^3 + EE t^4]
(* どうぞ *)
Integrate[f[t], {t, 0, 1}]

>>896
> グラフィック描画はどのようにすればよいのでしょうか？
> (VCFEで簡単なプロットをさせてみようとしてもエラーが出て描画できないのですが，実装されていないのでしょうか？)
そう？
一番最初に設定されている式は結構時間かかるから待って見てごらん。
Plot3Dのサンプルにある下の式とか問題なく表示できたよ。
Plot3D[Sin[x y], {x, 0, 4}, {y, 0, 4}, PlotPoints -> 40, Mesh -> False,
FaceGrids -> All, AxesLabel -> {"Length", "Width", "Height"}];
Plotのサンプルは続けてやらないとダメダヨ。
t = {Sin[x], Sin[2x], Sin[3x]};
Plot[t, {x, 0, 2π}];
pp = Plot[Evaluate[t], {x, 0, 2π}];

> WMFファイルのオープン？とかさっぱりです・・・
たぶんWMFの知識とかソースが何をやっているか理解できないんでしょう（＾＾；
ソースに書いてある事をすべて説明する余裕はないので、
C言語の本などを読んで見てください。
言語の説明は板違いになるので、プログラム技術板でお待ちしています。
898 名前：837 mailto:sage [04/01/13 02:39]: ここでも出来た！
3 Sqrt[AA + BB x + CC x^2 + DD x^3 + EE x^4]
integrals.wolfram.com/
でも横が切れちゃってる……。
899 名前：844 mailto:sage [04/01/13 15:37]: >>837さん
デフォで入力してあるPlot3Dを実行してみると
OpenTemporary::noopen:
Cannot open C:\\DOCUME~1\\\214\212\340V\215G\216\:8180\\LOCALS~1\\Temp\\000001a01440.
Part::partd: Part specification $Failed \:301A 1 \:301B
is longer than depth of object.
Part::partd: Part specification $Failed \:301A 1 \:301B
is longer than depth of object.
General::stream:
$Failed is not a string, InputStream[ ], or OutputStream[ ].
Display::fname: Filename $Failed [[ 1 ]] is not a string.

とエラーが出てしまうんですよね・・・
多分WMFファイルが開けないって言ってると思うんですが・・・
(ソース等は何もいじってません)
900 名前：844 mailto:sage [04/01/13 16:13]: tempフォルダのパスに2バイト文字が入っていたのがいけなかったようで，
環境変数からパスを変更してみたところ表示できるようになりました

お騒がせしました・・・
901 名前：837 mailto:sage [04/01/13 20:19]: (,,ﾟДﾟ) ｶﾞﾝｶﾞﾚ!
902 名前：１３２人目の素数さん mailto:sage [04/01/14 14:00]: PolarPlotを実行しても、なぜか命令文と同じ文が表示されて、
グラフが描かれません。
どのようにすればグラフが描かれるようになるのでしょうか？
教えてください。お願いします。
Version5です。
903 名前：837 mailto:sage [04/01/15 01:36]: << Graphics`Graphics`
904 名前：902 mailto:sage [04/01/16 16:53]: >>903
さんくす
905 名前：１３２人目の素数さん mailto:sage [04/01/16 17:53]: バックグラウンドを黒にすることはできないの？
906 名前：837 mailto:sage [04/01/16 18:38]: Alt+E,F
書式設定→フォント設定→(FontColor || Background)
907 名前：１３２人目の素数さん mailto:sage [04/01/16 19:24]: さんくす。まじで。
908 名前：１３２人目の素数さん mailto:sage [04/01/17 20:33]: InverseLaplaceTransform[ (12s - 6)/(4s^2 - 4s + 37) , s , t]
を使って逆ラプラス変換を求めてみたところ次のようになりました。

3/2 * ( Ｅxp（αｔ）＋Ｅxp(βｔ））　[α＝1/2 - 3i　β＝1/2 + 3i]

でもじぶんで計算したら、頭についてる係数が、2/3 でなくて 6 になった。。。

mathematicaのほうが合ってるんやろか？
mathematicaの計算結果って絶対正しいんかな～

ちなみにバージョンは４でつ。
909 名前：844 mailto:sage [04/01/17 20:52]: いつもお世話になっております

ある評価をした際にグラフが2つ出力される場合，その2つをそれぞれ2つのピクチャボックスに表示させたいのですが，
どうすればよいでしょうか？

描画関連の関数はとりあえずVCFEのを流用して，改変してみたのですが，
1つ目のグラフをピクチャボックス1に表示させ，2つ目をピクチャボックス2に表示させようとすると，1つ目が表示され，その後2つ目が表示されるのですが
1つ目がなぜかクリアされてしまいます・・・

OpenWMFFile関数内のDeleteEnhMetaFileを実行しなければクリアされないかと思ったのですが，うまく行きませんでした

スレ違いかもしれませんが，VCFE関連ですのでご容赦ください
910 名前：837 mailto:sage [04/01/18 00:07]: >>908
逆ラプラス変換についてはわからないど素人ですが、
私のmathematica4で計算させましたところ、こんなんでました。
3/2 E^((1/2 - 3I) t) (1 + E^(6I t))

>>909
だって保存しているのが同じ名前だから
二つめのWMFFileが上書きされてしまうでしょう。

受け取ったWMFFileを保存するさいの名前を別のにしましょう。
911 名前：１３２人目の素数さん [04/01/18 00:56]: お願いがあります。
どなたかSech[at]のフーリエ変換をしてください。
FourierTransform[Sech[at], t, \[Omega]]
912 名前：837 mailto:sage おらおらおらおら [04/01/18 01:57]: In[6]:=
FourierTransform[Sech[t], t, ω]

Out[6]=
\!$\@\(π\/2$\ Sech[$π\ ω$\/2]\)

In[7]:=
FourierTransform[Sech[a t], t, ω]

Out[7]=
\!$\(\@\(π\/2$\ Sech[$π\ ω$\/$2\ a$]\)\/a\)
913 名前：１３２人目の素数さん mailto:sage [04/01/18 07:27]: フーリエ変換の係数は指定できます。ヘルプ参照。
914 名前：908 mailto:sage [04/01/18 13:06]: >>910
ごちゃごちゃ検算とかしてみたところ、ただのじぶんの計算ミスだった
みたいです。mathematicaマンセー。

ラプラス変換を使うと微分方程式が簡単に解ける。（宿題）

１．DSolveで解く
　　DSolve[　{4y''[t] - 4y'[t] + 37y[t] == 0, y[0] == 3, y'[0] == 3/2}　, y[t]　, t　]

２．InverseLaplaceTransformしてFullSimplifyしてみる
　　InverseLaplaceTransform[　(12s - 6)/(4s^2 - 4s + 37)　, s　, t　]
　　FullSimplify[%]

結果
３Exp(t/2) * Cos(3t)

ってかんじで2通りの方法で検算してみますた。
915 名前：908 mailto:sage [04/01/18 13:15]: (12s - 6)/(4s^2 - 4s + 37)を部分分数分解するところで

Ａ／（ｓ－α）　＋　Ｂ／（ｓ－β）とやったのがまずかったみたいです。
分母の最高次の係数を１にしないとダメぽです。

1/4 * [Ａ／（ｓ－α）　＋　Ｂ／（ｓ－β）]が正しい。

部分分数分解なんて大嫌いだ～～

ttp://www12.plala.or.jp/ksp/formula/mathFormula/html/node60.htmlの
”三つにわかれるもの 1”が参考になりますた。インターネット、マンセー
916 名前：911 [04/01/18 14:25]: >>912
感謝。
917 名前：844 mailto:sage [04/01/18 17:35]: >>837さん
レスありがとうございます
保存というのはどこで行っているのでしょうか？

fopen関数でメタデータが入ったtempファイルを読み出して
SetWinMetaFileBitsでメタファイルを作成しているんですよね？
その作成したメタファイルを保存しているのでしょうか？
自分の見たところどこでは保存しているのかちょっとわからないです・・・
918 名前：１３２人目の素数さん [04/01/19 12:59]: コピーに厳しいMathematicaより
コピーに優しいMapleだろ、やっぱ
919 名前：837 mailto:sage おらおらおらおら [04/01/20 02:09]: >>917
きっとテンポラリファイルを作る関数で作ってるよ。
fopenで開いているのであれば、fopenで渡す引数が生まれたところだ。

生まれた場所へ誰もが還るのよ……。

>>918
Mapleは>>894,897,898の問題は解けますか。
解けるのであれば速攻移行します。またMathLinkみたいな
プログラムがMapleを呼び出すインターフェイスはありますか？
920 名前：844 mailto:sage [04/01/20 21:01]: >>919さん
レスありがとうございます
たびたびすみません

tempファイルを作る関数ってありますでしょうか？
DISPLAYENDPKTが返ってきたときにtempファイルのフルパスが入ったOutputStringを取得し，
それをOpenWMFFile関数のfopenにおくっていると思うのですが。
ちなみにフルパスを確認したところ，tempファイルは2つ作られていて上書きされてはいませんでした
921 名前：837 mailto:sage おらおらおらおら [04/01/21 00:19]: >>920
たぶん、きちんと二つのWMFファイルをメモリ上に読みこんでいないのでは。

いちおう二つめのWMFファイルを手に入れる前に
新たにtempfileとして一つ目のファイルをコピーしてから、メモリ上に読みこんでおき、
そして二つめのWMFファイルも同じく
新たにtempfileとしてコピーしてからメモリ上に読みこんでおき、
描画メッセージが送られてきた時に表示するようにするのが言いと思います。
ピクチャボックスならHANDLEとして渡している部分がきっとあるでしょう。
そして必要なくなったらHANDLEを閉じて、読みこんだメモリを削除して、
コピーしたtempfile達も削除すると。
922 名前：１３２人目の素数さん mailto:sage [04/01/21 11:45]: science2.2ch.net/test/read.cgi/math/1074420890/323です。誘導されてきました。

mathematicaを使って解かなければならないのですが
------------
番号iのつけられた壺が無限個ある。（ i = 0, 1, 2, 3, ... ）
色だけが違う同じ大きさの赤玉と黒玉が無限個あり，0番目の壺にはnTotal個の玉が入っている。
その内nRed0個は赤玉である。
残りの壺の中にそれぞれnTotal個の玉を順次，次の手順で入れる。
［手順］
i 番目の壺に玉を入れるときには，( i - 1 ) 番目の壺から nTotal 回玉を復元抽出し，
抽出するたびに抽出された玉と同じ色の玉を 1 個選び i 番目の壺に入れる。

試行のたびごとに確率的に定まるある番号iFixがあり，その番号iFixから先の壺の中の玉はすべて同じ色になる。
このことを色固定と呼び，最初に色固定が起こる壺の番号iFixを色固定番号と呼ぶ。
また，色固定が起こった壺の中の玉の色fixColorを固定色と呼ぶ。
固定色fixColorと色固定番号iFixはnTotalとnRed0を与えると定まる確率変数である。
その分布をmathematicaによって調べよ。（玉入れ試行を多数回行い，その結果を統計的に整理せよ。）

ヒント：
・一回の玉入れ試行は，色固定が起これば，その壺で終了してよい。
・固定色fixColorが赤である確率pRedを求めよ。
・色固定番号iFixの平均meanIFixと分散varIFixを求めよ。
・nTotalは10 ～ 100程度を考えよ。
nRed0は0からnTotalまで変化させることができる。結果はnRed0にどのように依存するか調べよ。
-------------
ご教授よろしくお願いします。
923 名前：１３２人目の素数さん mailto:sage [04/01/21 18:33]: どこまでわかってるの？
924 名前：１３２人目の素数さん mailto:sage [04/01/21 19:00]: >>922
実験のフレームワークと1番目の解答まで書いてみよう。

(* GetBall : get a ball from the vase given as the args. *)
GetBall[Vase[nTotal_Integer, nRed_Integer]] :=
If[Random[Integer, {1, nTotal}] <= nRed, Red, Black]

(* Onestep : {n, n-th vase} -> {n+1, (n+1)-th vase} *)
OneStep[{nth_Integer, v : Vase[nTotal_, nRed_]}] := {nth +
1, (#[v] &) /@ Table[GetBall, {i, nTotal}] //
Vase[nTotal, Count[#, Red]] &}

(* Trial returns {iFix, fixColor} *)
Trial[nTotal_, nRed0_] :=
NestWhile[OneStep, {1, Vase[nTotal, nRed0]},
#[[2, 2]] != 0 && #[[2, 2]] != nTotal &] /. {Vase[nTotal, 0] -> Black,
Vase[nTotal, nTotal] -> Red}

(* pRed *)
pRed[trys_][nTotal_, nRed0_] :=
Count[Table[Trial[nTotal, nRed0], {trys}], {_, Red}] / trys // N
925 名前：924 mailto:sage [04/01/21 19:07]: 機能に影響は無いがちょっと修正
args -> arg
{i, nTotal} -> {nTotal}
926 名前：１３２人目の素数さん mailto:sage [04/01/21 23:27]: 正直、>>924さんのを見て驚愕しています。
mathematicaの基本的な使い方しか授業では習っていなかったもので....。
この問題も解法の見当がつかない次第でした。
明日大学で>>924さんのを入力してみます。ありがとうございます
927 名前：837 mailto:sage [04/01/22 14:00]: 面白いなぁ。おいらには答えられない。
928 名前：844 mailto:sage [04/01/27 17:14]: >>837さん
遅レスですがどうもありがとうございました
tempファイルをコピーしておいて，ピクチャボックスのプロシージャや
メタファイルハンドル2つ用意してやってみたらなんとかできました

また何かありましたらよろしくです＾＾；
929 名前：１３２人目の素数さん mailto:sage [04/01/27 20:42]: mapleｽﾚがないのでこちらで質問させていただきます。
mapleのワークシートの背景を変更することはできないのでしょうか。
白だと目が痛いです・・・。
930 名前：１３２人目の素数さん mailto:sage [04/01/27 22:43]: Mathematica4で、ListPlotでx軸が0から1までの2次元の折れ線グラフを書いているのですが、
x軸とy軸の交点の部分の0を表示させたいです。
それと、x軸の数字というか、刻みが0.2刻みで、0.2,0.4,0.6,0.8,1とグラフが刻まれているのですが、
これを0,0.5,1の、0.5刻みにしたいです。
どのようなパラメータを与えれば、これらができるのでしょうか？
よろしくお願いします。
931 名前：837 mailto:sage [04/01/28 00:45]: >>928
ｶﾞﾝ(　ﾟдﾟ)ｶﾞﾚ

>>929
ヘルプファイルを熟読した後、
Mapleのスレを建ててください。
もしくは数学ソフト一般（Mathematica以外？）のスレ。

>>930
(* x の目盛りだけを指定する *)
SetOptions[ListPlot, Ticks -> {{0, 0.5, 1}, Automatic}]
(* ↑だと０が表示されないので、原点を少しずらす *)
SetOptions[ListPlot, AxesOrigin -> {-0.0001, -0.0001}]

(* またこれでもいい *)
ListPlot[{{0.0, 0.0}, {0.1, 0.5}, {0.3, 0.3}, {0.5, 0.4}, {1.0, 0.6}},
PlotJoined -> True, AxesOrigin -> {-0.0001, -0.0001},
Ticks -> {{0, 0.5, 1}, Automatic}]

(* こんなOptionも発見した *)
tl = ListPlot[{5, 3, 4, 6}, DisplayFunction -> Identity];
Show[tl, DisplayFunction -> $DisplayFunction];
932 名前：１３２人目の素数さん mailto:sage [04/01/28 01:26]: 深さ優先探索
擬似コード

局所変数：g,start,edgeOrdering,adjList,checked

search(v)
checked[[v]]←True
for each u∈N(v) do
if checked[[u]] = False then
Append(edgeOrdering,{v,u})
search(u)

depthFirstSearch(G,start)
for each vertex v∈V do
checked[[v]]←False
search(start)
Return edgeOrdering

depthFirstSearchで定義される局所変数はsearchで使用する。
そのためdepthFirstSearchはBlockとする必要がある。
隣接リストを使った方が楽である。

上記の擬似コードのMathematicaのプログラムを教えて下さいm(_ _)m
933 名前：１３２人目の素数さん mailto:sage [04/01/28 01:57]: >>931
オオ！完璧です！
あいらぶゆー
934 名前：１３２人目の素数さん mailto:sage [04/01/28 16:22]: スミマセン、私もListPlotで点を連結して2Dの折れ線グラフを書いてるんですが、
複数のグラフを重ねてShowで表示させ、グラフ同士を比較したいのですが、
どれがどのグラフなのかがわかりにくくなってしまいます。
そこで、各グラフの右端、つまりx軸の最大値に対してPlotした点の右脇に、
グラフの番号を表示させて、グラフを重ねてもグラフ同士の区別がつくようにしたいのですが、
どうやればいいんでしょうか…。
親切な方、教えて下さいm(＿)m
935 名前：１３２人目の素数さん mailto:sage [04/01/28 22:36]: >>934
ていうか色をつければいいのでは。
936 名前：１３２人目の素数さん mailto:sage [04/01/28 23:02]: >>935
印刷するのでカラー使えないんですよ…。
いよいよ方法がないとなったら、MSペイントか何かで直接画像に数字書き込めばいいんで
そんなせっぱ詰まってもないんですけど…。
937 名前：１３２人目の素数さん mailto:sage [04/01/29 14:45]: >>936
じゃあ、Legendをヘルプで。
938 名前：１３２人目の素数さん mailto:sage [04/01/29 19:22]: >>870

ちなみにkeygenはネットに出回ってるよ。5.0でも使える奴が。
939 名前：１３２人目の素数さん mailto:sage [04/01/29 22:50]: オイラー閉路の問題なんですが、

アルゴリズム

入力：オイラー閉路を持つグラフＧ
出力：Ｇのオイラー閉路
手順：任意の頂点Ｓから出発して、次の規則に従いたどれる辺が無くなるまで辺をたどる
　　　（１）たどった辺はＧから除く
　　　（２）ほかにたどる辺が無い場合以外は橋をだどらない。

上のオイラー閉路を見つけるためのアルゴリズムを使った、アニメーションプログラムを教えて下さい！
940 名前：１３２人目の素数さん [04/01/31 00:14]: 有限体GF(2)上の拡大体GF(2^8)の256個(0,1,α,α^2,・・・,α^254)の元
を8次元ベクトル表示した表が欲しいのですがどうすればいいか分かりません。
教えてください。よろしくお願いします。
(αは原始元，既約多項式はX^8+X^4+X^3+X^2+1です)
941 名前：837 mailto:sage [04/01/31 01:25]: >>932
擬似コードってなんですか。

>>939
この辺かな
mathworld.wolfram.com/EulerianCircuit.html
mathworld.wolfram.com/KoenigsbergBridgeProblem.html

>>940
有限体がよくわからないので、とりあえずここで検索してみてください。
mathworld.wolfram.com/

手計算で求められるのであれば、式を書きこんでおいてください。
誰かがmathematica上のプログラムに変換してくれます。
942 名前：940 [04/01/31 02:28]: >>941
アドバイス感謝します。
α^0=(0,0,0,0,0,0,0,1) α^1=(0,0,0,0,0,0,1,0) α^2=(0,0,0,0,0,1,0,0)
α^3=(0,0,0,0,1,0,0,0) α^4=(0,0,0,1,0,0,0,0) α^5=(0,0,1,0,0,0,0,0)
α^6=(0,1,0,0,0,0,0,0) α^7=(1,0,0,0,0,0,0,0)
α^8はα^8+α^4+α^3+α^2+1=0より，α^8=α^4+α^3+α^2+1
α^8=(0,0,0,1,1,1,0,1)
α^9=α^5+α^4+α^3+αより
α^9=(0,0,1,1,1,0,1,0)
とずっと計算していくとα^255=1
となります。こういうのmathematicaの関数で計算できないですかねえ？
943 名前：837 mailto:sage [04/01/31 19:29]: α^ 8=α^4+α^3+α^2+1 =(00011101)
α^ 9=α^5+α^4+α^3+α^1=(00111010)
α^10=α^6+α^5+α^4+α^2=(01110100)
α^11=α^7+α^6+α^5+α^3=(11101000)
α^12=α^8+α^7+α^6+α^4
　　 =α^7+α^6+2*α^4+α^3+α^2+1=(11010000)+(00011101)
　　 =(11011101) || (11001101)

α^12の答えはどっちだろう。
どちらにしろビット演算の繰り返しっぽいですね。
944 名前：１３２人目の素数さん mailto:sage [04/02/02 11:37]: アニメーションGIFについて質問です．
　
　Export["filename",graphics,"GIF"]

で出力はするのですが，それをPOWER POINTに貼り付けても動いてくれません．
Windows　画像とFAXビューアでも「描画できませんでした」と表示されます．

Internet Explorerで開いたときはちゃんと動くのですが．．．（しかし，ループはしません）

どうすればPOWER POINTで動くようにできるのでしょうか？
ご存知の方よろしくお願いします．
945 名前：837 mailto:sage [04/02/03 00:23]: >>944
animation GIFのloop flagがたっていないだけだから、
手ごろなanimation gif editorで、flagを立てた上に
スレ違いなのでヘルプ読め
GIF アニメーションの図をスライドに挿入する
946 名前：944 mailto:sage [04/02/03 12:37]: >>837
スレ違いなのに答えてくれてありがとうございます．
今からエディターでやってみます．
947 名前：１３２人目の素数さん mailto:sage [04/02/05 23:23]: Mathematicaって他のソフトと違ってバグフィックスのアップデータないけど
バグらしいバグなんてないのかな？
それはそれですごいと思うけど。
948 名前：１３２人目の素数さん mailto:sage [04/02/05 23:42]: >>947
あるけど直さないだけ。おれも何個も見つけた。
949 名前：１３２人目の素数さん mailto:sage [04/02/07 16:10]: >>948
へ～
どんなのです？
再現性はありますか？
試してみたい。
950 名前：１３２人目の素数さん mailto:sage [04/02/08 15:07]: >>940
単にAlgebra`FiniteFields`を使えばいいのでは？
ヘルプ読んでみた？
951 名前：１３２人目の素数さん [04/02/08 15:27]: >>940

<< Algebra`FiniteFields`
fld = GF[2, {1, 0, 1, 1, 1, 0, 0, 0, 1}]
α = fld[{0, 1, 0, 0, 0, 0, 0, 0, 0}]
Table[α^n, {n, 255}]
952 名前：１３２人目の素数さん mailto:sage [04/02/08 21:24]: >>949
個人情報があれなので、俺(も発見した)が報告したわけではないやつを。
Integrate[D[f[x, y, t], x, t] + D[f[x, y, t], y, t], t]
が f^(0,1,0)[x,y,t] を返す(V4.2)
953 名前：１３２人目の素数さん mailto:sage [04/02/09 11:33]: f^(0,1,0)[x,y,t]+f^(1,0,0)[x,y,t]を返す(V3.0)
954 名前：mathematicant [04/02/09 16:10]: >952
>個人情報があれなので

つまり、御禁制品を使用していると・・・
955 名前：１３２人目の素数さん mailto:sage [04/02/09 16:24]: >>954
どういうロジックでそうなるんだ？
レポートしたから、おれが誰だか分かってしまうと言う意味だ。
956 名前：837 mailto:sage [04/02/10 08:10]: >>950,951
( ･∀･)つ〃∩ ｽｹﾞｪｰｽｹﾞｪｰｽｹﾞｪｰｽｹﾞｪｰｽｹﾞｪｰ
出てくるまでちょっと時間がかかるのと、bit列が左と右が逆なので注意。
957 名前：１３２人目の素数さん mailto:sage [04/02/13 18:48]: >>952
へ～
バグですね
ほかのバグありますか？
958 名前：１３２人目の素数さん mailto:sage [04/02/14 22:38]: Table[f[i], {i, 4}]　－>　{f[1],f[2],f[3],f[4]}
ですが、
中カッコを外した形式
f[1],f[2],f[3],f[4]
を得るにはどうすればいいですか？
959 名前：１３２人目の素数さん mailto:sage [04/02/15 01:23]: >>958
Sequence@@Table[f[i], {i,4}]

(* Example *)
Plus[Sequence@@Table[f[i], {i,4}]]
960 名前：１３２人目の素数さん mailto:sage [04/02/16 20:25]: うちのMathematica5.0ではこの式を評価すると確実に落ちるのだが、
バグだろうか。

7910*Csc[(4*Pi)/15]*Sqrt[((-(1 -
Cos[82000/797267]^2))*Sin[(4*Pi)/15]^2 +Sqrt[(1 -
Cos[82000/797267]^2)^2*Sin[(4*Pi)/15]^4 - (-1 +
Cos[82000/797267]^2)*Sin[(4*Pi)/15]^2*(Cos[82000/797267]^2
- Sin[(4*Pi)/15]^2)])/(Cos[82000/797267]^2 -
Sin[(4*Pi)/15]^2)]//FullSimplify
961 名前：１３２人目の素数さん mailto:sage [04/02/16 23:06]: >>960
うちの5.0でもカーネルが落ちる。
なんでかねぇ。
フロントエンドは落ちないけど。
962 名前： ◆BhMath2chk mailto:sage [04/02/17 16:00]: ＭｕＰＡＤ使用開始から三日。
963 名前：１３２人目の素数さん [04/02/19 11:16]: 教えてください。

Mathematica5.0を使っています。
Eigenvalueを使って固有値を計算させているのですが、結果が「#1」といった
引数を含む表現になってしまい，計算結果を利用して解析を進めることができません。

結果の中の引数を展開する方法をお教え願えないでしょうか？
964 名前：１３２人目の素数さん mailto:sage [04/02/21 00:36]: >>963
変数でも入ってたの？
何をやったか書いてよ。
965 名前：837 mailto:sage [04/02/22 15:46]: >>963
そうだそうだ。
どの式がうまくいかないのかinputの式をここに書いてくれ。
君の姿は見えないんだ。help meと聞こえてもどう助けていいかわからない。
966 名前：１３２人目の素数さん mailto:sage [04/02/22 17:12]: >>963
ふと君がなにをいってるのか推測できた。
5次以上の行列で、固有多項式の根がRootオブジェクトになったんだろうな。
あれはRoot[関数, 引数]というもので、全体は関数じゃなくて値。
数値が欲しいならEigenvalue[A] // Nとかしとけ。
967 名前：１３２人目の素数さん [04/02/26 05:25]: 多項式のGalois群を計算するMathematicaのパッケージってありますか？
Mapleだとgalois()って関数があるんだが。
968 名前：１３２人目の素数さん mailto:sage [04/02/28 09:15]: >>967
Ver.5のHelpを検索すると標準パッケージでついているようだが。
Galois fields, Algebra`FiniteFields` , Algebra`PolynomialPowerMod`
969 名前：１３２人目の素数さん [04/02/29 15:24]: >>938
keygenってlinuxでも動きますか？
970 名前：１３２人目の素数さん [04/02/29 19:54]: >Galois fields
Ver.4にもある。が整式のがろあ群を出してくれるかはしらん
ただの有限体じゃないの

>keygen
win用じゃ。linuxで動かないだろ、えみゅでも使え

MathematicaのStudent版、アカデミック版、Highschool版だかは
バージョンが二つ前のやつは、Math IDに対するパスワードを
発行してくれないようだな。アカ版はぢかうか
以後、自分で買うのはよして、大学のをつかうことにしよう。
971 名前：１３２人目の素数さん [04/02/29 20:11]: >>968
Galois field（= Finite field = 有限体）
≠Galois group
なんですが･･･｡代数の教科書読んでね｡
で、結局、Galois群を計算するパッケージってあるの？
（標準のものでも誰かが自分で作ったんでもいいんだが）
972 名前：968 mailto:sage [04/02/29 20:29]: >>971
すまんね、勉強不足で。
library.wolfram.com/infocenter/Articles/2872/
973 名前：１３２人目の素数さん [04/02/29 23:18]: >>972
情報サンクス。
だけど、ざっと見たところこれは、
Mathematicaを使って特定の多項式のGalois群を頑張って計算してみました、
って感じのもの見たいね・・・
Maple （俺が持ってるのはV5R4）だと
> galois(x^5 - 5*x + 12);
　+D5, 10, {(1 2 3 4 5), (2 5)(3 4)}
こんなふうに、任意の（7次以下の）多項式のQ上のGalois群をさくっと
出してくれるんだが・・・
974 名前：１３２人目の素数さん mailto:sage [04/03/01 00:04]: >>973
5次なら、
library.wolfram.com/infocenter/TechNotes/158/
にあるが。
MathematicaはMapleとは方向性が違うからな。
ライブラリにある機能しかつかわない人にはMapleのほうがいいだろう。
975 名前：１３２人目の素数さん mailto:sage [04/03/01 00:26]: >>973
www.google.com/search?num=100&hl=ja&ie=UTF-8&oe=UTF-8&c2coff=1&q=%22Galois+groups%22+mathematica+package&btnG=Google+検索&lr=
976 名前：１３２人目の素数さん [04/03/01 00:41]: >>975
googleで検索してもガイシュツのものしか見つからないんだが。
977 名前：１３２人目の素数さん mailto:sage [04/03/01 01:17]: 5のアイコンより4.2のアイコンの方がよくない？
978 名前：１３２人目の素数さん [04/03/01 19:17]: どっちでもいいから
せこいコピー風呂テクと止めれ

正規ユーザがばかをみてるぞ　Wolfram
979 名前：１３２人目の素数さん mailto:sage [04/03/01 21:16]: カオス理論の本読んでたらこのソフトの作者がでてきた。
意外なつながりだ。
980 名前：１３２人目の素数さん mailto:sage [04/03/03 02:38]: >>978
どういうこと？
981 名前：１３２人目の素数さん mailto:sage [04/03/06 09:51]: 次スレですよ～

〓 Mathematica 2 〓
science2.2ch.net/test/read.cgi/math/1078534285/
982 名前：１３２人目の素数さん mailto:sage [04/03/07 09:51]: ９８２。
983 名前：１３２人目の素数さん mailto:sage [04/03/08 09:51]: ９８３。
984 名前：１３２人目の素数さん mailto:sage [04/03/08 10:15]: 984
985 名前：１３２人目の素数さん [04/03/08 11:32]: age
986 名前：１３２人目の素数さん mailto:sage [04/03/08 15:49]: Maple >>>(超えられない壁) >> Mathemetica >>> Maxima >> mupad
987 名前：１３２人目の素数さん mailto:sage [04/03/08 22:23]: うめー
988 名前：１３２人目の素数さん mailto:sage [04/03/08 22:24]: うめー
989 名前：１３２人目の素数さん mailto:sage [04/03/08 22:25]: うめー
990 名前：１３２人目の素数さん mailto:sage [04/03/08 22:26]: うめー
991 名前：１３２人目の素数さん mailto:sage [04/03/08 22:27]: うめー
992 名前：１３２人目の素数さん mailto:sage [04/03/08 22:28]: うめー
993 名前：１３２人目の素数さん mailto:sage [04/03/08 22:29]: うめー
994 名前：１３２人目の素数さん mailto:sage [04/03/08 22:30]: うめー
995 名前：１３２人目の素数さん mailto:sage [04/03/08 22:31]: うめー
996 名前：１３２人目の素数さん mailto:sage [04/03/08 22:32]: うめー
997 名前：１３２人目の素数さん mailto:sage [04/03/08 22:33]: うめー
998 名前：１３２人目の素数さん mailto:sage [04/03/08 22:34]: うめー
999 名前：１３２人目の素数さん mailto:sage [04/03/08 22:35]: うめー
1000 名前：１３２人目の素数さん mailto:sage [04/03/08 22:36]: うめー
1001 名前：１００１ [Over 1000 Thread]: このスレッドは１０００を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef