Octaveやmaximaなどの数学系のフリーソフトのスレ

**１３２人目の素数さん** · 2012/01/07(土) 18:53:36.06

タダより貴重なソフトは無いね

**１３２人目の素数さん** · 2015/06/05(金) 22:39:39.70

maple18だとeq:=x^2+x+1 と表示される

**１３２人目の素数さん** · 2015/06/05(金) 23:59:45.12

Maple11では無理なんですかね。。

**１３２人目の素数さん** · 2015/06/09(火) 18:52:14.59

maple9でもeq:=x^2+x+1と表示された

**１３２人目の素数さん** · 2015/06/23(火) 22:07:59.33

maximaのsolveって複素数解もつ連立方程式は解けない？
solve([x^2+y^2=1,(x-3)^2+(y-5)^2=1],[x,y])
試してみたんだけど{}だった。(x-3)^2+(y-5)^2=25にしたら実数解は出たんだけど

**ご冗談でしょう？名無しさん** · 2015/07/13(月) 11:05:06.20

Maple11の環境です。

以下の行列の展開式でxyzの項だけ取り出したいのですがcollectコマンドが機能しません。
何が原因なのでしょうか?

with(LinearAlgebra);

A := Matrix(3, 3, {(1, 1) = a11, (1, 2) = conjugate(a21), (1, 3) = conjugate(a31), (2, 1) = a21, (2, 2) = a22, (2, 3) = conjugate(a32), (3, 1) = a31, (3, 2) = a32, (3, 3) = a33});
B := Matrix(3, 3, {(1, 1) = b11, (1, 2) = conjugate(b21), (1, 3) = conjugate(b31), (2, 1) = b21, (2, 2) = b22, (2, 3) = conjugate(b32), (3, 1) = b31, (3, 2) = b32, (3, 3) = b33});
Q := Matrix(3, 3, {(1, 1) = q11, (1, 2) = conjugate(q21), (1, 3) = conjugate(q31), (2, 1) = q21, (2, 2) = q22, (2, 3) = conjugate(q32), (3, 1) = q31, (3, 2) = q32, (3, 3) = q33});

eq0 := Determinant(x*A+y*B+z*Q); eq1 := collect(eq0, x); eq2 := collect(eq1, y); eq3 := collect(eq2, z);

**１３２人目の素数さん** · 2015/07/14(火) 02:17:58.43

Maple とか使わないから知らんけど
coeff を使うとか？

**ご冗談でしょう？名無しさん** · 2015/07/14(火) 06:05:02.60

お蔭様で上手くいきました。

**１３２人目の素数さん** · 2015/07/23(木) 22:29:53.37

octave触り始めたんですが、
mac版は3dのプロットがすごくカクカクなんですが
ほかもそうなんでしょうか？
4.0で劇的に変わったりします？まだmac版が無いっぽいですが。

**１３２人目の素数さん** · 2015/10/13(火) 14:37:03.79

octave 4.0 で2G以上のメモリを確保するにはどうすればいいですか？

**１３２人目の素数さん** · 2015/10/15(木) 07:17:33.91

>>185
windowsを捨ててmacかLinuxを使う。

**１３２人目の素数さん** · 2015/10/28(水) 20:25:38.03

MaximaってIEのjavascriptからライブラリとして呼び出せませんか。

**１３２人目の素数さん** · 2015/10/31(土) 00:38:15.70

点(x1,y1),(x2,y2)を通る直線を求めたいときsolveをつかって何とかうまいこと求められないでしょうか。
p: [x=x1,y=y1]
q: [x=x2,y=y2]
l: solve([a*x+b*y+c=0,p,q][a,b,c])
のような感じでa,b,cが求められるとよいのですが。

**１３２人目の素数さん** · 2015/10/31(土) 20:01:22.40

p: [x=x1,y=y1];
q: [x=x2,y=y2];

line: a*x+b*y+c=0;

l1: line,p;
l2: line,q;

solve([l1,l2],[a,b,c]);

こうかな

**１３２人目の素数さん** · 2015/10/31(土) 22:00:54.42

line: a*x+b*y+c=0;

get_line(p,q) := sublis(
append(
solve([sublis(p,line),sublis(q,line),c=1],[a,b,c]),
solve([sublis(p,line),sublis(q,line),a=1,c=0],[a,b,c]))[1],
line);

できた。

**１３２人目の素数さん** · 2015/10/31(土) 22:09:16.24

p=qの時の処理も入れなきゃな。

**１３２人目の素数さん** · 2015/10/31(土) 22:55:21.72

y=0も

**１３２人目の素数さん** · 2015/11/01(日) 00:16:16.25

座標 (a′, b′)の点と座標 (c′, d′)の点を通る直線は
(d′ ? b′)(x ? a′) + (c′ ? a′)(y ? b′) = 0

