MAIN { Main data { DataSet name `` file `Int.` rows d200 cols d3 col { Column type d1 N d200 nr d1 nc d1 min f0 max f6.28000021 } col { Column type d2 N d200 nr d1 nc d1 min f-1.623664021 max f0.376297593 eqn `(1- c1*cos(ikep(c1,x)))*(cos(ikep(c1,x))- c1)/ (1- c1*cos(ikep(c1,x)))` } col { Column type d2 N d200 nr d1 nc d1 min f-0.781658769 max f0.781661808 eqn `(1- c1*cos(ikep(c1,x)))*sqrt(1- ((cos(ikep(c1,x))- c1)/(1- c1*cos(ikep(c1,x))))^2)* (if(0,ikep(c1,x),pi)- .5)*2.0` } } data { DataSet name `` file `Int.` rows d1 cols d8 col { Column type d0 N d1 nr d1 nc d1 min f5 max f5 data l1 5 } col { Column type d2 N d1 nr d1 nc d1 min f-1.4214077 max f-1.4214077 eqn `(1- c1*cos(ikep(c1,t)))*(cos(ikep(c1,t))- c1)/ (1- c1*cos(ikep(c1,t)))- .02` } col { Column type d2 N d1 nr d1 nc d1 min f0.471375257 max f0.471375257 eqn `(1- c1*cos(ikep(c1,t)))*sqrt(1- ((cos(ikep(c1,t))- c1)/(1- c1*cos(ikep(c1,t))))^2)* ifs(0,drem(ikep(c1,t),2*pi),pi)- .02` } col { Column type d2 N d1 nr d1 nc d1 min f-1.381407738 max f-1.381407738 eqn `(1- c1*cos(ikep(c1,t)))*(cos(ikep(c1,t))- c1)/ (1- c1*cos(ikep(c1,t)))+ .02` } col { Column type d2 N d1 nr d1 nc d1 min f0.511375248 max f0.511375248 eqn `(1- c1*cos(ikep(c1,t)))*sqrt(1- ((cos(ikep(c1,t))- c1)/(1- c1*cos(ikep(c1,t))))^2)* ifs(0,drem(ikep(c1,t),2*pi),pi)+ .02` } col { Column type d1 N d1 nr d1 nc d1 min f0 max f0 } col { Column type d1 N d1 nr d1 nc d1 min f0 max f0 } col { Column type d1 N d1 nr d1 nc d1 min f0 max f0 } } data { DataSet name `` file `Int.` rows d1 cols d8 col { Column type d0 N d1 nr d1 nc d1 min f5 max f5 data l1 5 } col { Column type d0 N d1 nr d1 nc d1 min f-0.025 max f-0.025 data l1 -0.025 } col { Column type d0 N d1 nr d1 nc d1 min f-0.025 max f-0.025 data l1 -0.025 } col { Column type d0 N d1 nr d1 nc d1 min f0.025 max f0.025 data l1 0.025 } col { Column type d0 N d1 nr d1 nc d1 min f0.025 max f0.025 data l1 0.025 } col { Column type d0 N d1 nr d1 nc d1 min f0.600000024 max f0.600000024 data l1 0.600000024 } col { Column type d0 N d1 nr d1 nc d1 min f0.699999988 max f0.699999988 data l1 0.699999988 } col { Column type d0 N d1 nr d1 nc d1 min f0 max f0 data l1 0 } } data { DataSet name `` file `Int.` rows d2 cols d8 col { Column type d0 N d2 nr d1 nc d1 min f1 max f1 data l2 1 1 } col { Column type d1 N d2 nr d1 nc d1 min f0 max f0 } col { Column type d1 N d2 nr d1 nc d1 min f0 max f0 } col { Column type d2 N d2 nr d1 nc d1 min f-1.496039033 max f-1.278679132 eqn `(1- c1*cos(ikep(c1,t+ c2*ifs(row,0))))*(cos(ikep(c1,t+ c2*ifs(row,0)))- c1)/(1- c1*cos(ikep(c1,t+ c2*ifs(row,0))))` } col { Column type d2 N d2 nr d1 nc d1 min f0.382158816 max f0.590662062 eqn `(1- c1*cos(ikep(c1,t+ c2*ifs(row,0))))*sqrt(1- ((cos(ikep(c1,t+ c2*ifs(row,0)))- c1)/(1- c1*cos(ikep(c1,t+ c2*ifs(row,0)))))^2)* ifs(0,drem(ikep(c1,t+ c2*ifs(row,0)),2*pi),pi)` } col { Column type d1 N d2 nr d1 nc d1 min f0 max f0 } col { Column type d1 N d2 nr d1 nc d1 min f0 max f0 } col { Column type d1 N d2 nr d1 nc d1 min f0 max f0 } } plot { lsPlot name `` pfl d0 col C0.000000 0.000000 0.000000 -1.000000 c0 d1 c1 d2 c2 d3 c3 d2 flags d1 linw f1 lind d0 sym d0 syms f5 } pds d0 plot { shapePlot name `` pfl d0 flags d0 marg f0.050000001 linw f1 aro1 f30 aro2 f60 arol f8 } pds d1 plot { shapePlot name `` pfl d0 flags d0 marg f0.050000001 linw f1 aro1 f30 aro2 f60 arol f8 } pds d2 plot { shapePlot name `` pfl d0 flags d1 marg f0 linw f1 aro1 f30 aro2 f60 arol f8 } pds d3 view { PlotView2D axc C0.000000 0.000000 0.000000 -1.000000 xmac C0.000000 0.000000 0.000000 -1.000000 xmic C0.333323 0.333323 0.333323 -1.000000 ymac C0.000000 0.000000 0.000000 -1.000000 ymic C0.333323 0.333323 0.333323 -1.000000 xlab `` ylab `` linew f1 xspc f0.5 yspc f0.5 xmin d0 ymin d0 flags d2 xfl d0 yfl d0 xmal f4 ymal f4 xmil f2 ymil f2 xhil f0 yhil f0 xlol f0 ylol f0 xdig d0 ydig d0 mwid f5 mhig f5 lfont F14.000000 Helvetica sfont F10.000000 Helvetica legox f50 legoy f50 dlim R-2.000000 -1.500000 3.000000 3.000000 drec R5.000000 5.000000 376.000000 321.000000 dwin `292 269 388 363 0 0 1084 864 ` } vpl d0 vpl d1 vpl d3 vpl d2 var l3 0 6.28000021 2.069792986 var l3 0 1 0.623702407 var l3 0 1 0.256055355 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 var l3 0 1 0 dtool { DataTool mode d1 dsel0 d0 dsel1 d0 dsel2 d5 vsel d0 } winvm `1 292 631 533 133 0 0 1084 864 ` winwc `1 738 1 347 346 0 0 1084 864 ` winws `0 251 175 514 135 0 0 1084 864 ` winwm `0 644 1 333 332 0 0 1084 864 ` winin `1 769 1 310 391 0 0 1084 864 ` winps `1 292 1 447 269 0 0 1084 864 ` winnv `1 4 89 289 479 0 0 1084 864 ` winfit `0 3 239 226 300 0 0 1084 864 ` } {\rtf0\ansi{\fonttbl\f0\fswiss Helvetica;} \margl40 \margr40 \f0\b0\i0\ulnone\ql\fs24\fi0\li0\gray0\fc0\cf0\up0\dn0 This is a plot of a keplerian orbit, solved numerically in the time domain. The arrows represent fixed time offsets from the position of the 'planet'.\ \ Use the 't' variable to vary the orbital time.\ \ c1 varies the ellipticity of the orbit.\ \ c2 varies the separation of the arrows in time. }