【数式処理システム】 Maxima 【Common Lisp】

1 名前：１３２人目の素数さん mailto:sage [2008/09/02(火) 16:31:35 ] フリーの数式処理システムMaximaに関するスレです。 MaximaはCommon Lispで書かれており、GNU GPLライセンスのもと 配布されています。各種UNIX、Windows、OSXで動作します。 ■関連サイト □Maxima, a Computer Algebra System (英語) - 公式サイト 　ttp://maxima.sourceforge.net/ □Maxima Beginner's FAQ (英語) - 初心者向けのFAQ 　ttp://beshenov.ru/maxima/faq.html □Imaxima (英語) - Emacs用のGUIフロントエンド 　ttp://members3.jcom.home.ne.jp/imaxima/Site/Welcome.html □Maxima interest list (英語) - Maximaに関するメーリングリスト 　ttp://www.ma.utexas.edu/mailman/listinfo/maxima □Maximaで遊ぼう (日本語) - マニュアルの日本語訳や、インストールガイド等 　ttp://www.bekkoame.ne.jp/~ponpoko/Math/maxima/MaximaMAIN.html □数式処理システムMaximaで楽をしよう (日本語) - 初心者向けのMaxima導入編 　ttp://www1.bbiq.jp/kougaku/maxima.html ■過去スレ MAXIMA science6.2ch.net/test/read.cgi/sim/1011102458/ 　ttp://www.23ch.info/test/read.cgi/sim/1011102458/ (ミラー) Maxima part2 science6.2ch.net/test/read.cgi/sim/1218342679/ 2 名前：１３２人目の素数さん mailto:sage [2008/09/02(火) 16:41:48 ] >>1 スレ立て乙です 3 名前：１３２人目の素数さん mailto:sage [2008/09/02(火) 17:40:31 ] (Common Lisp処理系) sbcl (32bit/64bit) ネイティブコンパイラで高速cmuclから 　　　　　　　　　　　　分岐した処理系活発に改良している。 www.sbcl.org/ cmucl (32bit) ネイティブコンパイラで高速 www.cons.org/cmucl/ clisp (32bit/64bit) 低速なcommon lisp (処理系さまざまプラットホームで動く) clisp.cons.org Clozure CL (64bit) + 32bit PPC maximaは不完全だが対応 trac.clozure.com/openmcl ECL (32bit/64bit) ver 0.9I よりMaximaに正式対応 こちらもさまざまなプラットホームで動く ecls.sourceforge.net/ おまけ Maxima 普及委員会 www.cymric.jp/ 4 名前：１３２人目の素数さん [2008/09/02(火) 18:52:26 ] まずは基本から教えてくれ 必須コマンドと使用例よろすく 5 名前：１３２人目の素数さん mailto:sage [2008/09/04(木) 19:26:32 ] Ｒとの連携はlisp経由でできるんでしょ？ matlabは？ 6 名前：１３２人目の素数さん mailto:sage [2008/09/04(木) 19:28:45 ] Lispでプログラム書く気になれば何とでも連携できるんじゃね？ 7 名前：１３２人目の素数さん mailto:sage [2008/09/06(土) 23:06:36 ] 　　　　　　☆ ﾁﾝ　　　　　ﾏﾁｸﾀﾋﾞﾚﾀ〜 　　　　　　　　　　　　　　　　　　　　　　　　ﾏﾁｸﾀﾋﾞﾚﾀ〜 　　 　 　 ☆　ﾁﾝ　　〃　 ∧＿∧　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣ 　 　 　 　　ヽ　＿＿_＼（＼・∀・）　＜　MAXCELまだ〜？ 　 　 　 　 　 　 ＼＿／⊂　⊂＿ )　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿ 　 　 　 　 　 ／￣￣￣￣￣￣ ／| 　　　　　　　|￣￣￣￣￣￣￣|　 | 　　　　　　　|　　　日　立　　　|／ 8 名前：１３２人目の素数さん mailto:sage [2008/09/07(日) 02:45:58 ] >>7 つ www.maxell.co.jp/ 9 名前：１３２人目の素数さん [2008/09/08(月) 08:41:56 ] 懐かしいヌレだ 10 名前：１３２人目の素数さん mailto:sage [2008/09/13(土) 01:05:28 ] ほしゅ 11 名前：１３２人目の素数さん mailto:sage [2008/09/17(水) 17:12:04 ] 5.16.3をsourceからmakeしたらmake checkでベッセル関数の値が期待したのと違う値になるという警告でる 12 名前：１３２人目の素数さん [2008/09/17(水) 17:45:32 ] age 13 名前：１３２人目の素数さん mailto:sage [2008/09/17(水) 18:59:04 ] um 14 名前：１３２人目の素数さん mailto:sage [2008/09/22(月) 23:08:59 ] maxima-5.16.3.exeがウィルスに感染してたようなんだけど、誤作動かな？ ちなみにアンチウィルスソフトはca インターネットセキュリティ 15 名前：１３２人目の素数さん [2008/09/23(火) 21:13:28 ] 保守 16 名前：１３２人目の素数さん mailto:sage [2008/09/23(火) 21:19:18 ] f(x)をx>=0ではx^2、x<0ではx^3みたいな範囲によって変わる関数に定義して、 それをグラフにプロットするにはどうすればよいのでしょうか 17 名前：１３２人目の素数さん mailto:sage [2008/09/23(火) 21:36:58 ] >>16 f(x):=if x >= 0 then x^2 else x^3; plot2d(f(x),[x,-1,1]); 18 名前：１３２人目の素数さん mailto:sage [2008/09/23(火) 22:06:10 ] >>17 できないかと思ったけど最新verにしたらできました ありがとうございます 19 名前：１３２人目の素数さん mailto:sage [2008/09/27(土) 03:01:41 ] imaximaの1.0bでinline plotできるって聞いたんでやってみたんだけど、 これ結構重いな 20 名前：１３２人目の素数さん mailto:age [2008/09/28(日) 11:43:14 ] msysでソースをコンパイル出来るのでしょうか？ USBメモリにいれて持歩きたいんですけど 21 名前：１３２人目の素数さん [2008/09/30(火) 06:38:53 ] wxmaxima のボタンを日本語にしてみたんですが, 誰かやってましたか? ttp://www.cs.tohoku-gakuin.ac.jp/~otofuji/ETC/wxMaxima/ 22 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 07:28:53 ] >>20 msysってなんや？ 23 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 12:17:14 ] ぐぐれかす 24 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 12:19:33 ] >>23 おしえるためにきいてるからぐぐるきなんてないよ ボケ だったら、自分でやれ 25 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 12:26:38 ] >>24 おまえの敵は一人じゃないぞｗｗ 26 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 12:27:36 ] >>24 判らないなら教えなくても良いのだよ 27 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 13:25:26 ] >>26,27 だったら、おまえら協力してやれよ。下衆には関わるつもりないから。 28 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 13:51:20 ] >>27 関わるなよｗ 29 名前：１３２人目の素数さん mailto:sage [2008/09/30(火) 13:51:59 ] つか、そもそも MSYS 知らん奴に答えられる問題じゃなかろうにｗ 30 名前：20 mailto:sage [2008/09/30(火) 20:20:47 ] 仮想レジストリソフトなるもので自己解決しましたのでもういいです 31 名前：１３２人目の素数さん mailto:sage [2008/10/13(月) 19:24:10 ] imaximaをインストールできずに困っております。 ご教示いただければ、幸甚です。 (1) ubuntu 8.04をインストール (2) セキュリティのアップデートを済ませる。 (3) 日本語版セットアップからemacs-env-ja, latex-env-ja, latex-extra-jaを導入 (4) maxima_5.16.3-1_i386.debとwxmaxima_0.7.6-1_i386.debを導入。 (5) imaxima-imath-0.99.tar.gzを解凍。sudo ./configure; make; make installを実施。 (6) sudo texhash実施。また、sudo kpsewhich breqn.styを実施。/usr/share/texmf-texlive/tex/latex/breqn/breqn.styとして、breqnにパスが通っていることを確認。 (7) 自分のディレクトリの.emacsに以下の記載をし、パスを通す。 (push "/usr/local/share/emacs/site-lisp" load-path) (push "/usr/local/share/maxima/5.16.3/emacs" load-path) (autoload 'imaxima "imaxima" "Frontend for maxima with image support" t) (autoload 'maxima "maxima" "Frontend for maxima" t) (autoload 'imath-mode "imath" "Imath mode for math formula input" t) (setq imaxima-use-maxima-mode-flag t) (8) ここまでやってから、emacsを動かして、M-x imaximaと読み込ませますと、 /usr/local/lib/maxima/5.16.3/binary-clisp/lisp.run: error while loading shared libraries: libsigsegv.so.0: cannot open shared object file: No such file or directory Process maxima exited abnormally with code 127 とエラーがでますがコンパイルされている、lisp.runのどこをいじれば良いのかわかりません。そもそもこのライブラリがどこで読み込まれるべきかも併せて、ご教示いただければ幸甚です。 32 名前：１３２人目の素数さん mailto:sage [2008/10/13(月) 20:25:46 ] mixiとのダブルポストご苦労様ｗ 33 名前：１３２人目の素数さん mailto:sage [2008/10/14(火) 07:05:08 ] 2chのシミュ板とダブルポストなら、結構気にしますが、mixiとならば、（一応）クローズド環境と完全なオープン環境との違いですので、ダブルポスト可ではないかと私が判断しました。 34 名前：31 mailto:sage [2008/10/15(水) 21:03:35 ] 何とか動きましたので、mixiの方に内容を記載しました。 35 名前：１３２人目の素数さん [2008/10/25(土) 15:49:50 ] age 36 名前：１３２人目の素数さん [2008/11/15(土) 00:20:40 ] maximaってフーリエ変換できますか？ 37 名前：１３２人目の素数さん mailto:sage [2008/11/15(土) 01:00:14 ] できないわけないだろう 38 名前：１３２人目の素数さん [2008/11/15(土) 11:55:17 ] maximaには標準でフーリエ変換が用意されてますか？ 自分で変換の式かかないとダメ？ 39 名前：１３２人目の素数さん mailto:sage [2008/11/15(土) 17:09:19 ] 無いわけないだろう 40 名前：１３２人目の素数さん mailto:sage [2008/11/16(日) 14:56:35 ] 開発がなんか最近低調になっちゃったみたいだね。そろそろ乗り換え時期かな 41 名前：１３２人目の素数さん mailto:sage [2008/11/16(日) 15:07:08 ] 何に乗り換えるん？ 42 名前：１３２人目の素数さん mailto:sage [2008/11/16(日) 15:28:50 ] 開発者が死んでたのか・・・ 43 名前：１３２人目の素数さん mailto:sage [2008/11/16(日) 15:31:51 ] それはそんな最近の話じゃないぞ。 ML見てると最近また活発だという印象。 44 名前：１３２人目の素数さん [2008/11/22(土) 17:59:52 ] 12月13日（土）に福岡大学で開催されるオープンソースカンファレンス2008 Fukuoka において、「はじめてのMaxima」の著者、横田博史さんによるセミナーが予定されています。 お近くの方は是非お越しください。 www.ospn.jp/osc2008-fukuoka/ 関連情報 www.knoppix-math.org/wiki/index.php?OSC2008Fukuoka 45 名前：１３２人目の素数さん mailto:sage [2008/11/22(土) 22:18:35 ] maximaを使って片対数グラフを作成したいのですが， 方法が分からずに困っております． どなた様かご教授願えませんでしょうか． 46 名前：１３２人目の素数さん mailto:sage [2008/11/23(日) 04:17:28 ] 最後の出力結果のLogをとれば良いだけじゃね？ 47 名前：１３２人目の素数さん mailto:sage [2008/12/07(日) 13:49:41 ] １００変数の関数を x1,x3,x5...についてのみ微分したい という場合はx1,x2..を全部手で書くしかないのでしょうか 48 名前：１３２人目の素数さん mailto:sage [2008/12/07(日) 15:38:43 ] 書くプログラムを書けばよいのでは？ 49 名前：１３２人目の素数さん mailto:sage [2008/12/08(月) 03:10:33 ] wxMaximaの使い方がやっとわかった SHIFT+ENTERだったのね 50 名前：１３２人目の素数さん mailto:sage [2008/12/08(月) 23:49:22 ] >>45 2D plotでgnuplotを選んでオプション欄に set logscale x あるいは set logscale y 51 名前：１３２人目の素数さん [2008/12/13(土) 00:58:35 BE:1748874479-2BP(808)] >>49 うん。 でも、勝手に挿入される閉じカッコを無視して式末へジャンプする方法がわからないんだ。 ＿|￣|○ 52 名前：１３２人目の素数さん mailto:sage [2008/12/13(土) 01:40:37 ] ∫[0.2π]log(5-2cosx)dxを計算しようとして integrate(log(5-4*cos(x)), x, 0, 2*%pi); という入力に対してIs 4*sin(x)^2+4*cos(x)^2-4*cos(x)+1 positive or negative?と返ってきたのだがどういうことです 53 名前：１３２人目の素数さん mailto:sage [2008/12/13(土) 02:45:12 ] 4*sin(x)^2+4*cos(x)^2-4*cos(x)+1　は正なのか負なのか？と聞いておられるのだ 54 名前：１３２人目の素数さん [2008/12/17(水) 02:15:13 ] %rというのが出てきたんですが、これって何の定数ですか？ 55 名前：１３２人目の素数さん mailto:sage [2008/12/18(木) 22:30:04 ] >>1 待望のMaximaスレッドを建ててくれてありがとう！ 56 名前：１３２人目の素数さん [2008/12/21(日) 02:07:35 ] sqrt(2)*sqrt(3) の結果として √6 が欲しいんだけど、 √2√3 になってしまう・・・。 √6にするにはどうしたらいいんだ　ofz 57 名前：１３２人目の素数さん mailto:sage [2008/12/21(日) 02:10:25 ] sqrt(2*3) 58 名前：56 mailto:sage [2008/12/21(日) 02:20:50 ] >>57 ありがとう。 検索したら、見つかった。 sqrt(2)*sqrt(3) rootscontract(%) これで√6 になる。 59 名前：１３２人目の素数さん mailto:sage [2008/12/21(日) 18:41:24 ] a*x+a+b*x+b があって、a*x+a を a*(x+1)に置き換えるにはどうしたらいいんでしょう？ 以下のような感じで処理したいのですが・・・ a*x+a+b*x+b a*x+a = a*(x+1) a*(x+1)+b*x+b ←置き換え b*x+b = b*(x+1) a*(x+1)+b*(x+1) ←置き換え x+1 = c a*c+b*c factor(%)を使えば簡単に因数分解できるのは知ってますが、手で 計算するように順を追って処理したいのです。 60 名前：59 [2008/12/21(日) 20:47:34 ] また自己レス f1 : a*x+a+b*x+b ratsimp(f1,a,b) で、b*(x+1)+a*(x+1)が得られた。 rat(...)でも同様のことができるっぽい。 ただ、 f2 : ratsimp(f1,a) →a*(x+1)+b*x+b f3 : ratsimp(f2,b) →b*(x+1)+a*x+a なんで →b*(x+1)+a*(x+1) にならないんだ・・・orz 61 名前：１３２人目の素数さん mailto:sage [2008/12/22(月) 00:29:22 ] >>60 Mathematica使えば？ 62 名前：59 [2008/12/22(月) 01:15:32 ] >>61 年収400万の末端工員の俺にそんなこと言われても困るんだが orz 貧乏で馬鹿だけど、少しは勉強しようと思って始めたんだ。 そのへん察してくれ。 63 名前：１３２人目の素数さん mailto:sage [2008/12/22(月) 01:52:44 ] 紙と鉛筆がお勧め 64 名前：１３２人目の素数さん mailto:sage [2008/12/22(月) 01:55:22 ] Maximaは仕事で使っていいんだっけ？ライセンス的に 65 名前：１３２人目の素数さん mailto:sage [2008/12/22(月) 01:56:29 ] うんいいよ 66 名前：１３２人目の素数さん [2008/12/22(月) 14:52:32 ] 遊びで使ってもいいんだよ 67 名前：59 [2008/12/23(火) 16:37:56 ] >>63 俺、計算中によく±間違えたり、いつの間にか項が 増えたり減ったりするから紙は合わない。 wmMaxima を最新版にアップグレードしたら使い方が大分違ってて あわてた orz 適応能力が低い俺にとっては結構迷惑な話だｗ さて、今日はどんな計算をして遊ぼうかな・・・。 68 名前：59 [2008/12/23(火) 23:32:27 ] また躓いた orz チェーンルールが再現できない depends(x,u); depends(y,x); f1:diff(x,u); f2:diff(y,x); f3:f1*f2; →('diff(x,u,1))*('diff(y,x,1)) orz 教えてください・・・。 69 名前：１３２人目の素数さん mailto:sage [2008/12/27(土) 17:40:32 ] e^xをMaximaで計算させようとして limit((1+x/n)^n,n,inf) としたのですが Is x positive, negative, or zero? と言われてしまいます。こういった場合、どうすればいいのでしょうか？ どうやってxがゼロだとか負だとか正だとか指示すればいいのでしょうか？ というか、ある程度のパターン数までは自動で場合分けやってくれないんですかね？ 70 名前：１３２人目の素数さん mailto:sage [2008/12/27(土) 18:14:47 ] positive; 71 名前：１３２人目の素数さん mailto:sage [2008/12/27(土) 18:20:09 ] 3パターンしかなくてしかもすべて結果が同じなんだからスルーしてほしいよ・・・ 72 名前：１３２人目の素数さん mailto:sage [2008/12/27(土) 18:29:31 ] いちいち対話的に聞いてくるんじゃなくて \{ 1/x (x=0) ∞　(x≠0) みたいな感じで場合分け済みの結果を出してくれるとありがたいね 73 名前：１３２人目の素数さん mailto:sage [2009/01/02(金) 15:23:36 ] Maximaのlimit計算はあまり賢くない Is x-1 positive, negative, or zero?pos; Is x-2 positive, negative, or zero?pos; Is x-3 positive, negative, or zero?pos; ・・・ みたいなことになる場合も多々あるｗ 74 名前：１３２人目の素数さん mailto:sage [2009/01/03(土) 00:56:24 ] 0から2πで積分しようとしてるのにIs x-π positive, negative, or zero?とかきかれても困る 75 名前：１３２人目の素数さん mailto:sage [2009/01/03(土) 19:24:18 ] >>74 0からπ、πから2πというふうに2回に分けて計算しろということだろう 76 名前：１３２人目の素数さん mailto:sage [2009/01/04(日) 01:13:17 ] 積分区間をわざわざ分けるのか。よく分からん 77 名前：１３２人目の素数さん mailto:sage [2009/01/14(水) 12:11:03 ] eclで動くということはcから呼べるの？ 78 名前：１３２人目の素数さん [2009/01/18(日) 11:47:49 ] >>69 assume(x>0); してから limit((1+x/n)^n,n,inf); では駄目？ 79 名前：１３２人目の素数さん mailto:sage [2009/01/28(水) 13:52:25 ] >>31 mixi見れないのでここに書いてもらえませんか？ 80 名前：１３２人目の素数さん mailto:sage [2009/01/28(水) 20:55:22 ] 「PDAで数式処理してる奴いる？」 pc11.2ch.net/test/read.cgi/mobile/1138982038/l50 で知ったのだけど、Xcasってどうよ？ ttp://en.wikipedia.org/wiki/Xcas ttp://www-fourier.ujf-grenoble.fr/~parisse/giac.html 81 名前：１３２人目の素数さん [2009/01/30(金) 02:39:33 ] age 82 名前：１３２人目の素数さん mailto:sage [2009/02/01(日) 00:07:26 ] >>31 (1) 作成者様からのご回答 shell上でmaxima 5.16.3は正常に動作していますか？terminalから % maxima とやってちゃんと起動するかどうか試してみて下さい。このエラーメッセージを見る限り、imaximaというよりはmaxima自身がうまくインストール出来ていないように見えます。 (2) 私の返事その１ ご指摘の通り、wxmaxima.sourceforge.netから落としたdebパッケージがインストールされておりませんでした。(8.04 Hardyテスト済みと書いてありましたが。。。） 仕方なく、ubuntuのsynapticインストールでは5.13.0しかないため、そちらにダウングレードしましたが、今度は/usr/share/maximaができましたが、emacsフォルダがなく、今度はload-pathを考えなくてはなりません。 だめなら、5.15.0にアップグレードしてみて、様子を見てみます。 ありがとうございました。 (2) 私の返事その２ 稼働ご報告。 結局、wxmaxima.sourceforge.netからダウンロードする、maxima-5.16.3のdebパッケージの内部の依存関係がおかしかったようです。 5.15.0のdebパッケージを探してきて、インストールすると、imaximaが動きました。 どうもありがとうございました。 追伸。 alien --scriptでrpmパッケージをdebパッケージに変換しようとしましたけど、なかなかうまくいかないものです。 83 名前：79 mailto:sage [2009/02/02(月) 16:01:22 ] ありがとうございます 84 名前：１３２人目の素数さん mailto:sage [2009/02/06(金) 15:43:27 ] maximaの内部表現 (MLIST SIMP) ((MPLUS SIMP) $X ((MTIMES SIMP) 3 $Y))) を普通のlispの表現 (+ x ( * 3 y)) に変換する関数はないのでしょうか 85 名前：１３２人目の素数さん mailto:sage [2009/02/06(金) 22:31:00 ] maximaのグラフ(gnuplot)をpsやeps,tgifなどの諸形式で保存したいのですが、 可能ですか？ 86 名前：１３２人目の素数さん mailto:sage [2009/02/07(土) 08:45:29 ] >>85 たとえばこれでepsが吐き出されるよ。 plot2d([x^2], [x,-5,5], [gnuplot_preamble, "set size ratio 1; set zeroaxis;"], [gnuplot_term, ps], [gnuplot_out_file, "C:/temp/plot2d.eps"])$ 87 名前：１３２人目の素数さん mailto:sage [2009/02/08(日) 09:20:15 ] >>84 trans4.lisp の tr-lisp->lisp がそれっぽく見えるけどよくわからない 88 名前：１３２人目の素数さん mailto:sage [2009/02/16(月) 14:50:33 ] デルタ関数を作りたいのですが、Lispが出来ません。 作ってくれる人はいませんか？ 89 名前：１３２人目の素数さん mailto:sage [2009/02/20(金) 14:55:32 ] 行列のsvdがないみたいなのですが 自分で作るしかないのでしょうか 90 名前：１３２人目の素数さん [2009/02/20(金) 16:16:29 ] plot3dを使って直線l: x=2y=z　のグラフを表示させたいのですが、 plot3d([[x, 2*y]], [x,-5,5], [y,-5,5]); としても、うまくいきません。教えてください。 91 名前：１３２人目の素数さん mailto:sage [2009/02/20(金) 18:01:13 ] >>88 d.hatena.ne.jp/niitsuma/20080328/1235120364 92 名前：１３２人目の素数さん mailto:sage [2009/02/24(火) 00:00:18 ] >>89 dgesvd()は？ 93 名前：１３２人目の素数さん mailto:sage [2009/02/24(火) 22:17:57 ] maxima-5.17.0をダウソ、使い方をググルもそもそもInput欄などない。 はい、そこで終了。 94 名前：１３２人目の素数さん mailto:sage [2009/02/25(水) 15:08:39 ] >>92 linearalgパッケージが壊れているのか うちでは読み込めない 95 名前：１３２人目の素数さん mailto:sage [2009/02/26(木) 08:14:44 ] >>94 間違った load(lapack)が動かなかった しかし最新版にしたら動いた 96 名前：１３２人目の素数さん [2009/03/05(木) 17:21:50 ] Maxima 5.17.1 を使っています。(WinXP) ・ Option → Preferences でフォントサイズを指定 (Save Preferences) しても、 その時は有効なのに、次回起動したときはディフォールトに戻ってしまう。 ・(%i1) x : 1/3 $ (%i2) float(x); (%o2) 0.33333333333333 (%i3) fpprec : 6; (%o3) 6 (%i4) float(x); (%o4) 0.33333333333333 こういうものなんですか？　以前は出来たような記憶があるんですが。 有効数字の桁数を指定するには、どうしたらいいのでしょうか。 97 名前：１３２人目の素数さん mailto:sage [2009/03/06(金) 19:13:40 ] お！Maximaスレあったのか シミュ板には無くなってるからもう無いのかと思ったわ >>96 fpprintprec : 6$ fpprecは内部で計算する桁数で、表示する桁数はfpprintprec 98 名前：96 [2009/03/06(金) 20:00:12 ] >>97 ありがとうございました。やってみました。 (%i1) a:2/3 $ (%i2) float(a); (%o2) 0.66666666666667 (%i3) fpprintprec; fpprintprec : 6; (%o3) 0 (%o4) 6 (%i5) float(a); (%o5) 0.6667 あれ、有効数字の桁数ではなくて、print prec なのかー。 99 名前：98 [2009/03/06(金) 20:08:33 ] (%i1) fpprintprec; (%o1) 0 (%i2) a:2/3*1e-9; (%o2) 6.6666666666666664E-10 (%i3) float(a); (%o3) 6.6666666666666664E-10 (%i4) fpprintprec : 6 $ (%i5) float(a); (%o5) 6.666667E-10 ナンかヘンです。 100 名前：１３２人目の素数さん mailto:sage [2009/03/06(金) 20:54:58 ] >>98 どうヘンなのか言ってくれないと答えようがない マニュアルによるとfpprintprecはbigfloat型数値のための変数となってるが これで解決するか？ x:1/3;for i:2 thru 16 do block([],fpprintprec:i,print(i," ",float(x)," ",bfloat(x))); 101 名前：１３２人目の素数さん mailto:sage [2009/03/06(金) 23:40:18 ] >>100 ありがとうございます。 x:2/3*1e-6 $ for i:2 thru 16 do block([], fpprintprec:i, print(i, " ", float(x), " ", bfloat(x))); 2　 6.67E-7　 6.6b-7 3　 6.667E-7　 6.66b-7 4　 6.6667E-7　 6.666b-7 5　 6.66667E-7　 6.6666b-7 6　 6.666667E-7　 6.66666b-7 7　 6.6666667E-7　 6.666666b-7 8　 6.66666667E-7　 6.6666666b-7 9　 6.666666667E-7　 6.66666666b-7 10　 6.6666666667E-7　 6.666666666b-7 11　 6.66666666667E-7　 6.6666666666b-7 12　 6.666666666667E-7　 6.66666666666b-7 13　 6.6666666666667E-7　 6.666666666666b-7 14　 6.66666666666667E-7　 6.6666666666666b-7 15　 6.666666666666666E-7　 6.66666666666666b-7 16　 6.666666666666666E-7　 6.666666666666666b-7 とりあえずは bfloat でなく、float 型の方に注目しています。 ・有効数字の桁数の指定が思ったようにできない。自分で関数を作れ！ということなのか。 ・数値の丸め方がおかしい。例「2　6.67E-7　6.6b-7」 　あ、四捨五入の指定のしかたは、別にあるんでしょうか。 全然別件ですが、"," の後にはスペースを１つ入れる、というのが通常なんですが。 「Maxima is a computer program for doing mathematics calculations, symbolic manipulations, numerical・・」これだって、そうなんですが。 102 名前：101 mailto:sage [2009/03/06(金) 23:52:01 ] 言い忘れです。 (%i14) x=2/3 $ fpprintprec : 4 $ float(x); (%o16)　 0.67 「0.67」って４文字だ、と言い張るのでしょうか。 103 名前：101 mailto:sage [2009/03/06(金) 23:58:43 ] MAXIMA は数式処理プログラムであるからして、数値演算は２の次３の次、 というのはわかります。 しかし、近似多項式の演算とかにも使いたいのです。つまり多項式を想定し、 数値ペアを与えて、solve で各次の係数を求めたりする用途です。 結果をうまく丸めてくれないと、ちょっと困ります。 104 名前：１３２人目の素数さん mailto:sage [2009/03/07(土) 00:53:37 ] >>101-103 疲れる人だな Maximaなんかやめて、Octaveを調べた方がいいと思うわ >bfloat でなく、float 型の方に注目しています。 せっかくfpprintprecはbigfloat型数値のための変数と書いたのに・・ マニュアル読んでないだろ 105 名前：１３２人目の素数さん mailto:sage [2009/03/07(土) 08:53:18 ] >>104 ありがとうございます。 説明書のキーワード検索機能が弱いように思います。 bfloat にしても丸め方がおかしい。>>101 fpprintprec が、なぜ float の表示桁数に影響するのか。>>101 Option → Preferences でフォントサイズを指定 (Save Preferences) しても、 その時は有効なのに、次回起動したときはディフォールトに戻ってしまう。>>96 さようなら。 106 名前：１３２人目の素数さん mailto:sage [2009/03/09(月) 11:55:55 ] to_lisp() してから浮動小数演算してみるとlispのデフォルト浮動小数計算と違うのが動いてるのがわかる 107 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 04:08:48 ] >>101 その bfloat のヤツ。小中学生の頃から、そんな書き方をしてはいけません！ と言われて育ったぜ。 108 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 04:19:28 ] カシオ計算機とか、インテル・マイクロソフトの数値演算は さすがにスゴイぞ。文句のつけようがないほどやってる。 109 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 04:21:48 ] ガウスさんだっけ。数値演算ガンガンやって・・ 数値演算をバカにしては、なりません。 110 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 04:37:37 ] www7b.biglobe.ne.jp/~river_r/alcan/index.html 多桁演算、やってみてます。 111 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 04:52:02 ] >>110 実は C で多桁演算プログラムを書いて、うわ、間違ってる、とか やってた記憶があります。人間って間違えるものなのだな。 今ではそれほど後悔していない。 112 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 04:59:34 ] どうやって間違いを見つけたか。e の計算とかπの計算をやらかすのさ。 １００万桁だって同じはずダロ。違うとしたら、たぶんオレが間違い。 113 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 06:47:02 ] シュリニヴァーサ・ラマヌジャン（Srinivasa Aiyangar Ramanujan、1887年12月22日 - 1920年4月26日）はインドの数学者。 すげ。 114 名前：１３２人目の素数さん mailto:sage [2009/03/11(水) 08:29:26 ] fpprintprec周りの話が、MLでも出てましたね。 あと、 maxima.sourceforge.net/docs/manual/en/maxima_10.html には 　fpprintprec is the number of digits to print when printing an ordinary 　float or bigfloat number. って載ってますね。 現状ではfloatに対する挙動はmanual通りではないけど。 115 名前：１３２人目の素数さん mailto:sage [2009/03/25(水) 01:19:15 ] ev(e, x:1/x) みたいなのを関数化しようとして f(e, x):=ev(e, x:1/x); と定義してもうまくいきません。どうすればいいのでしょう？ 116 名前：盛り上げたい mailto:sage [2009/03/25(水) 08:26:34 ] :=では引数は評価されないので defineを使えばいけるんじゃね？ f(x):=x+1; g(x):=f(x); define(h(x),g(x)); ev(e,x:1/x)の動作例がいまいちわからんが define(f(e,x),ev(e,x:1/x)) てな感じのイメージ 117 名前：115 mailto:sage [2009/03/25(水) 20:32:45 ] >>116 f(x)→f(1/x)みたいな変換をする為に、 ev(f(x), x:1/x)でうまくいきそうなのですが、これを関数化するとこで つまずいてます。 f(x)として与える式の中に、xが全て陽(explicit)に現れてるのなら、 多分、ev()を使わずとも、subst()でx→1/xと置換すればいいのだと思うのですが (ちょっと自信が無いです)。 118 名前：盛り上げたい mailto:sage [2009/03/25(水) 21:17:41 ] それならお望みのものはこれだな reciprocal(func) := func(1/x); 例えばこう使ってくれ f(x):=x$g(a):=1+a+a^2$h(x):=sin(x)$ define(rf(x),reciprocal(f)); define(rg(x),reciprocal(g)); ひとまとめにこうしてしまうのも可 map(lambda([i],reciprocal(i)),[f,g,h]); 高階関数で調べると幸せになれる 119 名前：１３２人目の素数さん [2009/04/24(金) 06:01:37 ] ほっしゅ 120 名前：１３２人目の素数さん [2009/04/24(金) 11:34:46 ] 対数正規分布関数(変数t,パラメータμ,σ)を区間[0,t]で積分して得られる関数に対して実験データを非線形フィットさせてパラメータを決定したいのですが、 手持ちのソフトでは難しく、うまくいきません。 アルゴリズムとしてはlevenberg-marquardt法を用いようと思っているのですが、 maximaではこういったことは可能でしょうか？ 質問わかりづらかったら申し訳ありません。 121 名前：１３２人目の素数さん [2009/05/09(土) 17:42:46 ] Windows7で動かない。 orz 勇んでVistaから乗り換えた俺涙目wwww 122 名前：１３２人目の素数さん [2009/05/11(月) 17:00:39 ] インストールしたのは5.18.2なんですが、開いたwxMaximaには0.8.2と表示されます。 この、0.8.2って何の意味なんですか？ Maxima初心者なので教えてください。よろしくお願いします。 123 名前：１３２人目の素数さん mailto:sage [2009/05/11(月) 17:38:28 ] maximaが本体でwxMaximaはGUI 124 名前：１３２人目の素数さん mailto:sage [2009/05/16(土) 20:44:46 ] 保守 125 名前：１３２人目の素数さん mailto:sage [2009/05/19(火) 08:23:55 ] maximaでは演算子はどう表現するのでしょうか？ たとえば、量子力学で頻繁に使う昇降演算子l_+はつぎの性質を持ちます。 l_+ Y[l,m] = [(l-m)*(l+m+1)]^0.5 * Y[l,m+1] ただし、Y[l,m]は球面調和関数 昇降演算子に限らず、自分で演算子を定義し、使用する方法が知りたいです 126 名前：１３２人目の素数さん mailto:sage [2009/05/19(火) 08:39:44 ] powを演算子化したサンプル infix("pow",100,120)$(a pow b):=a^b$a pow b pow c; 100とか120ってのは左と右への束縛力(結合性みたいなもの)を表す(maxは200) 単項演算子ならinfixではなくprefixやpostfixを使う あまり賢くないので元の式の意味がわからんが適当に改良してみてくれ 127 名前：１３２人目の素数さん mailto:sage [2009/05/19(火) 20:02:57 ] ありがとうございます。 128 名前：１３２人目の素数さん mailto:sage [2009/05/23(土) 23:22:00 ] 関数を返す関数はどう記述されますか？ g(x)を受け取ると、g(x-1)を返すf(x) f(g(x)) := g(x-1); のような感じです。 もしかすると、lispを勉強しないと駄目ですかね 129 名前：１３２人目の素数さん mailto:sage [2009/05/24(日) 09:16:51 ] 自己解決しました。 シミュ板の過去スレのまとめが参考になりました。 ryoma3.hp.infoseek.co.jp/maxima/GO_MAXIMA_coments.pdf f(x,_h) := block([_a],_a:funmake(_h,[t]), subst(x-1,t,ev(_a))); g(x) := 2*x f(x,g); <- 2(x-1) 130 名前：１３２人目の素数さん mailto:sage [2009/05/29(金) 12:32:21 ] ベクトルの外積をするために load(vect); してから cross(v1,v2) としても外積がでないのですが 何か間違っているのでしょうか 131 名前：１３２人目の素数さん mailto:sage [2009/05/29(金) 18:11:31 ] matrixをlistに変換する関数はないのでしょうか listに対するtransposeはないのでしょうか 132 名前：１３２人目の素数さん mailto:sage [2009/05/30(土) 14:49:07 ] >>130 自分で定義する cross( v1, v2 ):= [ v1[2]*v2[3] - v1[3]*v2[2], v1[3]*v2[1] - v1[1]*v2[3], v1[1]*v2[2] - v1[2]*v2[1] ]; 133 名前：１３２人目の素数さん mailto:sage [2009/05/31(日) 12:27:13 ] maximaのfactorの出力結果ってアクセス出来ないのでしょうか 例えばfactor(84);%[1]などとやってもうまくいきません メビウス関数を定義したくて 「平方数で割り切れるかの判定」を、 素因数分解の結果を使って行おうと思ったのですが・・ それか、他になにか良い代替案はありますかね 134 名前：１３２人目の素数さん mailto:sage [2009/06/01(月) 05:17:45 ] load("vect") するとmatrixの掛け算で変なエラーが出ることないないですか？ nxnの行列と1xn　のベクトルの掛け算で 行列のサイズ合わないよ とか 135 名前：クリップボード [2009/06/01(月) 16:29:52 ] >>133 part関数を使ったらどうでしょうか。 　　　　(%i1) factor(84); (%o1) 2^2 3 7 part(%o1,1) 4 part(%o1,2) 3 part(%o1,3) 7 136 名前：クリップボード [2009/06/03(水) 16:15:47 ] >>131 matrixをlistに変換するには、演算子matrixを演算子listで置き換えればよいので substpart関数が使える。 　例えば、 　　　m:matrix([4,3],[2,1])$ substpart("[",m,0); [[4.3],[2,1]] となる。 137 名前：１３２人目の素数さん mailto:sage [2009/06/03(水) 18:03:05 ] なるほど 138 名前：１３２人目の素数さん [2009/06/04(木) 15:00:57 ] ベクトルの次元数を上げたり下げたりはどうやりますか？ (%i7) v0 : [1,2,3]; v1 : [v0,1]; v2 : [1,2,3,4]; v1 . v2; (%o7) [1,2,3] (%o8) [[1,2,3],1] (%o9) [1,2,3,4] incompatible dimensions - cannot multiply -- an error. To debug this try debugmode(true); (%o8)を、[1,2,3,1]にしたいです。 139 名前：１３２人目の素数さん [2009/06/04(木) 15:11:21 ] (%i1) v : [vx,vy,yz]; x : [xx,xy,xz]; v . x; (%o1) [vx,vy,yz] (%o2) [xx,xy,xz] (%o3) xz*yz+vy*xy+vx*xx (%o3)を、「v . x」に変形させることは出来ますか？ 140 名前：クリップボード [2009/06/04(木) 16:35:46 ] >>138 v0:[1,2,3]から[1,2,3,1]をつくるには、endcons関数で可能。 　　　endcons(1,v0); [1,2,3,1] 　　因みに、先頭に要素を入れるには、cons関数を利用する。 　　　cons(1,v0); [1,1,2,3] また、(%o8) [[1,2,3],1]から[1,2,3,1]をつくるには、flatten関数を 　　利用する。 　　　flatten(%o8); [1,2,3,1] 141 名前：１３２人目の素数さん mailto:sage [2009/06/04(木) 17:26:01 ] reshapeを標準関数にいれてほしい 142 名前：変態ロリーマン ◆ghclfYsc82 mailto:sage [2009/06/05(金) 05:52:39 ] 元々はメモリ管理がアカンのやけど、そやけど関数を増やせば増やす程、 更に動作がトロくなりますな 143 名前：クリップボード [2009/06/05(金) 15:34:36 ] >>141 求める関数と同じかわかりませんが,areshape関数を作ってみました。 　　reshape(A,m,n):=block([nr,nl,lis,a,b,c,s], nr:length(A), a:transpose(A), nl:length(a), if m*n # nr*nl then (print("error ",m,"*",n,"should be equal to the number of matrix elements"),return(done)), s:[], b:flatten(substpart("[",a,0)), for i thru n do(c:rest(b,-m*(n-i)), s:endcons(c,s), b:rest(b,m) ), transpose(apply(matrix,s)))$ 　　　 　　実際に適用してみると、 　　　　　p:matrix([1,4,7,10],[2,5,8,11],[3,6,9,12]); [ 1 4 7 10 ] [ 2 5 8 11 ] [ 3 6 9 12 ] reshape(p,2,6); [ 1 3 5 7 9 11 ] [ 2 4 6 8 10 12 ] となる。 144 名前：１３２人目の素数さん mailto:sage [2009/06/06(土) 23:22:07 ] >>140 ありがとうございました。 おかげで出来ました。 145 名前：１３２人目の素数さん mailto:sage [2009/06/11(木) 20:59:04 ] これバグかな？下の式は正常に動かないけども sublist([1,2,3],evenp); sublist([1,2,3],lambda([x],mod(x,2)=0)); 146 名前：１３２人目の素数さん mailto:sage [2009/06/11(木) 21:18:35 ] 自己解決。こうかけば良いようだ sublist([1,2,3],lambda([x],is(mod(x,2)=0))); 147 名前：１３２人目の素数さん mailto:sage [2009/06/16(火) 10:52:37 ] MaximaでZpやF2上での計算って出来ますか？ 148 名前：１３２人目の素数さん mailto:sage [2009/06/16(火) 23:33:33 ] >>147 なかったら自分で作るとか。 149 名前：１３２人目の素数さん mailto:sage [2009/06/17(水) 20:19:26 ] lisp ,@ に相当する演算子はないんだね 150 名前：１３２人目の素数さん mailto:sage [2009/06/17(水) 22:47:35 ] 最小インデックス計算したいだけなのに，なんでうまくいかないんだろう min_index(lst):=block(local(k,mi),mi:1, for k thru length(lst) do( disp([k,mi,lst[k]]), if lst[mi] < lst[k] then mi:k ), disp([mi,lst[mi]]), return(mi)); mir:min_index([3,1,5,0,10,60]) ４が答えなのに６になる 151 名前：１３２人目の素数さん mailto:sage [2009/06/17(水) 23:11:24 ] 全然デバッグはしてない mi(lst):=block([mini:1,minv:lst[1]], for i:1 thru length(lst) do if lst[i]<minv then block(minv:lst[i],mini:i), return(mini)); 152 名前：１３２人目の素数さん mailto:sage [2009/06/23(火) 05:27:36 ] drowでファイル名指定してepsファイルに出力する方法がわからない 153 名前：１３２人目の素数さん mailto:sage [2009/06/24(水) 01:14:07 ] どなたかご存じの方がいらっしゃれば教えていただけないでしょうか。 u(x,y,z,t) = u0(x,y,z,t) + \varepsilon u1(x,y,z,t) (u,u0,u1はすべてベクトル) p(x,y,z,t) = p0(x,y,z,t) + \varepsilon p1(z,y,z,t) (p,p0,p1はすべてスカラー) などと摂動展開を行い、これを基礎方程式に入れて線形的な関係式を求めたいと考えています。 #具体的には連続の式、Navier-Stokes方程式、Poissonの関係式から #音波の波動方程式をMaximaで導出してみたいと考えています。 任意の座標系でできればいいのですが、まずはデカルト座標だけでもできればありがたいです。 以下デカルト座標系であることを前提にして書かせていただきます。 ここで皆さんに伺いたいことは以下の3点です。 1. 摂動パラメータ\varepsilonのようにx,y,z,tに依存しない定数を指定するコマンドを教えていただけないでしょうか。 具体的な数値(1とか2.5など)ではなく、「なんらかの定数値をとる\varepsilon」というように指定できるコマンドがあるとありがたいのですが... 2. VECTパッケージではある変数pがx,y,zに依存することを前提にしているようですが、pがx,y,z以外の変数(たとえばt(時間))にも依存するように指定することは可能でしょうか。 3. 摂動パラメータ\varepsilonについて整理する(MathematicaでいうとCollectコマンドになると思います)コマンドがあるのでしょうか。 MaximaでVECTをロードするとベクトル解析関係の計算を(成分に分解しないで)シンボリックに計算できるようなので 見通しよく計算できると考えています。 MathematicaにせよReduceにせよ、成分に分解しないとベクトル解析の計算ができないので大変見通しが悪くなってしまうのが残念です。 よろしくお願いします。 154 名前：１３２人目の素数さん mailto:sage [2009/06/27(土) 22:32:23 ] maximaで区分関数（？）を表すことってできないのでしょうか？ 例えば、下のような連続の式をグラフで表示させたいのですが。 ・y(x):=2*x　（x<0） ・y(x):=3*x　（x>=0） 155 名前：１３２人目の素数さん mailto:sage [2009/07/01(水) 22:02:03 ] ものすごい基本的なことなんですけれども、 行列の１つの要素を抜き出した時に１×１の行列扱いになって困っているのですが、単なる数字として扱う方法はありませんか？ 例えば A1: ([a11, a12, a13], [a21, a22, a23],[a31, a32, a33]); A2: A[1][1]*A1; みたいな事をやりたいんですけれど、どうやらA[1][1]の部分が行列扱いになっているようでエラーが出ます。 数字でなく式で扱っているのが悪いのでしょうか？出来れば式のままで整理したいのですが・・・ 156 名前：１３２人目の素数さん mailto:sage [2009/07/02(木) 18:59:25 ] それはA1がリストだからじゃないか？ Maximaではリストと行列は全くの別物 A1 : matrix ([a11, a12, a13], [a21, a22, a23],[a31, a32, a33]); 157 名前：１３２人目の素数さん mailto:sage [2009/07/02(木) 23:25:44 ] すいません、書き込みの際にmatrixは書き忘れました。 実際には>>156に書いてあるように宣言していますが、１要素を抜き取ってもそれは数でなく行列扱いになっています。 158 名前：１３２人目の素数さん mailto:sage [2009/07/03(金) 04:46:09 ] >A2: A[1][1]*A1; Aには何が入ってるの？ 159 名前：１３２人目の素数さん mailto:sage [2009/07/03(金) 04:56:11 ] A1の間違いですね、申し訳ないです。 具体的な例を書くつもりが逆に混乱させているようですいません。 とにかく、行列を設定したあとに要素を抽出して、その要素を行列やリストでないただの数字として扱いたいのです。 160 名前：１３２人目の素数さん mailto:sage [2009/07/03(金) 07:51:30 ] >A2: A1[1][1]*A1; この書き方だと、行列A1のすべての要素にA1[1][1](=a11)を掛けた行列がA2に代入される。 >とにかく、行列を設定したあとに要素を抽出して、その要素を行列やリストでないただの数字として扱いたいのです。 A1[1][1],A[1][2]等でアクセスできる。というか、上の例ではしてるのでは？ あと、具体的なインプット、アウトプットを書いてくれると分かりやすい。 161 名前：１３２人目の素数さん mailto:sage [2009/07/12(日) 17:34:19 ] すいません、質問です。 文字が入った行列の行列式を計算しました。 計算自体はうまくいったのですが、途中計算を出すためにはどうすればよいですか？ だれか教えてください。 162 名前：１３２人目の素数さん [2009/07/14(火) 22:24:45 ] A:matrix([a11,a12,....],[a21,a22,....],....) っていちいち書くのが面倒なんですが、自動で設定することはできないんでしょうか？ 163 名前：１３２人目の素数さん [2009/07/16(木) 00:31:40 ] macファイルをbatchコマンドで読み込むと$が無視されて数式がジャンジャン表示されるんだけど、$を有効にするにはどうしたらいい？ 164 名前：１３２人目の素数さん mailto:sage [2009/07/17(金) 05:03:17 ] 初歩的な質問で申し訳ないですけど、 例えば、 Y = X * Y + A っという方程式を、 Y = A / (1-X) のように 左辺にY でまとめて表示させるにはどうしたら良いでしょうか？ できれば、右辺も、X と定数で分けたいです。 複雑な式の場合に、Y = ○X ＋ △ ってまとめたいのですがそんな事が可能？ 165 名前：１３２人目の素数さん [2009/07/18(土) 13:46:33 ] ∬(D)f(x,y)dxdy D:x^2+y^2<=1 といった重積分をするとき integrate(integrate(f(x),y,sqrt(1-(x^2)),-1*sqrt(1-(x^2))),x,1,-1) と累次化した形で入力しているのですが計算してくれません この入力に何か間違いがありますか？ それともやり方自体が間違っているのでしょうか？ 166 名前：165 [2009/07/18(土) 14:32:57 ] すいません 積分範囲が逆でした integrate(integrate(f(x),y,-1*sqrt(1-(x^2)),sqrt(1-(x^2))),x,-1,1) と入力しましたが これでもできませんorz f(x)は任意の式です 167 名前：１３２人目の素数さん mailto:sage [2009/07/18(土) 14:35:31 ] y? 168 名前：165 [2009/07/18(土) 14:48:00 ] すいませんf(x,y)です integrate(integrate(f(x,y),y,-1*sqrt(1-(x^2)),sqrt(1-(x^2))),x,-1,1) 実行すると Is (x-1)*(x+1) positive, negative, or zero?と聞かれて positiveと入力すると Defint: Lower limit of integration must be real. -- an error. To debug this try debugmode(true); って返されるんですよね f(x,y)は１でやってます 169 名前：１３２人目の素数さん mailto:sage [2009/07/18(土) 15:11:57 ] integrateの引数3番目に数式指定するとどうなるの？ 170 名前：１３２人目の素数さん mailto:sage [2009/07/18(土) 15:18:43 ] (x^2-1)が正だとルートの中が負になるので、 Maximaは積分範囲は実数にしてくださいと言っている。 negativeと答えるべき。 171 名前：１３２人目の素数さん mailto:sage [2009/07/23(木) 12:31:05 ] maximaのlistとlispのlistを変換する関数ないのでしょうか 172 名前：１３２人目の素数さん [2009/07/24(金) 00:04:26 ] 10x10正方行列のdeterminantを行列成分で表示させようとしたらメモリが足りん的なこと言われたんだけどどうしたらよいの？ 173 名前：１３２人目の素数さん mailto:sage [2009/07/24(金) 00:11:40 ] あきらめる 174 名前：１３２人目の素数さん mailto:sage [2009/07/24(金) 12:07:49 ] このスレは、明らかに回答者の数が不足している 175 名前：１３２人目の素数さん mailto:sage [2009/07/24(金) 12:47:17 ] オレも質問者の側だ。回答者に成れるくらいレベルアップしたいが。 ぽんぽこさんの本以外でマトモな資料ってある？ ブルーバックスのはオレにはクソだった。 176 名前：１３２人目の素数さん mailto:sage [2009/07/28(火) 05:51:56 ] ブルーバックスのはアマゾンでも評判が悪い 177 名前：１３２人目の素数さん mailto:sage [2009/07/29(水) 16:12:14 ] 次のバージョンはいつ出 translateのバグいつ治るんだろ 178 名前：１３２人目の素数さん mailto:sage [2009/08/11(火) 16:37:44 ] 一番バグが少ないのはどのバージョン？ >>177 みたいなバグがないなら古いバージョンでも使ってみたい 179 名前：１３２人目の素数さん [2009/08/12(水) 23:52:33 ] >>164 solve(y=x*y+a,y); 180 名前：１３２人目の素数さん mailto:sage [2009/08/27(木) 00:16:32 ] >>179 遅くなりましたが、ありがとうございます。かなり助かりました。 181 名前：１３２人目の素数さん mailto:sage [2009/08/28(金) 15:39:57 ] 本家がmediawiki置かせてくれるサーバーどっかにないか探してるね 182 名前：１３２人目の素数さん [2009/09/14(月) 02:48:02 ] a x + b などと掛け算を省略できないのは仕様ですか？ それと構造化言語しか使ったことないんですが、block([x,y], a,b,c)をつかうのでなくlisp風で書く（例えば階乗などを）とどうなるんでしょうか？ 183 名前：１３２人目の素数さん mailto:sage [2009/09/14(月) 23:07:28 ] 1234[100桁]5678 のように省略されているものを、すべて表示するためにはどうすればいいですか？ 184 名前：１３２人目の素数さん mailto:sage [2009/09/15(火) 01:03:44 ] >>183 そういう表示をする実例（計算式）を出してくれ。 185 名前：クリップボード [2009/09/15(火) 12:22:10 ] >>183 set_display(ascii)ですべて表示できます。 元にはset_display(xml)で戻ります。 186 名前：183 mailto:sage [2009/09/15(火) 21:51:21 ] ありがとうございます、助かりました 187 名前：１３２人目の素数さん mailto:sage [2009/09/23(水) 11:11:08 ] Σ[x=1,∞]x(2x-1)/3^x の計算をさせたいのですが、どう表記すればいいでしょうか。 sum(i(2*i-1)/3^i, i, 1, inf); だと駄目っぽいんですよね…。 お願いします 188 名前：１３２人目の素数さん mailto:sage [2009/09/23(水) 13:06:20 ] nusum(i*(2*i-1)/3^i, i, 1, n); limit(%, n, inf); 189 名前：１３２人目の素数さん mailto:sage [2009/09/23(水) 13:56:16 ] >>188 >>187 ではありませんが、勉強になりました。有り難うございます。 190 名前：１３２人目の素数さん mailto:sage [2009/09/23(水) 13:56:31 ] n : 30 でセットですけど、 n = 30 だとどういう意味なんですか？ 結果はブーリアンとして扱われるんですか？ 191 名前：１３２人目の素数さん mailto:sage [2009/09/24(木) 00:49:56 ] >>188 ありがとうございます。 直接、infは駄目なんですね。 192 名前：187,191 mailto:sage [2009/09/24(木) 01:03:49 ] そもそもの表記が間違ってましたね、気づきませんでした。 i(2*i-1) ではなく、 i*(2*i-1) とするべきだったんですね。 193 名前：１３２人目の素数さん mailto:sage [2009/09/25(金) 11:21:16 ] 微分演算子のみを書いておくことって出来ますか？ヤコビアンの確かめで 使いたいんですが、ある式を、たとえば sin(t)*d/dr+1/r*d/dt とかって感じで表現できますか？ 194 名前：１３２人目の素数さん mailto:sage [2009/09/27(日) 22:31:46 ] emacsでmaximaモードってないんでしょうか？ imaximaのことじゃなく、maximaの文法に沿ってインデントしてくれるモードです。 195 名前：１３２人目の素数さん mailto:sage [2009/09/28(月) 22:24:16 ] お前が書けばある 196 名前：１３２人目の素数さん [2009/09/29(火) 00:59:03 ] linsolve([eq1,eq2,eq3],[x,y,z,]) とかってした場合の x = a*b/y みたいな解が出てくるのですが、 変数リストの文字をすべて定数に置き換えて表示させる方法ってありませんか？ 上記の例だと簡単ですが、もっと式、変数の数が増えてくると 解けてるのかどうかの確認すらできません。 よろしくお願いします。 197 名前：１３２人目の素数さん mailto:sage [2009/10/05(月) 15:02:05 ] wxmaximaのどこが便利なのかわからない 198 名前：１３２人目の素数さん mailto:sage [2009/10/27(火) 12:19:24 ] octaveを使っている人はいるのかな。スレッドがない。化学板にもない。 199 名前：１３２人目の素数さん mailto:sage [2009/10/27(火) 13:32:23 ] ⊃ ttp://pc12.2ch.net/test/read.cgi/unix/1086016412/ 200 名前：１３２人目の素数さん mailto:sage [2009/11/20(金) 02:12:57 ] wxMaximaを使ってある関数g(x)がある区間x=x1〜x2でとる値を書き出す方法を教えてください。 xとg(x)をデータリスト(dat)に格納して load(numericalio); write_data(dat,"hoge.dat"); で書き出すことは分かったのですが格納の仕方がわかりません。 参考ページ ttp://cosmo.phys.hirosaki-u.ac.jp/wiki.cgi/maxima?page=%BE%EF%C8%F9%CA%AC%CA%FD%C4%F8%BC%B0%A4%CE%BF%F4%C3%CD%B2%F2%CB%A1 dat:rk( x, x, 1, [t, 0, 2, 0.1] );の部分です。 宜しくお願いします。 201 名前：115 mailto:sage [2009/11/20(金) 22:44:39 ] 結局at()を使って r(e,x):=at(e,x=1/x); みたいにやる事になりそうです。 202 名前：１３２人目の素数さん mailto:sage [2009/12/02(水) 22:44:30 ] wxMaximaのテキスト入力等に日本語を使用する事はできますか？ xml出力でテキスト・タイトル等が残せるのに気がついて、 更に日本語が使えると非常に便利だと思うんですが… 203 名前：１３２人目の素数さん mailto:sage [2009/12/06(日) 19:10:26 ] 行列 4*4の行列式を計算して算出しました。 このときに変数を用いたのですが、それに変数に値を代入して計算を行うことができますか？ 204 名前：１３２人目の素数さん mailto:sage [2009/12/06(日) 22:23:11 ] 質問を重ねてすみませんが、 コマンドプロンプトの方は起動すると自動的に (%i1)　とでて、式を入力し、Enterを押すと計算が行われるのですが、 GUI操作の方であるwxMaxima 0.8.3を起動して文字を入力しても (%i1)が出てきません。 自分で　 (%i1)　計算式　Enter を入力しても駄目です。 GUI上ではどのようにして計算式を入力すればいいのでしょうか？ 205 名前：１３２人目の素数さん mailto:sage [2009/12/06(日) 22:33:32 ] 数式入力後　Shift+Enterで実行だったのですね。 スレ汚しすみません。 206 名前：１３２人目の素数さん mailto:sage [2009/12/07(月) 18:14:50 ] axiom ってどうなんだろう maximaから乗り換えたとかいう書き込みも時々みる だけどmaximaのlispで関数書ける便利さは捨てれないなあ 207 名前：１３２人目の素数さん mailto:sage [2009/12/12(土) 10:28:20 ] >>206 maxima言語(？)のみで完結せずに、lispに頼ってるとこは、 あんまり美しくないなぁとか思ってるんですけど…駄目？ 208 名前：１３２人目の素数さん mailto:sage [2009/12/14(月) 10:33:45 ] ?princ("test"); のように先頭に ? をつけるると lisp の関数が呼び出せるようですが、xx-yy のように - を含む関数の場合、呼び出す方法は あるのでしょうか？ そのままかくと - が演算子として認識されてしまいます。 209 名前：１３２人目の素数さん [2009/12/16(水) 07:36:13 ] 5.20.1あげ 210 名前：１３２人目の素数さん [2009/12/17(木) 20:27:10 ] 2次元実数上で、領域の包含関係を判定するソフトを作りたいのですが とりあえず、一次不等式と円の内側を表す不等式　　(x-a)^2 + (y-b)^b <= rなど） を処理できるものを作ろうと思ってまつ。 Maximaでこれらの不等式制約の処理（交点を求めて利用するとか）できますか？ おねがいします。 211 名前：１３２人目の素数さん mailto:sage [2009/12/23(水) 11:58:26 ] >>208 バックスラッシュでエスケープするそうです。 ttp://maxima.sourceforge.net/docs/manual/en/maxima_3.html >Hyphen -, asterisk *, or other special characters in Lisp symbols >must be escaped by backslash \ where they appear in Maxima code. >For example, the Lisp identifier *foo-bar* is written ?\*foo\-bar\* in >Maxima. 212 名前：１３２人目の素数さん mailto:sage [2009/12/25(金) 01:20:53 ] >>211 うまくいきました。ありがとうございます。 マニュアルに載っていたんですね。 213 名前：１３２人目の素数さん [2010/01/03(日) 03:07:15 ] 外部にファイルを書き出させたいのですが、（tex()の結果のみを外部ファイルに出力） どのように書けばよいのでしょうか。 検索してみてもsave()が出てくるだけで，直接ファイルに書き出す方法を探しています。 214 名前：１３２人目の素数さん mailto:sage [2010/01/04(月) 21:20:19 ] >>213 with_stdout (f, expr_1, expr_2, expr_3, ...) は、どうでしょうか。 215 名前：１３２人目の素数さん mailto:sage [2010/01/04(月) 21:35:03 ] >>213 tex(exp(1), "foobaa.txt"); とか。 ttp://maxima.sourceforge.net/docs/manual/en/maxima_9.html >destination may be an output stream or file name. When destination is >a file name, tex appends its output to the file. The functions openw >and opena create output streams. 216 名前：１３２人目の素数さん mailto:sage [2010/01/08(金) 23:26:43 ] >>210 これはできそうでできない。というか、LISPで最初から書けばできるといういみではできるがな 217 名前：猫は淫獣 ◆ghclfYsc82 mailto:age [2010/01/18(月) 00:56:50 ] ココでちょっとしたメッセージや ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ ★★★「小沢氏は主張を通して検察に対して徹底的に対抗すべし」★★★ 小沢先生、頑張って下さい。私は最後まで味方になります。 猫 218 名前：１３２人目の素数さん [2010/01/25(月) 22:09:05 ] すみません、複雑な式をfloat表示させると、 ( 3.14133*(10^25) * x^2 + 2.12412*(10^24) * x^3 + 2.12412*(10^23) * x )/ ( 6.124123*(10^23) * x^2 + 5.124612*(10^25) * x^3 + 4.124312*(10^24) * x ) といった風に、分子分母で明らかに通分出来そうな 10の指数が現れて、とても困っています。 これらを通分するにはどうすればよいのでしょうか？ どうかよろしくお願いします。 219 名前：クリップボード [2010/01/30(土) 11:34:42 ] >>218 分子分母を通分するには,cancellingの機能を有するrat,ratsimp,factor等の関数を 利用すればよい。 220 名前：１３２人目の素数さん [2010/01/31(日) 17:00:33 ] もう３０年以上も経つ古い設計のままのシステムだから、 完全にリフォームしなおすか、さもなければ安楽死させてはどうだろうか？ 221 名前：１３２人目の素数さん mailto:sage [2010/02/01(月) 15:14:59 ] 数値計算関係をgsllあたりに統合してほしい 222 名前：１３２人目の素数さん mailto:sage [2010/02/18(木) 17:59:58 ] 今、俺に一番必要なのがこのソフトだと悟った 俺の逝くべき道が見えてきた これから質問することが多々あると思うので どうか宜しくお願いたもう 223 名前：１３２人目の素数さん [2010/02/18(木) 19:56:30 ] 質問です shikou: 10000$ mieta: 0$ for i: 1 thru shikou do( x: random(9.0), y: random(9.0), if abs(x - y) < 1.0 then mieta: mieta + 1 )$ kakuritu: float(mieta / shikou)$ print("simulation=", kakuritu)$ print("theory=", float(17/81))$ とwxMaxima 0.8.3に打ち込んだんですけど どうやってシミュレーションを始めたらいいんですか？ Start Animationとか書いてある三角形のマークはグレイになってて押せません 224 名前：223 mailto:sage [2010/02/18(木) 20:40:29 ] うわ、自己解決しました wxMaximaの方じゃなくて xmaximaの方にコピペするだけかよ・・・ てか、じゃあ、wxMaximaはなんのためにあるんだよ・・・

---

---

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

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