Now simplified program code available as starting point for your own tries. Pascal sources and executables running under Linux and DOS.
Linux programs * DOS program * Windows32 program
The program runs on all machinae with 80 bit numeric processor. The compiler must support these 80 bits (`Turbo Pascal' and `Free Pascal' do it). Don't reduce to 64 bits.
The Linux package contains a program for 128 bit machinae.
As example, and for comparison, see already achieved results as simple diagrams exclusively in ASCII text format. There are also the reference values as well from literature as the equivalent integration constants to insert into the initial conditions. See also the listed results (the same in German).
The results consist in coincidences of the inserted integration constants at most stable solutions (regarding the chaotical development during the computation) with the known particle values. The program has been designed to test many values of the integration constants in order to find any dependences of the chaotical development on the integration constants. As well, the program counts the computation steps from one radius to the next (with the difference dr). More information is stored by the old DOS versions attached to the works.
The old results should be reproducible on an i486. Unfortunately, the new fast processors have not more the full 80 bit support. That made me begin test series, which make possible statistical evaluation. If you take interest in that, contact me.
The clean solution consists in higher precision (128 bit or more), but the use of high-precision libraries does not more guarantee that the entire algorithm may be error-free. (First tries have shown that my newest machina is overtaxed.) The most realistic solution should be the use of a 128 bit machina.
Please note the copyright.