wikiより

**１３２人目の素数さん** · 2015/11/02(月) 16:34:07.48

octave でfopenするときって、ファイルはカレントディレクトリにあればいいんですか？

**１３２人目の素数さん** · 2015/11/02(月) 21:25:02.77

maximaで円の連立方程式を解かせようとしたら、
以下のようになりました。

(%i68) solve([(x-2)^2+(y-2)^2=8,x^2+y^2=8],[x,y]);

(%o68) []

なぜ解が求まらないのでしょうか。

**１３２人目の素数さん** · 2015/11/08(日) 21:43:09.21

Windows7でMaximaをいれたところsin(x)の微分や式表示ですら
バグって応答なしになってしまいます
どうすればいいですか

**１３２人目の素数さん** · 2015/11/08(日) 22:45:21.53

バージョン落としてインストールしてみて

**１３２人目の素数さん** · 2015/12/23(水) 10:18:03.67

>>195
Maximaのバグらしいから悩むのはやめよう
ttp://sourceforge.net/p/maxima/bugs/2736/

**１３２人目の素数さん** · 2016/01/03(日) 02:41:09.75

maximaをバージョンアップした時によくバグが入っていることが多い。
使わさせてもらってる立場であるから言い難いけど、
どういう管理をしてるのかな？ってちょっと思う。

**１３２人目の素数さん** · 2016/01/19(火) 23:23:37.33

Mathematicaとの格差ががが

**１３２人目の素数さん** · 2016/01/21(木) 21:33:46.68

octave 4.0 windows版って2GB以上のデータ使えるようになった？

**１３２人目の素数さん** · 2016/01/23(土) 13:13:20.83

Mathematicaもバグ多いけどな。

**１３２人目の素数さん** · 2016/01/23(土) 16:55:16.81

(ﾟДﾟ)ｴｯ?

**１３２人目の素数さん** · 2016/01/23(土) 21:17:06.10

今のmathematicaは至れり尽くせり。高機能になるとバグも増えるのは分かるがそれでもmaximaよりは
ずっと少ない。mathematicaで連立方程式が解けないとか間違った解を出すなんてのはありえないよ

**１３２人目の素数さん** · 2016/01/24(日) 15:23:16.85

maximaの結果が正しいかどうか検算って虚しいな

**１３２人目の素数さん** · 2016/01/24(日) 15:57:05.84

結果の保証は誰もしてくれない。

**１３２人目の素数さん** · 2016/01/25(月) 17:45:55.19

Mathematicaも、使い込むと色々と出てくるわけよ。
去年はバグ報告の回数が２桁になったw

**１３２人目の素数さん** · 2016/01/25(月) 20:16:20.22

Mathematicaは更新がはやい。10.3.1じゃとっくに解決してる
商用だから当たり前だが
いくらフリーだからって>>198レベルのを放置は管理する気ないだろ
フリーでも火狐やUbuntuはしょっちゅう更新するってのに

**１３２人目の素数さん** · 2016/01/26(火) 00:23:11.57

mathematicaの話はここでしなくていいよ。
Expand[x*(x-1),x]とかEと大文字にしてしかもsift+returnが
うざいし、UIの外観を変えるためのプログラミングもうざい。
mathematicaはlispに似てると聞いたけどどうなんだろ？
texへの文書の変換もややこしくない？>mathematicaの人

maximaはlispを使って動いてますよね。emacsをmaximaのUIに
使うこともできるんでしょうね。

**１３２人目の素数さん** · 2016/02/15(月) 03:05:33.74

maximaで括弧が表示されなくなったのですが解決方法が分かりません

例えばX^2+2*X+1　の因数分解をさせると、
(x+1)^2ではなく、X+1^2と表示されてしまいます。

**１３２人目の素数さん** · 2016/02/15(月) 07:45:31.21

最新ですか

**210** · 2016/02/15(月) 17:06:06.16

>>211
バージョンは最新みたいです

**１３２人目の素数さん** · 2016/02/21(日) 11:48:43.23

>>210-211
Windowsの5.36.1(SBCL)版のコマンドラインではちゃんと動くっぽい。

(%i1) display2d:false;
(%o1) false
(%i2) factor(X^2+2*X+1);
(%o2) (X+1)^2

実は表示だけの問題だったら(この例ならexpand()で元に戻るのなら)、
気にしないのも手だなｗ

>>210
現象の切り分けは出来てますか？
コマンドライン版Maxima起動直後でも再現するか？
それとも特定のGUIでのみ再現するのか？

**210** · 2016/02/24(水) 16:48:35.06

>>213
返信遅くなって申し訳ないです。

