Model the orbits of Jovian planets, using the same simple symplectic-integrator.
Thanks to Mark C. Lewis for suggesting this task.
Useful symplectic integrators are freely available, for example the HNBody Symplectic Integration Package.
How to implement
We ask that contributed programs not only give the correct result, but also use the same algorithm to calculate that result.
Each program should:
ndiff -abserr 1.0e-8 program output N = 1000 with this output file to check your program output has the correct format, before you contribute your program.
Use a larger command line argument (50000000) to check program performance.