gnuplot equations

Plot of fX(t), fY(t), fZ(t) where,

a = mean radius (AU)

e = eccentricity

Om = longitude of ascending node (deg)

om = argument of periapse (deg)

i = angle of inclination (deg)

b = a*(1-e*e) = angular momentum

ex(t) = b*cos(2*t)/(1+e*cos(2*t))

ey(t) = b*sin(2*t)/(1+e*cos(2*t))

r1(t) = ex(t)*cos(om)-ey(t)*sin(om)

r2(t) = ex(t)*sin(om)+ey(t)*cos(om)

fX(t) = r1(t)*cos(Om)-r2(t)*sin(Om)*cos(i)

fY(t) = r1(t)*sin(Om)+r2(t)*cos(Om)*cos(i)

fZ(t) = r2(t)*sin(i)

x, y, z = heliocentric coordinates (AU)