If you are satisfied with a rough estimate to highlight key components, consider creating a lookup table for all planets over several years. For example, calculating the orbital period of Neptune which is 164 years can allow you to determine the positions of each planet every month within that timeframe to compile a manageable table size. To visualize the changes in orbital speed more accurately, finer resolution is necessary. Once the calculations are complete, constructing an animation to showcase the planetary positions is the next step.
The computations involved in this process are complex and detailed. Rather than going through each calculation here due to its length, you can find a comprehensive explanation here, along with a sample program coded in QBasic.
The fundamental steps include:
1. Determining the position of the planet in its orbit
2. Calculating the number of days since the elements' date
3. Finding the mean anomaly based on the Mean Longitude and daily motion
4. Obtaining the true anomaly using the Equation of Centre
5. Establishing the radius vector of the planet
After referencing this position to the Ecliptic, ascertain the heliocentric ecliptic coordinates of the planet.
Once you have the heliocentric coordinates, convert them to your specific frame of reference. The provided page demonstrates this transformation for geocentric coordinates, but you will need to adapt it accordingly for your needs. Incorporate these transformed coordinates into your table.
You may opt to run the calculations in real-time, offering more flexibility but potentially impacting the frame rate. It might be beneficial to conduct some trial and error experiments in this regard.
A special thanks to Keith Burnett (the author of the linked page) for providing the detailed information summarized above.