... The Programmer God ...
A simulation universe hypothesis at the Planck scale

I am putting the simulation programs online. These are basic versions but can be easily modified to experiment with, if you get some interesting results, or improvements, please let me know.

Gravitational orbits

This gives a simple 2-body orbit, 1 point rotates around a center mass. After 1 orbit, the calculated data is compared with the simulation data. To continue the orbits, comment out line 234 //i = 65535; and //j = 65535;
The algorithm itself is independent of this orbit, it just rotates points around each other, and so you may assign to each point any (x, y) co-ordinates. The number of points is limited only by cpu power and memory allocation for the arrays, you may have to adjust those settings.
download Gravitational-orbitals-2body.c
download Gravitational-orbitals-nbody.c

Python is too slow so the C version will save the data onto a file to be plotted by this python script. With many points the simulation can take days to run. It is supposed to emulate Planck time operation (to answer the question - could gravity be a Planck scale event) and so no optimisation techniques are used.
download orbital_graph.py

Atomic orbits

This program models the electron as it transtions from the n=1 orbital shell to the n=max_n orbital shell. For reference I have included experimental data for the n=2, 3, 4 shells, the n=1 shell is the most precisely measured. The wavelength of the electron is set to 1 and so the H variable is used to remove the wavelength from the experimental data. This version is very basic, it assumes the barycenter is always at (0, 0) which is not realistic, nevertheless the results are surprisingly close. The python version is just a copy of the C version. If yoou want a high max_n then you may have to increase the array sizes, they are set quite small.
download H-orbital.c
download H-orbital.py