07/10/11 20:12:37
>>808
octaveだと数値計算であって、式では解いてくれないじゃないかー
810:名無しさん＠５周年
07/10/12 05:57:17
>>809
確かにそうだね。>>805みたいに式を解くにはmaxima、mathematicaだね
811:名無しさん＠５周年
07/10/24 05:34:04
5.13.0 for Win です。
func(t) := (sin(t)^4*z^4*R^4+2*sin(t)^2*z^4*H^2*R^2+z^4*H^4)
/(4*sin(t)^2*H^4) -z^4/(4*sin(t)^2)
を
integrate( func(t), t, 0, 2*%pi )
とすると 0 になってしまいますが、
一旦 factor して
integrate( factor( func(t) ), t, 0, 2*%pi )
だと正しく
z^4*R^2*(%pi*R^2+4*%pi*H^2)/(4*H^4)
と出ます。
何が悪いのでしょうか？
812:名無しさん＠５周年
07/10/24 22:45:50
>>811
その関数のfunc(0)の値がアレだからじゃないすか？
813:811
07/10/25 00:55:54
>>812
どもです。
factor せずに不定積分すると
z^4/(4*tan(t))
という項が入るため t = 0 に出来ないようですね。
factor すると上記の項がまとまるので大丈夫と。
814:名無しさん＠５周年
07/11/02 01:59:02
教えて下さい。多項式近似の数値係数を求めようとしています。
このようなヤツを実行させると
URLﾘﾝｸ(aikofan.dee.cc)
・・結果が得られますが、重要かつ必要な部分
URLﾘﾝｸ(aikofan.dee.cc)
だけをファイルに書き出したいと思います。どのようにしたらいいのでしょうか。
今はコピペと手動編集をしています。バカっぽいですね。
説明書は読んだのですが、なかなか目的に合う関数・コマンドが見つかりません。
よろしくお願い致します。
815:名無しさん＠５周年
07/11/02 03:13:48
ファイルを操作する関数が判らないってこと？
君のプログラミングスキルが判らんから、どこから説明したらいいやらなのだが。
816:名無しさん＠５周年
07/11/02 04:10:47
>>814
エディタは何を使ってるの？
817:814
07/11/02 09:07:48
>>816
普通のテキストエディターです。
>>815
＞君のプログラミングスキルが判らんから、どこから説明したらいいやら
こんな C のコードで MAXIMA 用のプログラム (コマンド) を作っています。
fprintf(ofp, "display2d:false $ ratprint:false $ fpprec:8 $\r\n");
fprintf(ofp, "fa(x) := a5*x^5 + a4*x^4 + a3*x^3 + a2*x^2 + a1*x + a0 $\r\n");
int temp = TEMP_INIT;
for (i = 0; i < MAXI; i += 1) {
fprintf(ofp, "print(\"/* temp = %d */\") $\r\n", temp);
fprintf(ofp, "solve([\r\n");
for (j = 0; j < MAXJ; j += 1)
fprintf(ofp, "fa(%12.8g)=log(%g)%s\r\n", d[i][j], r[j], (j == MAXJ - 1) ? "" : ",");
fprintf(ofp, "],[a5,a4,a3,a2,a1,a0]) $\r\n");
fprintf(ofp, "bfloat(%%);\r\n\r\n");
temp += TEMP_STEP;
}
ratprint:false は効果ありなのに、keepfloat:true は無効なんですね。よろしくお願いします。
818:名無しさん＠５周年
07/11/02 22:28:16
display2d:false $ fpprec:8 $
fb(x) := b6*x^6 + b5*x^5 + b4*x^4 + b3*x^3 + b2*x^2 + b1*x + b0 $
ss:solve([ fb(-5)= -7260.3313,
fb(5)= -5974.0679,
fb(15)= -4943.6644,
fb(25)= -4136.0326,
fb(35)= -3513.7263,
fb(45)= -2993.4723,
fb(55)= -2566.961 ],[b6,b5,b4,b3,b2,b1,b0]) $
b1:bfloat(ss[1][1]);
ってな事じゃないよね？
819:名無しさん＠５周年
07/11/02 22:29:51
書き間違えた
display2d:false $ fpprec:8 $
fb(x) := b6*x^6 + b5*x^5 + b4*x^4 + b3*x^3 + b2*x^2 + b1*x + b0 $
ss:solve([ fb(-5)= -7260.3313,
fb(5)= -5974.0679,
fb(15)= -4943.6644,
fb(25)= -4136.0326,
fb(35)= -3513.7263,
fb(45)= -2993.4723,
fb(55)= -2566.961 ],[b6,b5,b4,b3,b2,b1,b0]) $
a1:bfloat(rhs(ss[1][1]));
820:814, 817
07/11/03 07:36:28
>>819
ありがとうございます。
display2d:false $ ratprint:false $ fpprec:8 $
fb(x) := b6*x^6 + b5*x^5 + b4*x^4 + b3*x^3 + b2*x^2 + b1*x + b0 $
cof:solve([ fb(-5)= -7260.3313, fb(5)= -5974.0679, fb(15)= -4943.6644,
fb(25)= -4136.0326, fb(35)= -3513.7263, fb(45)= -2993.4723,
fb(55)= -2566.961 ],[b6,b5,b4,b3,b2,b1,b0]) $
cof;
for i:1 thru 7 do
print("b", 7-i, "=", bfloat(rhs(cof[1][i]))) $
これの print() 出力をファイルに書き出したい、というのが質問でした。
その他ヒントがありましたら、よろしく。
821:820
07/11/03 07:46:19
あっ、気付きました。これでいいんですね。
for i:1 thru 7 do print(bfloat(cof[1][i])) $
822:820
07/11/03 09:37:38
追記。例えば gnuplot の説明書を見ると、
Examples:
To launch an interactive session:
gnuplot
To launch a batch session using two command files "input1" and "input2":
gnuplot input1 input2
などと書いてあって、shell コマンドとか system() で gnuplot に実行すべき
入力ファイルを指定して起動できます。
MAXIMA の場合も同じようなことができるとは思うのですが、具体的にどうしたら
いいのか、説明が見つかっていません。どこを探したらいいのでしょうか。
よろしくお願いします。
823:822
07/11/03 13:31:04
小出しにせずに全部ぶっちゃけますと、
URLﾘﾝｸ(aikofan.dee.cc)
これの処理：
hum-SC-2.log → [t20.exe] → 1.txt
1.txt → [xmaxima.exe] ＋ コピペ操作 → 2.txt
2.txt → [t21.exe] → 3.txt
3.txt → [xmaxima.exe] ＋ コピペ操作 → 4.txt
4.txt → [t22.exe] → 5.txt
"[ ]" に入っているのはプログラム名。
t2?.exe は MSYS + MinGW 環境で t2?.c から makefile のように作成しました。
MAXIMA 実行のところで "コピペ操作" (手動) をしています。これをなくせれば、
全部 shell とか batch で済むのに・・・。 いろいろヒントを、お教え下さい。
(hum-SC-2.log から一挙に MAXIMA で 5.txt 類似が出来そうな気もするのですが・・)
824:名無しさん＠５周年
07/11/03 17:26:51
c:\test\
というフォルダがあるとして、
stringout("c:\\test\\out.txt",a1)$
とか？
825:823
07/11/03 17:34:49
ありがとう。ここも見てみます。
URLﾘﾝｸ(cosmo.phys.hirosaki-u.ac.jp)
826:名無しさん＠５周年
07/11/03 17:51:04
>>>824 これですか。
URLﾘﾝｸ(www.bekkoame.ne.jp)
ありがとう。一般人にはむづかしすぎる説明文ですが、まずはやってみます。
827:名無しさん＠５周年
07/11/03 22:36:19
一般人ｗｗｗ
828:名無しさん＠５周年
07/11/03 23:26:11
>>827
いやね、ここでいう「一般人」とは MAXIMA に関して、
それほど詳しくない人のこと。
全ての人が MAXIMA 関して、とても詳しい、なんてことは、ないだろ。
↑
変な言い方。一応は反省する。
829:名無しさん＠５周年
07/11/19 11:03:26
830:名無しさん＠５周年
07/11/24 15:43:31
URLﾘﾝｸ(www.interq.or.jp)
はじめてでもできる Maxima のインストール（Windows XP 編）
を見て maxima-5.13.0a.exe をインストールしたんですが
スタートメニューから「Command line Maxima」を実行すると
Error in SYSTEM::SET-UP-TOP-LEVEL [or a callee]: :ROOT is an illegal pathname name.
Fast links are on: do (use-fast-links nil) for debugging
Broken at SYSTEM:STAT. Type :H for Help.
1 (Abort) Return to top level.
dbl:>>
と表示されてしまいます。エラーメッセージの単語をキーワードにして検索してみても
こういう症状が出た人はないようです…… どうすれば動くか分かる人はいませんか？
831:名無しさん＠５周年
07/11/24 21:26:37
>>830
とりあえずは、
wxMaxima
を使っていれば問題ないとは思うが、
illegal pathname ってことは、どうせ何かのパスが漢字なんじゃないかと邪推。
XPにログオンしているusernameに漢字を使っているとか、そんなことはないすか？
832:830
07/11/25 08:10:04
まだ悪戦苦闘してます。
>>831
最初wxMaximaを起動したんですが、「Maxima本体と通信できなくてタイムアウト」ってエラーが出たので、
あれ？と思ったらコマンドラインからも起動できなかったのです。XPだからWin9xの問題でも無いし、本体で
エラーが出るからFirewallの問題でも無いし、DEPの問題でもないし（← boot.iniで無効にしてある）。
デフォルトのインストール先は C:\Program Files\Maxima-5.13.0 ですが、アプリケーションをインストール
するフォルダやユーザー名は昔からすべて（いわゆる）半角で、日本語や空白を含めない主義です。
一度アンインストールして1つ前の旧バージョンをインストールしてみたんですがこれも駄目。
他のアプリケーションの影響を受けるってことがあるのかな？ 環境変数の問題でもなさそうだし……
833:830
07/11/25 08:33:00
batファイルの内容を見て流れを追ってみたけど、特に変なパス（空白入りとか）は設定されてないみたい。
winとunixのパスの違い（￥と／）でもないよなぁ……
エラーメッセージの単語を手がかりにinfoファイルをgrepしたり、リファレンスのchmファイルを探ってみたけど出てこない。
何かそれ以前の問題なのか…… orz
834:名無しさん＠５周年
07/11/25 10:40:03
>>832
>Maxima本体と通信できなくて
ファイアーウォールとかのセキュリティソフトとか、なんとかだったりして？
835:名無しさん＠５周年
07/11/25 12:11:36
>>834
それは最初にブロックするかどうかダイアログが出るのですが、許可してます。
（ファイアウォールソフトは内蔵のWindowsファイアウォールのみ）
それにファイアウォール云々というのはREADMEにも書いてある通り wxMaximaを
使う場合に関係あることなので、本体のコマンドラインのプログラムが単体で
起動できないのはやっぱりおかしいのです。
話は逸れますが、自分はPCの設定はあまりいじらず、できるだけデフォルトの状態で
使うようにしています。常駐ソフトも気を付けてなるべく入れず、最小限に。
デスクトップ改造ソフトなんてもってのほか。何が駄目なのか……
836:名無しさん＠５周年
07/11/27 23:22:13
>>830 です。相変わらずMaximaが起動できない原因はさっぱり分からないのですが、代わりに
URLﾘﾝｸ(www.knoppix-math.org) KNOPPIX/Math Project
という数学系のソフトが詰まった1CD Linuxを見つけたのでこれを試してみようと思います。お騒がせしました。
837:名無しさん＠５周年
07/12/16 18:52:39
wxMaximaもxMaximaも使えません・・。
xMaxima:
Starting maxima timed out. Wait longer?
wxMaxima:
Not connected to maxima!
ためしにコマンドプロンプトからmaxima.batを起動してみると、
C:\Program Files\Maxima-5.13.0\bin>maxima.bat
Unrecoverable error: Segmentation violation..
This application has requested the Runtime to terminate it in an unusual way.
Please contact the application's support team for more information.
何ですかこれは？まったく使い物になりません。
ファイアーウォール(ZoneAlarm)、安置ウイルス（Avast Home Edition）は無罪です。多分。
838:名無しさん＠５周年
07/12/16 19:34:09
5.13.99rc1を試してみるとか
839:名無しさん＠５周年
07/12/16 20:04:30
おなじです
840:名無しさん＠５周年
07/12/16 21:28:29
>>837
まぁキミの環境が運悪くダメなんだろう、としか言えない。
だいたい、普通に使えているヒトの方が圧倒的なので。
841:名無しさん＠５周年
07/12/16 22:06:56
>>837
URLﾘﾝｸ(www.interq.or.jp)
再起動しても立ち上がらない場合は、データ実行防止（DEP）が Maxima をブロックしている可能性がありますので、ブロックを解除してみてください。手順は以下のとおりです。
これはどうでしょう？
842:名無しさん＠５周年
07/12/21 23:27:50
factor(x^4 + y^4);
を実行して因数分解をさせてみたけど、できないみたい。
人間はできるんだけどね。
843:名無しさん＠５周年
07/12/21 23:35:00
>>842
それいったら、x^2+y^2だって出来ないよｗ
844:名無しさん＠５周年
07/12/21 23:50:58
いや、x^4 + y^4 は
(x^2+y^2 - 2^0.5*x*y) (x^2+y^2 + 2^0.5*x*y)
に因数分解できるんだけどね。
845:名無しさん＠５周年
07/12/22 00:20:00
>>844
それって因数分解っていうかなぁ？言うのか
846:名無しさん＠５周年
07/12/22 00:23:36
>>843
それなら x^2+y^2だって、(x+yi)(x-yi) って書いてくれなきゃヤダなヒトですか？
847:名無しさん＠５周年
07/12/24 15:44:02
ラベルの低さにgkbr
848:名無しさん＠５周年
07/12/24 17:37:06
Windows 版
maxima-5.14.0.exe
(2007-12-22 08:43)
URLﾘﾝｸ(sourceforge.net)
出ています。
849:837
07/12/24 22:09:31
動きました。
ZoneAlarmのせいだったのかも・・・
850:名無しさん＠５周年
07/12/24 23:57:00
>>848
どこが変ったか、簡単に詳しく
851:名無しさん＠５周年
07/12/25 22:28:28
>>848
GNU common lisp ベースになってますね！
852:名無しさん＠５周年
07/12/26 00:02:30
Windows版は前からそうじゃなかったっけか？
853:名無しさん＠５周年
07/12/26 21:10:19
既出でしょうが、感動した１行、
フーリエ級数（方形波を正弦波で近似）の描画。
plot2d (sum(sin((2*n-1)*x)/(2*n-1), n, 1, 16), [x, -5, 5])$
gnuplot で苦労していた処理が、
たった１行で・・・嬉しくてつい涙が出てしまった。
以上、Windows 版 wxMaxima にて。
854:名無しさん＠５周年
07/12/27 00:44:26
>>853
やってみた！
nをどんどん増やしていくと方形に近づくのが分かるね
855:名無しさん＠５周年
07/12/27 03:24:42
>>853
ギブス現象（ フーリエ級数の奇妙な振る舞い ）
方形波のように不連続点を持つ周期関数をフーリエ級数で表した場合に生じる特有の現象。
これを簡単に確認できたので、感動したのさ！
856:名無しさん＠５周年
07/12/28 07:49:11
>>851
だったっけ?
Windows版もClispベースじゃなかったっけ?
(GCLはANSI規格準拠じゃなくなったらしいんで、コード互換性の問題で変わったと思う。
Linux版もClispベースだった)
もっとも「ソースからビルド」するんだったら好きな処理系選べると思う。
今回初めてCMUCLでビルドした。
速くなったかどうかは分からんw
857:名無しさん＠５周年
08/01/01 00:29:13
(%i1) solve([x+y=1,x-y=0],[x,y]);
(%o1) [[x=1/2,y=1/2]]
とした後に
(%i2) x+y;
(%o2) 1
となるようにxとyに%i1の結果を代入したいのだけどどうしたらいいですか？
今は当然
(%o2) x+y
となります。
858:名無しさん＠５周年
08/01/01 10:41:06
>>857
(%i2)first(first(%o1))+last(first(%o1));
(%o2) y+x=1
となるはず。
基本、MaximaはLisp系なんで、Lispと同じように扱ったりする。
元々の
(%i1) solve([x+y=1,x-y=0],[x,y]);
(%o1) [[x=1/2,y=1/2]]
の2行めはLispで言うと
'((1/2 1/2))
の様なリストを表している。つまり、Common Lisp(CL)で言うと、
(+ (car (car '(1/2 1/2))) (last (car '(1/2 1/2))))
のような操作をすればよい。
リストのリストから先頭要素を取り出すのがCLだとcarなんだけど、それがMaximaでは
firstとして定義されている。
詳しくは
URLﾘﾝｸ(www.bekkoame.ne.jp)
参照。
859:名無しさん＠５周年
08/01/01 10:42:32
>>858
×リストのリストから先頭要素を取り出すのがCLだとcar
○リストから先頭要素を取り出すのがCLだとcar
メンゴ
860:名無しさん＠５周年
08/01/01 11:41:29
CLにもfirstはあるよ。
861:名無しさん＠５周年
08/01/01 13:48:20
>>860
知ってるけど使わない人多いだろうからw
862:857
08/01/01 16:42:55
>>858
出来ました。ありがとうございました。
863:名無しさん＠５周年
08/01/03 17:50:49
>>857
>>858 の資料 p57 のように ev 関数を使って
(%i1) solve([x+y=1,x-y=0],[x,y]);
(%o1) [[x=1/2,y=1/2]]
(%i2) ev(x+y, %o1[1]);
(%o2) 1
としても良いのでは？もしくは
(%i3) %o1[1][1]+%o1[1][2];
(%o3) y+x=1
864:名無しさん＠５周年
08/01/08 16:07:28
Maximaで行列Mのi乗計算(何かと言うと、確率遷移行列)させて
iがある閾値を超えるiを手っ取り早く求めたいのですが、
二分法みたいな奴って組み込まれてますか？
865:864
08/01/08 17:42:56
結局if then elseとかwhile do とか作ってMaximaでやってみてます。
やってみたら遅すぎて泣けてきましたが(多分行列の1024乗とか
せずに答を出す方法があるんだろうな‥‥‥)
866:名無しさん＠５周年
08/01/09 18:42:38
てす
867:名無しさん＠５周年
08/01/11 09:01:58
m
868:名無しさん＠５周年
08/01/15 12:14:21
-
869:名無しさん＠５周年
08/01/15 12:16:40
-
870:名無しさん＠５周年
08/01/15 12:17:21
-
871:名無しさん＠５周年
08/01/15 12:17:41
-
872:名無しさん＠５周年
08/01/15 12:18:02
-
873:名無しさん＠５周年
08/01/15 12:18:32
-
874:名無しさん＠５周年
08/01/19 12:02:43
tan(y) + e^(-y) = x
のグラフを書くにはどうしたらよいですか？
（陰関数のグラフ？）
875:名無しさん＠５周年
08/01/23 13:16:01
ｔ
876:名無しさん＠５周年
08/03/09 09:59:19
>>874
最近触ってないから忘れたけど、陰関数用のパッケージは付属してたよ。
877:名無しさん＠５周年
08/03/09 10:02:00
>>858
HyperSpecではfirstはcarと同じに定義されてるよ。CLではどちらも使えるよう
になってる。ちょっと古い標準的なlispの教科書にはfirstで統一されてるもの
があるよ。
878:名無しさん＠５周年
08/03/21 14:55:25
cで既に書いてあるプログラムと連携したいのだけど
これ使えばいいのかな？
URLﾘﾝｸ(common-lisp.net)
879:名無しさん＠５周年
08/03/25 12:36:04
数値計算８割、たまに気が向いた時だけ数式処理という使い方の予定なのですが
matlabの方がましなのでしょうか？
880:名無しさん＠５周年
08/03/25 13:37:16
両方試してみればいいじゃん
881:名無しさん＠５周年
08/03/26 09:39:10
>>879
SAGEというのがあってな。先日、数学会で話を聞いた。面白そうだった。
でもウチは PPC 版の Mac なんだよ。
882:名無しさん＠５周年
08/03/26 22:57:46
>>879
MAXIMA と OCTAVE と両方使っとけ
883:名無しさん＠５周年
08/03/30 10:29:39
URLﾘﾝｸ(d.hatena.ne.jp)
maximaからcのプログラムを呼ぶ方法
884:名無しさん＠５周年
08/03/30 10:36:22
[x1,x2,x3...]
と添字つき変数を並べたいのですがlispでdefmacro使うのが正しいやりかたなのでしょうか？
pythonにはmacroがないのでsageの場合には同じことはできないということでいいのでしょうか？
885:名無しさん＠５周年
08/04/05 21:26:22
axionのスレはどこかにないのでしょうか？
axionはlispとboot(pythonの前身？)の混合で書かれてるみたいですが
maximaみたいにto_lisp();でlispで関数かけたりするんでしょうか？
886:名無しさん＠５周年
08/04/06 14:56:08
ワトソン君、参考書を網羅的に列挙してみてくれたまえ。
887:名無しさん＠５周年
08/04/07 07:35:33
URLﾘﾝｸ(axiom-wiki.newsynthesis.org)
888:名無しさん＠５周年
08/04/17 07:59:55
allroots()って、fpprecでは精度を変更出来ないのでしょうか？
889:名無しさん＠５周年
08/04/20 19:52:39
Vistaの上でマキシマ動きますか？
890:名無しさん＠５周年
08/04/20 20:07:01
>>889
動かしてみて動かなかったら、動かないんだよ。きっと。
891:名無しさん＠５周年
08/04/20 21:56:47
>>889
データ実行防止を解除しないとダメみたいだよ
892:名無しさん＠５周年
08/04/21 00:58:05
データ実行防止って？
893:名無しさん＠５周年
08/04/21 21:21:34
>>892
URLﾘﾝｸ(windowshelp.microsoft.com)
894:名無しさん＠５周年
08/04/21 21:26:15
↓Windows Vista対応
URLﾘﾝｸ(shop.kodansha.jp)
895:大変さん
08/04/22 17:15:02
日本語化ってないのかな？
896:名無しさん＠５周年
08/04/22 21:29:00
>>895
HELPを？
897:名無しさん＠５周年
08/04/26 15:51:42
>>881
これ？
URLﾘﾝｸ(www.sagemath.org)
Maximaに計算の下請けを出してたりしてるけど、そう云った部分をどこ迄減らして
いくかが楽しみですね。＞さげます
898:881
08/04/26 18:12:38
>>897
> >>881
> これ？
> URLﾘﾝｸ(www.sagemath.org)
それ！
で、sage 起動後に MAXIMA に切り替える事もできる。あと sage から
間接的に MAXIMA を呼ぶときには、注意が必要。例えば
f( x ) = a * x + b
を積分するとき MAXIMA は a, b について何も言わなくてもイイが sage
では何かを言わないとイケナイ。詳しくはチュートリアル読んでください。
GAP などもバンドルしてるので便利だとは思います。
899:名無しさん＠５周年
08/04/26 22:12:33
>>881
さげます？
900:名無しさん＠５周年
08/04/27 01:51:33
いやそれあげてますから
901:名無しさん＠５周年
08/04/29 15:38:25
5.15.0記念あげ
902:名無しさん＠５周年
08/04/30 17:37:50
今の理系の学科では普通に
maxima,spice,matlab,scilabを使っているのですか？
903:名無しさん＠５周年
08/04/30 23:55:17
>>902
matlabはアカデミックライセンスで。
scilabは趣味だな。
学校だったらscilabよりはoctaveを薦めていそうだ。
904:名無しさん＠５周年
08/05/01 06:42:31
>>903
僕の職場では開発設計なんですが、
工学部卒のはずなのに数式に拒絶反応を示すか
計算過程見せても理解してもらえません。
過渡現象すら全く理解してるひとが無く
電子課卒がSPICEすら使えないですよ。
情報工出身の僕ですらspice使えるのに
無力さを感じるんですが
世の中は全般的に数学が嫌いなんですか？
905:名無しさん＠５周年
08/05/01 08:27:22
すれ違いだろどっかいけクズ
906:名無しさん＠５周年
08/05/01 12:09:48
>>905
いや、喪前のようなバカは氏ねばいいんだよ。
氏ね
907:名無しさん＠５周年
08/05/01 14:07:24
>>904
はい。残念ながら世の中は全般的に数学が嫌いです。
908:名無しさん＠５周年
08/05/01 20:32:46
あなたは思いやりのあるうさぎの船に乗りますか？
それとも意地悪たぬきの泥船に乗りますか？
909:名無しさん＠５周年
08/05/02 10:53:15
>>904
昔いた職場がそんな感じだ
理解してくれる上司にだけアピールしとけば給料だけはあがった
910:名無しさん＠５周年
08/05/03 11:12:34
>>879
知人は出来るだけ複数使うようにしていると言ってた。
911:名無しさん＠５周年
08/05/04 17:36:28
◆こんなソフトウェアつくってください～Part10～◆で
>760 名無しさん＠お腹いっぱい。 [sage] Date:2008/02/12(火) 01:43:53 ID:Q+FVK/ZI0 Be:
>MAXIMAの命令、微分diff(f(x),x)、積分integrate(f(x),x)、方程式の解solve()、行列の演算matrix([a,b],[c,d])を実行し、結果をExcelに返すExcelのアドインをつくってください。
>例えばExcel側でセルA2に'diff(sin(A1),A1)と打ち、セルA3を指定すると、MAXIMAで処理され、A3に=cos(A1)が入ってくるというものです。
>TeXmacsがMAXIMAとプロセス間通信をしているそうで、Excelでも可能なのではと考えました。
>具体的には、Excel側入力の先頭の『'』を省き、末尾に「;」をつけ、MAXIMAにバックグランド処理をさせ、アウトプットの先頭に『=』をつけて指定セルに返します。
>A1、R1C1はMAXIMAが関数とみなしますが、$、[、]の入った$A$1、R[-1]C[-1]は工夫が必要です。
てなお願いをして思いっきり無視されたんですが、難しいのですかね？
Excelで微分積分が出来たら便利だと思うのですが・・・。
どなたかお願いできませんか？
912:名無しさん＠５周年
08/05/04 17:48:15
Excelってとこでやる気減退する人が多いと想像する
913:910
08/05/04 18:51:58
>>910
複数使うと時々結果が違ってて面白いとも言ってた。
一種類じゃ危険。
914:名無しさん＠５周年
08/05/04 19:38:23
行列Rをtの関数R(t)にするにはどけんすればよかとですか？
915:名無しさん＠５周年
08/05/05 15:42:13
>>914
全ての成分を t の関数にすればいい？
916:名無しさん＠５周年
08/05/11 02:32:36
64bitについてコメント：
64bit linux上でテストをしてみた。対象は 32bit版のcmucl 19eと 64bit版のsbcl
1.0.16.11 を使ったmaxima 5.15.0
5%ほど64bit版sbclの方が速かったです。
32bit環境がないので32bit環境でのものと比較はできてない。64bitt環境とは体感
速度でも速さの違いは感じるので32bit環境とはもう少しさがあるかもしれない。
利用したもの
time maxima -b maxima-5.15.0/tests/rexamples.mac > /dev/null
cmucl 32bit
4.46user 0.76system 0:05.23elapsed 99%CPU (0avgtext+0avgdata
0maxresident)k
0inputs+0outputs (0major+108275minor)pagefaults 0swaps
sbcl 64bit
4.55user 0.46system 0:05.02elapsed 100%CPU (0avgtext+0avgdata
0maxresident)k
0inputs+0outputs (0major+38911minor)pagefaults 0swaps
時間は一例ですが、何度かの試行をやってもほぼ同じ値になっていた。
環境は Core2Duo E8750 2.66GHz + 4GByte linux 2.6.24 カーネル x8664版
他の環境と比べてみてね。
917:名無しさん＠５周年
08/05/11 02:35:42
32bit版バイナリを動かすためには32bitバイナリが動くライブラリを入れれ
ば使えるようになります。補足です。
918:名無しさん＠５周年
08/05/13 19:09:02
立方体って定義できますか？
919:名無しさん＠５周年
08/05/13 19:14:59
どう定義しても、それが君にとっての立方体だ。定義というのはそういうもんだろ。
920:名無しさん＠５周年
08/05/14 10:13:47
検知力テスト（2008年5月12日）
Most Effective Antivirus Tools Against New Malware Binaries URLﾘﾝｸ(mtc.sri.com)
Rank Detects Missed Product
1st 95% 77 Ikarus
2nd 93% 133 AntiVir (優)
3rd 92% 148 BitDefender (優)
4th 92% 152 Webwasher-Gateway
5th 90% 185 AVG (優)
6th 88% 212 CAT-QuickHeal
7th 87% 238 F-Secure
8th 87% 241 Norman
9th 86% 248 Kaspersky (良)
10th 83% 311 ClamAV
11th 83% 321 Sophos
12th 83% 322 Microsoft (可)
13th 82% 331 Avast (可)
14th 80% 368 eTrust-Vet
15th 80% 369 VirusBuster
16th 79% 382 TheHacker
17th 79% 387 DrWeb
18th 78% 416 F-Prot
19th 77% 429 AhnLab-V3
20th 75% 468 Symantec (不可)
21st 73% 511 Rising
22nd 72% 529 VBA32
23rd 70% 566 Fortinet
24th 70% 568 Panda
25th 69% 572 McAfee
26th 67% 623 NOD32v2 (Joke!)
921:名無しさん＠５周年
08/05/14 13:08:27
すみません．説明不足でした．
まだはじめたばかりなんで上手く説明できませんが，クラインの壷などのように
立方体（四角柱）を表示させるための数式（？）ってありますか？
922:名無しさん＠５周年
08/05/15 04:26:18
>>921
参考のために、そのクラインの壺を表示させるための数式ってのを示してくれないか。
そうしたら、立方体を示すこともできるんじゃないかな。
923:名無しさん＠５周年
08/05/15 16:02:08
maximaについて調べたときに出てきた「クラインの壷」と書かれていたのがこの式です
plot3d([5*cos(x)*(cos(x/2)*cos(y)+sin(x/2)*sin(2*y)+3.0)-10.0,-5*sin(x)*(cos(x/2)*cos(y)+sin(x/2)*sin(2*y)+3.0),5*(-sin(x/2)*cos(y)+cos(x/2)*sin(2*y))],[x,-%pi,%pi],[y,-%pi,%pi]);
これと同じような方法で立方体を表示させることができるのかという点で今悩んでいます
補足として私はwxMaximaを使用しております
924:名無しさん＠５周年
08/05/15 20:05:24
四角い平面を数式化することもできないの？それずらすだけだよな？
925:名無しさん＠５周年
08/05/16 15:40:20
平面は下のようにしたんですがそこからz軸方向に平行移動させるには
どのように修正すればいいのか・・・・・
plot3d([x,y,0],[x,0,10],[y,0,10]);
上の式がそもそも間違っているかもしれませんが。
926:名無しさん＠５周年
08/05/16 15:50:37
plot3d([x,y,1],[x,0,10],[y,0,10]);
とかすりゃ、z軸方向に移動するだろう。
複数面を同時に描くのは難しいかも。（plot3dの制約）
927:名無しさん＠５周年
08/05/16 16:01:10
ということは平行移動の軌跡みたいなことが表示できない限り立方体を
表示させることは難しい感じですね
ベクトル表記で表示させたりとかは可能でしょうか？
928:名無しさん＠５周年
08/05/17 13:29:56
>>923
URLﾘﾝｸ(aikofan.dee.cc)
ここに3Dグラフをうｐっした。
凄杉
929:名無しさん＠５周年
08/05/21 15:46:58
>>928
URLﾘﾝｸ(maxima.sourceforge.net)
930:名無しさん＠５周年
08/05/21 20:13:00
>>929
(∩´∀｀)∩ﾜｰｲ
見ました。メビウスの輪も書けるとはスゴス
931:名無しさん＠５周年
08/05/22 22:55:42
>>911-912
Mathematica Link for Excelみたいなのが欲しいって事なのかな？
Maxima Link for OpenOffice.org/StarSuite(仮称)なら需要あるのかも。
932:名無しさん＠５周年
08/05/23 01:29:04
Maxima Link for Koffice/gnumeric
933:911
08/05/24 17:26:01
>>931
>Mathematica Link for Excelは、1,000 種類を超える関数やオプションを
>Excel 上で利用できます。
>関数やオプションは、Mathematica ファンクション・ウィザードを使うことで、
>インタラクティブに探し求めることができます。
>Excel 関数の書式から数式処理を実行したり、複素根を求めたり、
>方程式をプロットしたりといった様々な演算を素早く実行できます。
まさにこんな感じです。そんなに多機能でなくても良いので
『Maxima Link for Excel（仮称）』というフリーソフトがあったらいいなと。
>>Excelってとこでやる気減退する人が多いと想像する
なんですかね・・・。
934:名無しさん＠５周年
08/05/24 20:21:03
>>933
maximaのプログラムを書く人たちって、linux/bsd系だろうから、
その時点でexcelとのリンクは頭にないかもしれないね。
gnumericあたりならpythonでプラグインが作れたりするから、
作ろうとする人はいるかもね。swankサーバーからmaxima が
使えるようになるといいとは思うけどな。
935:名無しさん＠５周年
08/05/25 21:01:58
L d/dt I(t)+ R I(t)+(1/C) q(t)=0
をmaximaで解けますか?
936:名無しさん＠５周年
08/05/25 22:02:53
(%i9) set_display('none)$
(%i10) df:L*'diff(I,t,2)+R*'diff(I,t,1)+1/C=0;
(%o10) 'diff(I,t,1)*R+'diff(I,t,2)*L+1/C = 0
(%i11) ode2(df,I,t);
Is R zero or nonzero?nonzero;
(%o11) I = %k2*%e^-(t*R/L)-(t*R-L)/(C*R^2)+%k1(%i12)
937:名無しさん＠５周年
08/05/25 22:26:58
>>936 (∩´∀｀)∩ﾜｰｲ
もうとけたの。すごーい。
3種の解が求まるはずなんですが、これは
特性方程式の判別式D>0の解でつね。
938:名無しさん＠５周年
08/05/25 22:54:11
>>937
あ、間違えてたｗ
(%i12) df:L*'diff(I,t,2)+R*'diff(I,t,1)+1/C*I=0;
(%o12) 'diff(I,t,1)*R+'diff(I,t,2)*L+I/C = 0
(%i13) ode2(df,I,t);Is C*(C*R^2-4*L) positive, negative, or zero?p;
(%o13) I = %k1*%e^(t*(sqrt(R^2/L^2-4/(C*L))-R/L)/2) +%k2*%e^(t*(-sqrt(R^2/L^2-4/(C*L))-R/L)/2)
positive, negative, or zero?
の答え方で３つに分かれます。
939:名無しさん＠５周年
08/05/27 06:07:15
>>938
ｔｈｘ
940:名無しさん＠５周年
08/06/02 16:11:09
>>931
common lispのpyffiにOppen Ofiiceとの接続デモがある
excelじゃないけど、こっちじゃだめなの？
941:名無しさん＠５周年
08/06/08 01:21:22
clispからＯＬＥ呼び出せれるんなら、簡単につくれそうだけど。
もちろん、私は、使わないから作らない。
942:名無しさん＠５周年
08/06/14 08:15:21
論理式は扱えないですか？
943:名無しさん＠５周年
08/06/21 09:52:41
URLﾘﾝｸ(code.google.com)