2010年7月15日木曜日

Maximaとの連携成功!



朝の短い時間でクラインの壺(Klein bottle)ができるとは思わなかった。
しかし、色々試行錯誤した結果、Maximaにより"クラインの壺"(Klein bottle)をつくることができました!!!
凄い感動・・・



これは逆さまにしたもの。

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],['grid,40,40]);

によってMaxima(wxMaximaは使いこなせませんでした・・・)でgnuplotと連携させる事により、こういった複雑な図形をプロットする事ができるとわかり、かなり感動しています。
と言ってもまだトポロジーを理解しているわけではないので、クラインの壺がどのような図形であるかを掴みきれていないのですけど。

しかし、なにはともあれ、プロットということに関してかなり希望が持てるようになった瞬間でした!









連携させないで直接出力させようとすると、こんな感じになります。


set size square
set isosample 128,128
set parametric
set hidden3d
set xyplane 0
set pm3d at s
set contour base
set view 30,60
set cntrparam cubicspline
set cntrparam points 3
set cntrparam levels incremental -1, 0.1, 1
unset tics
x(u,v) = 5*cos(u)*(cos(u/2)*cos(v)+sin(u/2)*sin(2*v)+3.0)-10.0
y(u,v) =-5*sin(u)*(cos(u/2)*cos(v)+sin(u/2)*sin(2*v)+3.0)
z(u,v) = 5*(-sin(u/2)*cos(v)+cos(u/2)*sin(2*v))
splot [-15:15] [-15:15] x(u,v) , y(u,v) ,z(u,v) with lines notitle
という感じになります。

Geomviewに目をつけましたが、インストールがうまくいかなかったので、
今度にしようと思います。
-----------------------------------------------------------------



wxMaximaでの出力に帰宅後成功しました。出力結果はこんな感じです。
まだ、調整が出来ていないのでかなりノイズが多いですがご容赦下さい。

MacPortsも一旦uninstallしてから再インストール。

0 件のコメント: