Now simplified program code available as starting point for your own tries.
Pascal sources and executables running under Linux and DOS.
* 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
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.