コマンドライン版では正しく表示されるようですが、
GUI版だと、アスキー形式または形式なしで出力しない限り表示できないらしく、
xml形式での表示ができない状態でした。

xmlに問題があるのでしょうか。

**213** · 2016/03/04(金) 12:00:13.14

これですよね？
ttp://maxima.osdn.jp/maxima_39.html

>>214
GUIもいろいろあるぞ(XMaxima、wxMaxima、TeXemacs…)。

xmlじゃないと絶対に駄目ではなくって、アスキー/形式なしでも支障ないのなら
今は諦めてGUIを作ってる人の対応待ちかなぁ。

**210** · 2016/03/06(日) 17:55:09.82

>>215
ありがとうございます。

慣れれば別段問題なさそうなので、、
とりあえずアスキー形式で使おうと思います。

**１３２人目の素数さん** · 2016/03/08(火) 12:50:04.99

>>195

(%i1) solve([x^2+y^2=8,(x-2)^2+(y-2)^2-(x^2+y^2)],[x,y]);
(%o1) [[x=1-sqrt(3),y=sqrt(3)+1],[x=sqrt(3)+1,y=1-sqr(3)]]

**１３２人目の素数さん** · 2016/03/29(火) 23:12:57.08

octaveをhomebrewでmacにインストールできないんだよな
なんでなんだろ
el capitanでbrew でインストールしてるときにエラーがでる

**１３２人目の素数さん** · 2016/03/30(水) 23:26:28.23

ウルフラムアルファすげくね？
スレ違い？

**１３２人目の素数さん** · 2016/04/03(日) 21:51:08.90

>>219
Mathematicaと比べて、微妙に機能制限してるのは…
商売上、仕方がないんでしょうね。
RealDigits[]使おうとしたら使えなくてちょっと困ったｗ

**220** · 2016/04/03(日) 22:48:30.63

「ラズパイ買う金も無いのかよｗｗｗｗ」と云われると…その通りなんだけど
orz

**１３２人目の素数さん** · 2016/04/08(金) 19:26:54.83

maximaについて質問です。ver.11.08.0です。
f(x):=diff(g(x),x);
subst(g(x)=X(x)+Y(x),f(x))-diff(X(x),x);
と入力すると、出力は
d/dx(X(x)+Y(x))-d/dx(X(x))
となるんだけど、この出力を得た上で更に答えを
d/dx(Y(x))
に持っていくにはどうしたらいいですか？

**１３２人目の素数さん** · 2016/04/10(日) 15:34:32.46

>>222
現行のMaximaは5.38.0が最新みたいだけど

**１３２人目の素数さん** · 2016/04/10(日) 16:34:14.23

>>222
積分して微分したらできた。とりあえず

integrate(%, x);
diff(%,x,1);

wxMaxima-15.08.2 maxima-5.37.2

**222** · 2016/04/11(月) 10:17:52.96

>223
wxMaximaだった。。無知でした
>224
ありがとう。これで当面の目的は達せられそう

**１３２人目の素数さん** · 2016/04/22(金) 22:17:42.21

ウルフラムアルファまじすげぇええええぇ

Input interpretation:

(1+2+3+...+n)^2=1^3+2^3+3^3+...+n^3

Result:

True

**１３２人目の素数さん** · 2016/05/11(水) 14:59:10.17

wxMaximaでsin^2(M*Δk/2)/sin^2(Δk/2)を-3<Δk<3で描きたいです。Mは後で自由に変えて、例えばM=50,100を同じグラフ上に描きたいです。

**１３２人目の素数さん** · 2016/05/11(水) 18:02:15.67

>>227
sin(M*x/2)^2/sin(x/2)^2;

wxplot2d([at(%, M=50), at(%, M=100)], [x, -3, 3]);

**１３２人目の素数さん** · 2016/05/29(日) 08:50:11.06

64ビット版のoctaveが使いたいなら、私がビルドしたものをあげてある。
公式チームに64ビット版問題なく使えるからあげてくれと申し入れているがあげてくれないから、自分のサイトにあげた
http://www.tatsuromatsuoka.com/octave/Eng/Win/

ちなみに私は、ここはじめてきたし、これからも来ないとおもうので
なにかあったらunix板のoctaveスレに来てくれ。

**１３２人目の素数さん** · 2016/05/29(日) 10:13:14.99

>>208
SourceForge でbug ticketをcreateしてそこに書いたら？
Maximaづかいではないから（インスコはしてある）口ださんけど
常用している gnuplot，octaveはbugを見つけたらさっさと本家に挙げている。
2chで話題にするのは他の人バグをしらせるという意味はあるからいいけど，
もう少し踏み込んで本家にあげてほしい。
フリーウェアではユーザーからのフィードバックは大事だし，それがそのソフト
の発展につながるんだから。

