gunplotなどの立体ソフトについて

gunplotを使っているのですが、３Ｄでの惑星軌道が描けません。
よろしければそれに関するコマンドを教えていただけませんか。
また、gunplotのほかにもいい惑星軌道が書けるフリーウェア、あるいは安いソフトがあれば教えていただければ幸いです

No.3 ベストアンサー 回答者： masudaya 回答日時： 2007/10/26 16:32

では具体的に

set parametric
splot [0:4*PI] L/(1+e*cos(u))＊cos(u),L/(1+e*cos(u))＊sin(u),0

ではどうですか？前もってLとeは設定ください.
ちなみに私が使っているGnuplotはVer4.0です.

この回答へのお礼

　ありがとうございます。無事出力することが出来ました。
　ｔはそう使うんですね

お礼日時：2007/10/27 21:09

No.2 回答者： masudaya 回答日時： 2007/10/26 10:49

見づらくてすみません.

[t0:tend]
の区切り文字は”：”です.
”,”ではありません.

もう一度お試しください.

半直弦に関しては,私も良くは知りません.
軌道の方程式も楕円軌道とかで検索した結果です.
すみませんが,ご自身でお調べください.

この回答へのお礼

　ありがとうございました。ですが残念ながら何度計算を試みても無理みたいです。ｔ0の項目がどうしてもダメなようです。
　gunplotは諦めてほかのソフトを探すとしましょうか…

お礼日時：2007/10/26 15:13

No.1 回答者： masudaya 回答日時： 2007/10/24 16:50

表示する関数がパラメータtを用いて

x＝x(t),y＝y(t),z＝z(t)
とかけるときは,

set parametric
splot [t0:tend] x(u),y(u),z(u)

となります.パラメータは,u,vという風に用います.
t0,tendはパラメータtの初期値を最終値です.

たとえば,L は半直弦 (semi-latus rectum) として離心率eとすると惑星の軌道は

x＝x(t)＝L/(1+e*cos(t))＊cos(t)
y＝y(t)＝L/(1+e*cos(t))＊sin(t)
z＝0

となります.

この回答への補足

　すいません。補足をつけたせていただきます。
　[t0;tand]の項目なのですが、gunplotがその数値を受け付けてくれません。
　これはバージョンが古いせいなのでしょうか？
　それと半直弦に関しては、（軌道長半径）＊（１－離心率＊離心率）で入力してもかまわないのでしょうか？

補足日時：2007/10/24 19:53