maximaで領域の塗りつぶしをしたい。

load(implicit_plot)
implicit_plot([x^2+y^2=4,x+y=1],[x,-5,5],[y,-5,5])

とするとグラフAが得られますが、これをグラフBのように色分けしたいのです。

（１）plot2dとかimplicit_plotを使う場合
・maximaのplot2d、つまり陽関数表現で出力する場合には三項演算子を使って記述するテクニックが各所に掲載されています。　http://detail.chiebukuro.yahoo.co.jp/qa/question …　http://stackoverflow.com/questions/38589530/fill …　http://stackoverflow.com/questions/13394080/gnup …
・implicit_plotを用いた陰関数表現の場合にはついては記載がみつかりません。

（２）draw2dを使う場合
draw2d出力についても調べましたが、http://www.johnlapeyre.com/qinf/qinf_html/maxima …　によると、
"When filled_func is true, the region bounded by the function defined with object 'explicit' and the bottom of the graphic window is filled with fill_color."
であり、
○draw2d(filled_func=true,explicit(x^2,x,-5,5),explicit(x-1,x,-5,5));
×draw2d(filled_func=true,implicit(x^2+y^2=4,x,-5,5,y,-5,5),implicit(x+y=1,x,-5,5,y,-5,5));
要するに陽関数表現の時のみfilled_funcは有効になるようです。

陰関数表現されたグラフを色分けする方法ってないんでしょうか？
maxima自体のマクロとか、外部プログラムでも実現できれば助かるのですが。。

このあたりについて詳しい方がおられましたらよろしくお願いします。

No.5 ベストアンサー 回答者： jurupapa 回答日時： 2016/11/23 02:35

draw2d(x_voxel=100,y_voxel=100,region((x^2+y^2-4)*(x+y-1)<=0,x,-4,4,y,-4,4),implicit(x^2+y^2=4,x,-4,4,y,-4,4),explicit(-x+1,x,-4,4));

詳しくはMaximaマニュアル52章drawの中のグラフィックオブジェクトregionの定義を参照してください。
http://maxima.osdn.jp/maxima_52.html#IDX1891

この回答へのお礼

素晴らしい。感謝します！

同時に、maximaについて使いこなせていない部分が多々あると痛感いたしました。

お礼日時：2016/11/23 10:44

No.4 回答者： GBde 回答日時： 2016/11/22 19:59

頭の中で　(例えば　負)領域D の塗り絵 を　しない日はないでしょう.

　
http://userdisk.webry.biglobe.ne.jp/020/691/47/N …

12の倍返しの2*12次　代数曲線の　負領域に　塗り絵をしました。

　FAQ で　塗り絵したら　其の部分の面積を　求めずにはイラレナイ筈。
　　
　　　　　 どうぞ!

No.3 回答者： GBde 回答日時： 2016/11/22 18:30

https://people.richland.edu/james/spring15/m122/ …
に　塗り絵　在り。

この回答へのお礼

Shading a Region のdraw2dのオプションはexplicitなんですよ～。
implicitでできたらいいんだけども...

お礼日時：2016/11/22 19:05

No.2 回答者： GBde 回答日時： 2016/11/22 12:17

http://www.wolfram.com/mathematica/new-in-11/?sr …

この回答へのお礼

クラウド上のマセマティカでも描画できました！

http://www.wolframalpha.com/input/?i=RegionPlot% …

ホームユースでも年間1.5万円前後かかる(？)ってのも微妙～
RegionPlotについては、クラウドかラズベリーパイ導入するかにします。。

お礼日時：2016/11/22 19:17

No.1 回答者： GBde 回答日時： 2016/11/22 12:12

mathematica なら　↓で　瞬時に叶います;

RegionPlot[(x^2 + y^2 - 4)*(x + y - 1) <= 0, {x, -3, 3}, {y, -3, 3}]

この回答へのお礼

mathematicaのグラフィックスは一日の長がありますね！さすがと言わざるを得ません。

だけど、RegionPlotに相等するmaximaのマクロがあったら、新たに記法覚えなくて済むし、ありがたいんだけど...

悩むところです。

お礼日時：2016/11/22 17:44