**１３２人目の素数さん** · 2016/05/29(日) 10:22:34.99

>>218
私はMacづかいではないが，本家のhelp MLにいきな。
その手の質問はしょっちゅう出てくるし，親切なMacのoctaveユーザーが答えてくれる。

**１３２人目の素数さん** · 2016/05/30(月) 07:20:39.35

229-231です。もう来ないようなこと書いてすみませんm(__)m。書き方も横柄ですみません。
で，octaveは数学アプリケーションというより数値計算アプリケーションだからあるゆる科学，工学，経済学，統計学などをカバーします。 Maximaはoctaveよりはるかに数学色が強いですが，私のような数学が苦手な工学系の人間にとって便利なツールです。
数学板にあることは別に問題はありませんが，広い幅の人々を受け入れてくれるとよいなあと思っています。

**１３２人目の素数さん** · 2016/05/30(月) 15:35:59.42

マツオカくん
ここはユーザーメーリングリストじゃ
ないんだから突っぱねるような
レスはよしなさいw

**１３２人目の素数さん** · 2016/05/30(月) 17:59:28.35

>>233 だから、反省して232で謝ってんじゃないですか。本当にm(__)m。
まあ、仲良くしていきましょう。お願いします。m(__)m。

Mac上のoctaveのインストールに関しては以下を見てとなっているんですが、

http://wiki.octave.org/Octave_for_MacOS_X

結構難しそうですね。
Macでのoctaveでわからないことがあれば、やっぱり頑張ってOctaveのHelp MLにいくのが一番です。
親切なMacのパワーユーザーが教えてくれるます。そういう事例はなんどもみましたから。

**１３２人目の素数さん** · 2016/05/30(月) 18:23:02.17

>>234

octaveのMacへのインストールは、以下が参考になるかどうか。

http://octave.1599824.n4.nabble.com/octave-on-Mac-td4677113.html

**１３２人目の素数さん** · 2016/05/30(月) 19:36:06.22

今さらですが
>>147

3.6.4のMinGW版はインストーラを使わない。MSVC版もバージョンの違うoctaveとは共存できるはずなのだが。

3.6.4のMinGW版なら、フォルダを丸ごと削除していい。レジストリを使ってないから。
3.6.4のMSVC版なら、octaveのインストールフォルダにuninst.exeがあるから、これを実行すればいい。

でも、怒られるのは他の原因のような気がしますが、思い違いかもしれません。

**１３２人目の素数さん** · 2016/05/30(月) 19:39:19.85

>>194

C言語と同じです。カレントディレクトリにあれば、ファイル名だけでいいです。フルパス、相対パスがつかえます。

**１３２人目の素数さん** · 2016/05/30(月) 19:47:19.81

>>184
3Dのプロットでカクカクなのは、meshgridの切り方が荒いせいでしょう。バージョンはよほど古いものでなければ関係ありません。
ただし、あまり細かく切りすぎると描画が遅くなります。

以下の例は、いかがでしょう。

X=linspace(-pi/2,pi/2,50);
[XX, YY] = meshgrid (X);
ZZ=sin(XX.*XX+YY.*YY);
colormap(bone);
surfl(XX,YY,ZZ);
shading interp;

colormap(bone)は淡色のやや青みがかった灰色のカラーマップに設定する命令であり，shading interpは，メッシュ間の色を補完し滑らかに表現する命令です。

**１３２人目の素数さん** · 2016/05/30(月) 20:06:58.22

>>135
octaveが自前で搭載した表示機能ですが、現在はqtとfltkがあります。

octaveのguiとoctave --no-guiで起動した場合qtとfltKが、octave-cliで起動した場合は、fltkのみ使えます。

普通はqtを使うとおもいます。

どちらも、open GLを使っているので3dのプロットのレンダリングが高速です。

ただし、単精度実数なので、広いレンジの対数プロットや大きなオフセットがある場合のプロットでは問題があります。

その場合はgnuplot graphics toolkitを使う必要があります。

また、TeXモードを使うときもgnuplotを使う必要があります。

場合によって使い分けるのがよいとおもいます。

**１３２人目の素数さん** · 2016/05/31(火) 09:07:24.43

全レス少年

**１３２人目の素数さん** · 2016/05/31(火) 18:38:23.22

>>240 それは違う。全レスオヤジだ。でも、ちゃんと仕事してるぞ。

**１３２人目の素数さん** · 2016/06/01(水) 14:48:39.37

maximaで、方程式の解をグラフ化することはできませんか？

たとえば、x = a cos x の解をaを変えながら求めて、その解をaの関数としてグラフにする、という方法です。

解を求めるだけならfind_rootでできますが。

￥ ◆2VB8wsVUoo · 2016/06/01(水) 15:22:12.96