Plotfile data:
set term pbm small color
set view 0,0,1,1
set ytics -4,1,4
set xtics -2,1,6
set ztics -2,.5,2
set ticslevel 0
set grid
set size 1.000,1.30
set title ' Flight Day 0 Jan , 20 ' ,-10
set label "Phase: " at -1.92,3.7,0
set parametric
set samples 100
set xlabel 'x (AU) Vernal Equinox ->' 0,4
set ylabel 'y (AU)' -3,0
set zlabel 'z' ,0
set xrange [-2:6]
set yrange [-4:4]
set zrange [-2:2]
set nokey
rad = pi/180
set label ' Sun ' at 4.2,2.5,0 left
a_e =
e_e =
i_e =
om_e =
Om_e =
b_e = a_e*(1-e_e*e_e)
om_e = om_e*rad
Om_e = Om_e*rad
i_e = i_e*rad
ex_e(t) = b_e*cos(2*t)/(1+e_e*cos(2*t))
ey_e(t) = b_e*sin(2*t)/(1+e_e*cos(2*t))
r1_e(t) = ex_e(t)*cos(om_e)-ey_e(t)*sin(om_e)
r2_e(t) = ex_e(t)*sin(om_e)+ey_e(t)*cos(om_e)
f_eX(t) = r1_e(t)*cos(Om_e)-r2_e(t)*sin(Om_e)*cos(i_e)
f_eY(t) = r1_e(t)*sin(Om_e)+r2_e(t)*cos(Om_e)*cos(i_e)
f_eZ(t) = r2_e(t)*sin(i_e)
set label ' Earth ' at 4.2,2.8,0 left
a_w =
e_w =
i_w =
om_w =
Om_w =
b_w = a_w*(1-e_w*e_w)
om_w = om_w*rad
Om_w = Om_w*rad
i_w = i_w*rad
ex_w(t) = b_w*cos(2*t)/(1+e_w*cos(2*t))
ey_w(t) = b_w*sin(2*t)/(1+e_w*cos(2*t))
r1_w(t) = ex_w(t)*cos(om_w)-ey_w(t)*sin(om_w)
r2_w(t) = ex_w(t)*sin(om_w)+ey_w(t)*cos(om_w)
f_wX(t) = r1_w(t)*cos(Om_w)-r2_w(t)*sin(Om_w)*cos(i_w)
f_wY(t) = r1_w(t)*sin(Om_w)+r2_w(t)*cos(Om_w)*cos(i_w)
f_wZ(t) = r2_w(t)*sin(i_w)
set label ' Wild 2 ' at 4.2,3.1,0 left
set label ' STARDUST ' at 4.2,3.4,0 left
splot [t=0:pi] 0.03,0.03,0. w p 6 6, -0.03,0.03,0. w p 6 6, 0.03,-0.03,0. w p 6 6, -0.03,-0.03,0. w p 6 6 , 4.2,2.5,0 w p 6 6, , , w p 2 3, f_eX(t),f_eY(t),f_eZ(t) w l 2 1 , 4.2,2.8,0 w p 2 3, , , w p 3 3, f_wX(t),f_wY(t),f_wZ(t) w l 3 1 , 4.2,3.1,0 w p 3 3, , , w p 5 3, '/home/meyertr/public_html/3dplot/stardust.dat' using 4:5:6 w l 5 1 , 4.2,3.4,0 w p 5 3